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CHAPTER I 


INTRODUCTION 

1 . Purpose of the study. 

The analysis of the modal behavior of an elastic structure, 
i.e. the determination of its natural frequencies and correspond- 
ing normal modes is an important part of the actual design of the 
structure. The calculation of normal modes is important for several 
reasons. First, by examining the normal modes obtained and particu- 
larly the positions of the ?\odal lines, it may be possible to tell 
whether flutter is likely or not, and in any case such an examination 
will indicate what types of flutter should be investigated. It may 
also be possible to predict whether there is any likelihood of 
resonant response vibration, due to the proximity of the natural 
frequencies to the forcing frequencies of the power plant, gust 
fields, etc. Further, the normal modes are commonly used for the 
actual flutter and response calculations. Theoretically, any set 
of independent, complete deformation modes can be used as degrees 
of freedom in setting up the flutter equations, fe u t from practical 
considerations it is suitable to select those modes which give an 
accurate flutter speed when the number of chosen modes is small. 

This is more likely to be the case when the modes are related to 
the actual structure, as the normal modes are, than when arbitrary 
modes are chosen. An additional reason for using normal modes is 
that they can latter be compared with the resonance - test modes. 

It is then understandable that the modal analysis of elastic 
systems has been given extensive attention in the literature. The 
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list of the various methods developed for this purpose is much too 
large to be included here, therefore the next paragraph will review 
only the main conventional methods that have been used in modal 
analysis . 

2. Historical Background. 

The first half of this century has seen a large development 
of approximate methods for solving the eigenvalue problem for which 
analytic solutions were impossible. In 1921 Holzer [1] developed 
a step-by-step method for solving the differential equation of 
beams in pure torsion; this method was later extended to the case 
of flexural vibrations by Myklestad [2] in 1944. At about the 
same time Stodola [3] developed a method of iterative integration 
for solving the differential equation, while White [4] and Crout 
[5] focused their interest on finding solutions to the integral 
equations by numerical integration and collocation [4] or by 
collocation with assumed functions [5] . In 1939 Grammel [A] 
developed the Complementary Energy Method which was later intro- 
duced by Westergaard [7] and Reissner [8] . Alternatively, Duncan 
[9, 10] applied Galerkin and Rayleigh-Ritz method to the eigen- 
value problems with special polynomial functions used in the series 
expansions. Various matrix methods were developed in the process: 
matrix iteration [11], relaxation methods [12], etc. 

Some of the aforementioned methods failed to apply to the 
problems of complex elastic structures: For instance, in a Galerkin 

method the comparison functions depending upon the spatial coordinates 
must satisfy the boundary conditions associated with the problem 
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and this choice would be particularly difficult in the case of a 
structure such as an airplane. In any case, all the methods that 
have been mentioned so far suffer from the defects inherent in 
replacing a structure by a disc te mass system. Duncan has shown 
by considering the torsion of a uniform cantilever that the. error 
in the frequency, obtained by means of the discrete mass method is 
inversely proportional to the square of the number of degrees of 
freedom used. To obtain an accuracy of one per cent it is necessary 
to use six times as many degrees of freedom as there are desired 
modes. The accuracy for flexure may be a little better. If an 
accuracy to within two per cent u u be achieved, then at least 
four times as many degrees of freedom as there are desired modes 
should be used for either flexure or torsion. 

As a consequence, in the second half of this century most 
of the effort has been directed towards developing new methods 
in which the structures will no longer have discrete representations 
but will be approximated in a continuous manner. Hurty [13, 14] 
has developed a component modes approach for complex structures 
through the use of various types of modes: "rigid body modes", 

"constraint mod e s" and "normal modes", the latter being usually those 
of some (artificially) restrained component. Gladwell [15] speaks 
of ,T b ranch modes" which are again those of an artificially restrained 
component. Wittrick and Williams [37] present an algoritum for 

o 

computing the natural frequencies of elastic structures if their 
dynamic stiffness matrix K(w) corresponding to any finite set of 
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displacements is known. Walton and Steeves [16] present an inter- 
esting alternative analysis which deals directly with the equations 
of motion and constraint. Buckling problems, which are of course, 
mathematically analogous to free vibration problems, have been in- 
vestigated by Budiansky [17] and Reissner [18] . Most recently, 
Dowell has treated the free vibrations of a linear structure with 
arbitrary support conditions by a Rayleigh-Ritz method introducing 
the use of Lagrange Multipliers [19], and he applied a similar 
approach to the free vibrations of an arbitrary structure in 
terms of component modes [20] . The theoretical analysis presented 
in this study is developed from the works of Dowell aforementioned. 

3. Statement of the Problems. 

The intention of the present work is to point out the general- 
ity of the Rayleigh-Ritz Component Modes approach using Lagrange 
Multipliers, by illustrating this method for different types of 
problems . 

In some cases, the problems treated have already been con- 
sidered in the literature by ad-hoc methods, and their choice 
was intentional in order to check the method for accuracy. This 
is the case of the problems treated in Chapters II and V; in other 
cases the theoretical results obtained are checked with experimental 
results (e.g. rotor blade problems of Chapter II, beam-like con- 
figuration of Chapter III). 

Finally, as the method is being proved accurate, new prob- 
lems are investigated, as the oblique wing configuration of 
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Chapter IV. 

A more detailed statement of the problems is now discussed. 

The free vibrations of beams with non uniform properties are 
the object of Chapter II; in addition to the conventional dis- 
crete mass methods already mentioned, almost all the methods de- 
veloped for this type of problems have been ad-hoc methods. For 
instance, Cranch and Adler [23] offer an analytic solution for some 
classes of cantilever beams with rectangular cross section by 
means of Bessel functions, Conway and Dubil [24] used a similar 
approach for the study of truncated cone and wedge beams, and 
many other examples can be found through the literature. The 
theory presented in Chapter II is valid as well for an arbitrary 
variation in the beam properties along the span as for an arbitrary 
set of constraints, or boundary conditions, applying at the edges 
of the beam. For continuous beams with smooth non-uniform charac- 
teristics, an alternative method is developed in Appendix I. It 
is a perturbation - Galerkin method, particularly suitable for 
the determination of harmonics above the fundamental. 

The Rayleigh-Ritz component modes method using Lagrange multi- 
pliers is applied in Chapter III to the free vibrations of a 
flexible airplane with high aspect ratio components. The earliest 
theoretical treatments of this problem have been based upon the 
representation of the aircraft by a system of rigid masses 
elastically connected, as show the works of Traill-Nash [25, 26] 
and Hunh [27]; more recently Przemieniecki [28] used the sub- 
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structure concept with the boundary interaction forces as unknowns 
in the analysis. An interesting alternative method was presented 
by Schmitz [29] ; the equations of motion of the flexible airplane 
are written in the form of a nonself-adjoint system (including 
structural damping) , which is then treated by an eigenfunction ex- 
pansion. Although the analysis presented in Chapter III does not 
include structural damping, the method can be readily generalized 
to the non- conservative case. Appendix II provides in fact an 
analysis of the modal damping, also a Rayleigh Ritz component 
modes analysis using Lagrange Multipliers. 

The study of Chapter IV is motivated by the increasing interest 
developed m the past recent years for the concept of a yawed wing 
aircraft. In a recent article [30j , R. T. Jones from the NASA Ames 
Research Center suggested that an oblique wing configuration 
can be very suitable from an aerodynamic point of view, and 
encouraged in a more serious investigation of the problem. 

The concept was found interesting and resulted in a large 
number of research proposals from Universities and Companies. 

The Boeing Company [31] is already engaged in studies concerning 
the static aeroelastic phenomena associated with the oblique 
wing aircraft, and Stanford University has submitted a proposal 
for the investigation of the effects of asyaunetry on the 
dynamic stability of aircraft . As a first step in the aeroelastic 
analysis, the study of Chapter IV provides the structural material 
frequencies and normal mods for such configerations . 
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For straight-wing airplanes where the high aspect ratio 
assumption for the wing and the tail is not satisfactory, a theo- 
retical treatment is developed in the Appendix IV; the latter 
components are represented by plates with uniform properties. 

As a preliminary analysis, the study of the Capter V shows 
how successfully the method applies for the two-dimensional problems 
of plate vibrations. Two different examples are chosen for this 
purpose; the first example treats the problem of a point clamped plate 
and investigates the convergence of the natural frequencies towards 

those of a cantilivered plate, as given in Reference [35]. This 
analysis gives insight into the suitable values to be chosen for 
the parameters (initial nb. of modes used for the plate in x 
and y - directions, nb. of points of constraint along the plate) 
so that the natural frequencies obtained by the analysis of 
Appendix IV, are within a good degree of accuracy. 

The second example, which had a recent treatment in the 
literature [36] by the use of dual series expansions, is the 
problem of a simply supported plate with an internal rigid line 
of support. The present method is proven to be very satisfactory 
for such problems as shown by the convergence graphs obtained 
in the analysis. 
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CHAPTER II 

MODAL BEHAVIOR OF BEAMS WITH NON-UNIFORM PROPERTIES 
1 . Statement of the problem and basic assumptions. 

The present chapter deals with structures of the beam type, 
as illustrated below. 

4* 



We consider beams of length 21. along the x axis and we denote 
by E Young, modulus of elasticity and by I the structural mo- 
ment of inertia about the neutral axis , 

I = // Z^dy dz 

cross sectional area 

The bending stiffness El will be considered variable along the 
span; the variations can be due to the geometry of the beam 
(non constant thickness or width) , or to variations in the 
elastic properties (variable E along the span). The theoretical 
development is general and applies to all the possible cases 
mentioned above; numerical results are given in some particular 
cases . 

The basic assumptions of the present analysis are the following: 
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(1) The elastic axis of the unbent beam (as classically defined. 
Ref. (2)) is a straight line which bends into a plane curve during vi- 
bration. 

(2) The elastic axis coincides with the line of the centers of 
gravity of the cross sections, hence bending is uncoupled from twisting. 

(3) Initial stress in the axial direction is neglected although 

L 2 

conceptually it could be included. An additional term T / w' dx 

2 -L 

would appear in the potential energy; the details of the constraint 
conditions and effects on the distribution of internal shear and bending 
moment require more study. 

(4) The lateral displacements of the beam in vibration are 
small compared to its thickness and the slope is small compared 
to unity (this justifies a geometrically linear theory) . 

(5) The elastic constitutive equations are those of a linear 
isotropic material, in other words the stress - strain relation 
is given by Hooke’s law. 

(6) Shear deflection and rotatory inertia are neglected. 


2 . Theoretical background . 

The problem is considered in a variational formulation ; 
w (x,t) being the deflection of the beam at time t and station x, 
its kinetic and potential energies are: 

L 

r ~ 

T = I / (3w (x,t) y dm 

2 / at 

-L 

where dm is the element of mass per unit length. 
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u = 1 f EI(x)( 3 2w (x,t) ) 2 dx 
2 

J x 

-L 

When El has a complicated dependence on x, the last integral is 
very difficult to handle. 

The present method proposes a solution to the problem by 
conceptually disassembling the structure into sub-beam components 
of uniform bending stiffness. 

The rigourous justification of this approach is based on the 
mathematical approximation of a continous function defined on a 
compact set by step functions, in the following manner: 



For each of the substructures, the modal displacements 
are sought in a Rayleigh Ritz expansion and the continuity con- 
ditions are enforced by means of Lagrange multipliers. 

3. Derivation of the frequency equation . 

The following notations are now introduced: 

(EI)^: constant bending stiffness of the nth. component, n=l , N 
: half mass of the nth component 
: half length of the nth component 
dm^ 11 -* : element of mass/unit length relatif to the nth component 

wC ^ (x,t) : bending displacement of the nth component 
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: non dimensional frequency parameter of a free-free beam 
of uniform stiffness. 

<(> (x) •' corresponding normal mode (free- free beam) 

q^(n)(t): generalized coordinate of the nth component relative to 
the normal mode number i . 

The modal solutions are sought in the form: 


w (n) (x,t) = 


00 

£ 

i = 1 


( 1 ) 


q. (n) (t) *.(x) 

0) (n) 

for - i <_ x < t 

Subsequently, the total kinetic energy of the structure is: 

£(n) 

2 

(x , t ) dm (n ) (2) 

_f(n) 

By substitution of eq (1) and use of the orthonormal properties 
of the free-free beam modes, (2) can be written in the form: 


N 

T - — £ 

2 n=l 


T- ‘ ! “ H™ q.W(t) 

n=l i=l 1 


(3) 


Similarly, the potential energy of the entire beam is: 


£(n) 


U 


= j £ (EI) Cn) f ( ( - n - ) (x,t) ) 2 dx 
n=l J . . 3x^ 


_£(n) 


Substituti on of eq (1) and use of elementary beam equation yields 

(4) 


N 


£ ■ 


v 4 „ (n) r *>2 


(Ei) k: q . w (t)' 

n=l i=l 3 1 1 


t (n} ' 
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The be am- components are required to satisfy the following con- 
tinuity conditions : 

a) continuity of the bending displacements at the interfaces: 

w Cp) (f(P)), w (P + 1 )(_£(P +1 )} (5) where £^0® 

b) continuity of the slopes: 

l ( 6 ) 

3x W 9x 

At this point it has to be pointed out how convenient is the 
choice of the free-free beam modes as "trial" functions in the 
Rayleigh Ritz expansion. 

In addition to the fact that they allow a great flexibility 
in expressing the various constraint conditions, their choicewsatis- 
fies already the continuity of the second and third derivative of 
the w a t the connection points . 

In the case in which the initial beam is not free at both ends, 
corresponding additional constraints have to be added, as illustrat- 
ed in section 4.b) 

For the time being, the initial beam will be considered uncon- 
strained without loss in generality. 

Equations (5) and (6) can be written after suitable non-dimen- 


sionalization: 
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w( p ) (Jj s w(P + 1 )(-l) 

forP= 1, N-l 


(P) 


w c 


(P) 


( 1 ) 


1 

(P*D 


,-CP+i) 

w ' • 


(- 1 ) 


and further, by use of eq (1) : 


; OO | 

| * q t (P) Ct) ♦ i (l) - q A CP+1) (t) *..(-•!) = 0 ] (7) 




E 

i=l 


i p = i • 

q 4 (P) (t) »'•(!) - q i^^Ct) «E'(-1) = 0 | (8) 

i>Cp) t Cp + D 1 j 

Subsequently, the Lagrangian of the structure is: 


.N-l 


N 

df-i E 

" n-l 


E 

i=l 


M (n \ (n) (t)‘ 


- 2 Z 


(El) 


n=l i=l l 


(n) 


(n) 

3 


k. 4 q. (l °(t) 2 


N-l CO 

+ s 2 xn (q/^C^ + .Cl)- q. (n+1) (t) <J> i (-l)) + 

n=l i=l 

N-l CO 

E 1 Pn ( qj ^Ctj ^^d) - q 1 (n * 1) (t) 
n-l i = l £ Cn) t (n+1) 

where ^ n >P n are Lagrange multipliers. (9) 

It is interesting to note at this point that through the use 
of the Lagrange multipliers as they appear in (9) , the present 
method gives not only information about the modal behavior of the 
beam, but in addition provides results for shear forces and bending 
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moments . 

Indeed, it is known in Classical Mechanics that the Lagrange 
multipliers can be physically interpreted as being the generalized 
forces demanded for the maintenance of the corresponding constraint 
conditions. In the present case, the values of the and 

n=l,...N=i give respectively the distribution of the shear force 
and the distribution of the bending moment along the span of the 
beam at the points of constraint. 

The Lagrange's equations derived from (9) are: 

* M^’q ^(t) + (El) t 1 ) kiVi^Ct) - *^(1) “iij 0 

~? 1} 3 

• M (n) qi (n) (t) + (m Cr ° kiV^Ct) - ViU) + 


Vl V" 1 ) ” u n *V l) + Vl *'i(“ 1 ) * 0 

F n ) "7TS) 


n =2, N-l 


» 


M^'qC N )(;t) + (£Ij_ 

(N)3 

I 

*!hl ’ 0 

“IR-d 

f 


k i 4 q/^ct) + x^j^c-i) + 


00 


qi Cn) (t) ^(1) - q i Cn+1) (t) <f>.(-r) = 0 


0 


z 

i=l 


n=l,...N-l 
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» E q x {n) (t) #' t CD - q 1 ( "* 1) (t) t'.(-l) ■ 0 n»l,... N -l 

i.l j (n > J ( n+1 ) 

The last two equations are in fact the constraint equations 
expressing the continuity of the displacements and slopes. 

The following harmonic time dependence is assumed in the 
modal analysis: 

q . <"> t t) = e 


i 




wt 


x p ct) .. X CEIii 15 ^ 

( 1)3 


i 


(lit 


P (t) = p (El) 1 e 

n ( 1)2 

t 


£ T' 


:> ' 

— o)t 


A-," ' 3 C < ' ' 4 


and the following notations are introduced: 


m 

n 


= mC1) 
M (n) 


(El) 


(n) 


^n = 


M 


(n) 8n) 3 


(El) 


( 1 ) 


{f 1 ) 3 


= f (1) 

Y n Th) 


Fi^ 1 ) “ 
Fi( n ) = 



4 


n=2 . . .n 
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As a first step in solving the Lagrange’s equations, the 
generalized coordinates are expressed in terms of the Lagrange 
Multipliers from the equations of motion per se; 

( 10 ) 


( 11 ) 



Then (10), (11), (12) are used in the constraint equations 
in order to obtain a set of 2N-2 linear algebraic equations in 
terms of the unknown A n and This set of eigenvalue equations 

is now written for convenience in matrix form: 

(13) 

where A is the 2N-2 dimensional vector having for components A n 

(n=l, ...H-l) and (n=l, . ...N-l) and A (u) is the symmetric 

matrix of order 2^-2 defined by its elements A as follows: 

nm 
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r n=l : 


1 1 


2 2 


i=l B (1) _ (2) 


A = - T. 
1 ? 


i = l 


HLoi: Cl) <fc . (-1) 

2 1 ( 2 ) 1 

* • 

1 


A = 0 
1 m 


3< m <N-1 


‘l N = a ? X 2 ) 


F. 

i 


CD 


, £ sui^iC'lAii; 

i = l - . 1 


A = - 2 m Y a' 

1 N+l 


F. 


C2) 


A =0 
1 m 


N+2 < m < 2N-2 C.N>4) 


» £ = 2, . . . N-2 


A =0 
n m 


1 £ m <_ n-2 


A = - S ^-i CD *i(-l) 
n n-1 - -i “ ^ 1 

F .fn) 

i 


2 2 

A n n = E C m n _(D + m n+1 ^ (-1) ) 

U1 p. (n) F 

i 
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An n+l= ~ Z %+l hULtiLD 
i=l v (n+1) 

i 


A = 0 
n m 


n+2 < m < N+n 3 


A. - - i m n Yn ? i (1) 4 > 1 (-1) 

n N+n-2 i=l * 


00 


F. 

l 


N«-r J ( ™n T n iilllily * \ n iLA 


F. 

1 


(n) 


\ N+n“ “ m n +i Y n+ i 


F. 

1 


(n+1) 


A =0 
n in 


n+N+l <. m s 2N-2 


• n -N-l 

\ 

A|\j_ j m = 0 li m ^ N-3 


(N>_4) 


^N-l N-2 1 "Vj.i 1-j Cl)<t> . (-1) 

i_l 11 X i I 


F. 

l 


CN-1) 


^-1 N-l = 


1-1 F tN_1 ) ^ ^TN) 


F. 

l 


F. 

l 


N-l m = 0 


m £ 2N-4 


zm^izD ) 

(n+l) 

F. 
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N- 1 2N-3 


‘ViVi *1X12X1 ty 

F CN-1) 


Vi w-2 - . E , ( Vi y n-i *1x12x11? + - N Y N *<£ 4 Xfci? ? 

1=1 n CN-1) “ c (N) 


F. 

l 


* II..". N 


A N m = £ m N already defined for l£m<N-l 


N N 


^ CD 
i=1 f^ 


2 2 

+ m 2 y 2 £z£> ) 

n C2) 


N N+l 


Z IBoY 2 
i=l 


iLiCilililiil 

r ( 2 ) 


A N m = 0 


N+2<m< 2N-2 


* n = N+i, ... 2 N -3 ; p = n .N+i 


* = A already defined for l<m<N ■ N+l <n<2N-3 

n m mn 


4 n P = il/ Vp *iX« 4 llXii * V1V1 ll-tlXXd) ) 


F. 

1 


(p) 


F. 

1 


Tp+T) 
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'n p*l* ' . £ , %.l Y p+! hiniliid) 

1=1 ‘ F.'P* 1 

1 


A_ = 0 
n m 


p+2 <m< p+N-3 


'n N*p-2* - . E , Vp .♦ilill’ili) 

1=1 (P) 


F. 

1 


2 2 


4 n N.p-1= .1 t «tV Hi-S * VlVl ility ’ 

1 =1 r CP) ( 0 * 1 ) 


• F. 
x 


Si N+p“ " . £ 111 p+1 P +1 

i=l 


iOLiliLi^ 

F. (P+1) 


m - 0 


MfPfl < m < 2^2 


2 2 2 2 


^N-2 2N-2 ’ i=1 (VlVl +m N Y N— i^p]) 


F. 

l 


Equation (3) leads to the frequency equation 



(14) 


which is the necessary and sufficient condition for existence of 




- 21 - 


solutions p„ n = i N-l 

This equation cannot be solved analytically except maybe for 
very small values of N ; numerically, the fastest and simplest 
method to use is a trial and error approach. The left hand side of 
the equation (14) is plotted against the variable w ; if two successive 
values of w make the left hand side of (14) change sign, then one of 
the natural frequencies should be between these two oi’s . 

This is a quite straight forward operation except for the fact 
that care should be taken to decide whether the change in sign cor- 
responds to a valid solution (i.e. a natural frquency of the beam) 
or to a natural frequency of one of the components. Such types of 
situations are illustrated below: 



The accuracy of the natural frequency can be made as high as 
necessary, provided that the 2 successive w's are close enough and 
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that the values of the left hand side of equation (14) are close to 
zero. 

For each natural frequency the Lagrange multipliers are then 
determined from equation (13) up to multiplicity, for instance the 
ratios n=2,...N-l and u n n=l,...N-l 



can be determined. 

Subsequently the generalized coordinates for each component are 
determined from eqs . (10), (11), (12) and then the corresponding mode 
shape as given by equation (1) . 

4 . Numerical analysis. 

a) A simply supported beam of variable geometry. 

The present analysis treats the case of a stepped beam, simply 
supported at each end as shown in Figure 1. 

The variation in the bending stiffness is due to the variation 
in the thickness of the beam and the geometry of the problem suggests 
to consider two components. As mentioned previously, the theory 
extends to the cases where the initial beam is constrained at the 
ends; in the present problem two additional constraint conditions 
have to be reinforced, namely 

(-£) = 0 
w^ Uj = o 

For the sake of simplicity it is assumed that 1 = 1 and that 
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the cross - sections of the beam are rectangular, without great loss 
in the generality of the problem. The analysis is developed in two 
steps: the first step provides a check on the method for accuracy 

and convergence, the second step nrovides the functional dependence 
of the fundamental frequency with the ratio h^ . 

wCD 


In the particular case in which h^) = \ > the problem reduces 

jrrn 

to that of a uniform beam simply supported at the two ends for which 
the modal behavior is known analytically; in particular the first 
three non-dimensional frequency parameters are tOj = tt , =4tt , =9tt 
In a first analysis, the convergence of the natural frequencies of the 
beam towards the above values versus the number of modes used per com- 
ponent has been investigated, and the numerical results are presented 
below in tabulated form: 


Table 1: First natural frequency 


# modes used total # of 
per component modes 



percentage 

error 


3 

6 

11.94 

21.6% 

4 

8 

11.05 

12 % 

6 

12 

10.16 

3 % 

8 

16 

9.918 

0.5% 

10 

20 

9.898 

0.2% 
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Table 2: Second natural frequency 


# modes used 
per component 

total # of 
modes 

U) 

2 

percentage 

error. 

3 

6 

53.68 

36 . 7% 

4 

8 

45.01 

14.9% 

6 

12 

41.06 

4.7% 

8 

16 

39 . 835 

0.9% 

10 

20 

39.64 

0.4% 

Table 3: Third natural frequency 


# modes used 
per component 

total # of 
modes 

(Or* 

3 

percentage 

error 

3 

6 

151.89 

71 % 

4 

8 

103.93 

17 % 

6 

12 

94.16 

6.3 % 

8 

16 

90.25 

1.58% 

10 

20 

89.71 

0 . 87% 


The results given in Table 1 are plotted in Figure 2; it can 

be seen that the first natural frequency converges rapidly towards 
2 

tt and that the rate of convergence is very steep especially at the 
begining. Figure 3 represents a graph of the percentage error 
in versus the number oj: modes used per component; it can be seen 
by extrapolation that in order to reach the accuracy in corres- 
ponding to 10 modes per component, it is necessary to consider in 
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this case two additional modes for each segment. 

Figure 4 provides a picture of the rate of convergence in the 
third mode; in this case 14 modes are needed per component in order 
to reach an accuracy of 0.2%. 

This analysis is somewhat similar to the convergence analysis 
in the classical Rayleigh Ritz method; moreover all the numerical 
natural frequencies obtained are higher that the exact frequencies . 
Another important convergence parameter is N (the number of com- 
ponents into which the structure is conceptually disassembled); the 
rate of convergence versus N in the fundamental mode has been in- 
vestigated in the check case of the uniform beam. The results are 
shown in Figure 5; 3 free-free modes are used for each component. 
Although the graph shows rapid convergence, the results cannot be 
interpreted in a self-contained manner. The reason is that even 
with a fixed number of modes per component, the fact of increasing N 
has the effect of increasing the total number of modes considered 
and hence of improving the convergence. 

As a consequence, the total number of modes considered has to 
be chosen as a common parameter although it is not controlled direct- 
ly in the method. The comparative rate of convergence is shown in 
Figure 6; it can be seen that for the same total number of modes, 
better results are obtained when for a fixed number of components 
(N=2) the number of free-free modes per component gets larger. The 
reason in this case is that there is "no need" in fact to increase 
the number of components and thus making N>2 adds no physical improve 
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merit. From the numerical point of view, increasing the number of 
modes corresponds to simpler operations than increasing the number 
of components. By considering additional modes per component we 
are only increasing the number of terms in the series giving the 
Ai j (w), while an additional number of components correspond to an 
increase in the size of the matrix A (to) . 

Operations of the first type are handled by computer with higher 
accuracy than computations of determinants for instance, for large 
size matrices. 

It should be mentioned that the present convergence analysis 

applies only to this check case and therefore gives only an insight 

to how the convergence problem should be handled; general conclusions 

about convergence will be developed in the next paragraph. 

f 21 

The study of the case h w = 1 is completed by Figure 7 and 8. 

h™ 

Figure 7 shows the first mode shape and the orders of magnitude of 
the shear force per unit length and bending moment per unit length 
at the mid-point; Figure 8 shows similar results for the second mode. 
There is very good agreement with the classical sinusoidal shapes. 

This concludes the study of the check case; in the second part of 
the analysis the change in the fundamental frequency with the thick- 
ness ratio is investigated. 

The values of the non-dimensional frequency parameter oj^ are 

plotted versus the ratio H = in Figure <5. 

h' 2 > 

For ^ ^ the frequency parameter is based on the character- 
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istics o£ the first component, for < h^) is based on the 

characteristics of the second component. The reason for this choice 
will appear clearly as we proceed. 

The Figure shows that for > h^, as H increases the funda- 

9 

mental frequency decreases from the value 9.898 (i.e. it ) to zero. 

The first value correspond to the previous analysis, the second 
value correspond to the limiting case of a uniform beam simply sup- 
ported at one end and free at the other. In this limiting case the 
fundamental frequency is zero and corresponds to a rigid body rotation 
of the structure. The monotonic variation of with H checks our 
physical intuition; given the two beams illustrated below: 


it is expected that the resonant frequencies of the right-most beam 
be lower than the frequencies of the beam at its left. The configu- 
rations for which h^ ^ < h^ ^ correspond to the previous configurations 
by mirror- symmetry and hence the physical resonant frequencies should 

be the same for two beams of respectives thickness ratios Ho and 1 . 

Ho 

If we switch the "reference beams" on which the material frequen- 
cies are based, the same argument holds for the non-dimensional fre- 
quency parameters and this provides additional checking for the method 
and numerical work. This is the reason for which the fundamental 
frequency parameter is based on different components in Figure 9. 
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The Figure shows very good correspondence of the frequencies in the 
mapping of [1, »] into [1, o] . 

An additional checking of the quantitative results of Figure 9 

can be made by considering Figure 9 (bis) . The variation in based 

on h^ is represented as function of H = h^ and it is seen that 

h™ 

for H = 0 the obtained uniting value for to- is 15.6. This limiting 
case is that of a uniform beam simply supported at one end and clamped 
at the other, for which the first natural frequency has the exact 
value of 15.418; i.e. only 1.2% lower than the numerical value pre- 
dicted by the method. 

b) A non-uniform helicopter rotor blade (cantilivered beam). 

The choice of the present model is related to a more general 
investigation of the effects of distributed blade properties on the 
stability of hingeless rotors, done by Professor Dowell and Professor 
Curtiss at Princeton University for the NASA Ames Research Center. 

This project is of recent interest (1973) and includes extensive 
theoretical and experimental developments. 

In the first part of the investigation, the main interest has 
been focused upon a theoretical study of the structural modeling and 
determination of natural frequencies of a rotor blade with large 
changes in stiffness distribution. The classical Rayleigh Ritz method 
using an assumed Duncan polynomial shape (Ref. (38)) fails to give 
accurate results if the blade non-uniformity is severe; therefore 
the results of the method developed in the present chapter will be 
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compared to the results of a more accurate method developed by 
Professor Dowell, which is an initial value approach to boundary 
value problems (IVBV) . An experimental value of the fundamental 
natural frequency is also available. The non-uniform rotor blade 
under consideration has the following characteristics: 


Total 

length : 48 inches 



Span 

bending stiffness 

mass per unit 

length 

0 - 4 in 

3 x lO^lb x in 2 

0.52xl0" 3 lb x 

2,. 

sec /m 

4 -14 in 

24,400 lb x in 2 

0.66x10 4 lb x 

2,. 

sec /in 

14 -48 in 

3,276 lb x in 2 

0.35xl0‘ 4 lb x 

2.. 

sec /in 


The fundamental natural frequency of this blade predicted by 
different methods is given below: 

Rayleigh Ritz Myklestad IVBV Component Experiment 

method using (Ref. modes using 

Duncan polynomials (49) ) Langrange mult. 

“l 6.7 cps. 4.45 cps. 4.36 4.37 cps . 4.4 cps . 

cps . 

The IVBV and component modes method using Lagrange multipliers 
are seen to provide the most accurate results and shall be considered 
for further comparing the next two resonant frequencies . A summary 
of the numerical results is given below in tabulated form. 


Natural 

frequency 

Prediction of 
the component 
modes - Lagrange 

Prediction 
of IVBV 
(Ref. f49V> 

4 . 36 cps . 

Percentage 

disagreement 

fundamental 

4.37 cps. 

0.23% 

second 

24.6 cps . 

24.2 cps. 

0.16% 

third 

61 cps. 

59 . 3 cps . 

1.12% 
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The fundamental mode shape is shown in Figure 10 as predicted 
by both IVBV and CML methods; the figure shows good quantitative agree- 
ment - the largest deviation of 0.07% occurs at around 31 in. from 
the clamped end. Similar results are shown in Figures 11 and 12 for 
the second and the third mode. 

5 . Discussion on accuracy and convergence of the method. 

The question which most naturally arises when the method is 
to be applied to an engineering problem is that of the optimal choice 
of the number of components and of the number of modes per component, 
in order to obtain good numerical results with minimum computational 
effort. It is obvious that no general answer including all possible 
cases can be given to this question; furthermore the complexity of 
the problem is such that for most of the cases of interest there is 
no analytical answer. 

r 

The purpose of this paragraph is to develop a two-step procedure 
that can be called "modes -component selection" in order to provide a 
satisfactory answer to the above question. 

The procedure is illustrated on an arbitrary case. 

The stiffness distribution along the beam span is assumed to 
vary as follows: 



x 
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The method consists in: 

(1) Determining first an order of magnitude of the natural frequencies. 

For this purpose few components are considered (in order to keep 
the size of the eigenvalue matrix small) with no restriction on the 
number of modes per component. The number of components is directly 
related to the rate of change of EI(x); for instance region (1) can 
very weli be represented by only one component, region (3) could have 
similar representation within reasonable approximation, but more than 
one component should be used in region (2) . 

For the purpose of obtaining only orders of magnitude of the 
frequencies, two or three components should suffice in region (2). Thus 
in a first analysis a 4 or 5 components approximation is taken for 
the initial beam, which keeps the size of the matrix A of the order 
of 6 or 8. 

Additional saving in computer time can be made by noticing that 
the component in region (3) is much stiffer than the other components 
Thus most of the bending will occur in regions (1) and (2); region (3) 
is expected to have predominantly rigid behavior. As a consequence, 
less free-free modes need to be used for that component than for the 
others . 

The orders of magnitude of the natural frequencies are now 
determined, and we proceed to the second step. 

(2) Improving the solution in the neighborhood of each natural frequency. 

This is realized by first improving the modelling of the beam. 
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i.e. increasing the number of components in regions (2) and (3) for 
instance . 

A reasonable way of improvement would be for example to consider 
about six components in region (2) and two components in region (3) . At 
this point the size of the eigenvalue matrix is increased and thus 
more computer time is needed. A way of compensating for this dis- 
advantage is to make use of the information provided in the previous 
analysis in order to decrease the number of modes used per component 
to the strict minimum needed. 

The selection of the modes which are to be used for each com- 
ponent is made on the basis of the order of magnitude of the corres- 
ponding natural frequencies compared to the orders of magnitude 
found in the previous analysis. 

Such a selection can decrease the number of modes used quite 
considerably without any significant loss in accuracy. 

Mathematically, this selection can be justified by the fact that 
only the reasonably small ^ are to be kept in the analysis. Their 
inverses provide the main contribution for the expression of the 
A.. . (w) and hence for the frequency equation. 

It can be anticipated that if problems of non-uniform plates 
were to be treated, a more accurate solution would be obtained by 
increasing the number of components rather than the number of modes 
per component. The reason is that generally only approximate numeri- 
cal solutions are known for the modes of an unconstrained plate and 
that the precision decreases as the inodes get higher. Therefore it 
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is better to increase the total number of degrees of freedom for the 
plate by increasing the number of componens rather than considering 
higher (approximate) modes. 
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CHAPTER III 

MODAL BEHAVIOR OF A SYMMETRIC FLEXIBLE AIRPLANE 
WITH HIGH ASPECT RATIO BEAM-LIKE COMPONENTS. 

1. Statement of the problem and basic assumptions 

The type of structure considered in this chapter is a flexible 
airplane with a high aspect ratio wing. 

There are many various possible models which can represent such 
a structure, every one of these presents advantages and inconveniences. 
When the problem is modeled very close to reality, the components are 
represented by fairly complicated elastic structures for which the 
normal modes and frequencies are only known with a certain degree of 
approximation. These initial errors propagate up to the final result 
and there is no guarantee that the modal behavior is determined more 
accurately in such analysis. On the other hand if a simpler mathe- 
matical model is chosen to represent the structure, the natural fre- 
quencies and mode shapes can be determined more accurately but the 
resemblance between the physical and the mathematical problem may be 
slight. The choice between the two alternatives has been made on the 
basis of simplicity and clarity of the formulation and taking into 
account the fact that the errors due to physical simplifying assump- 
tions are easier to determine than mathematical errors propagating 
with the computations. Therefore a simple mathematical model has 
been chosen to represent the problem. The aircraft is considered to 
be made of substructures (fuselage, wing, tail) which have been apDroxi- 
mated by beams of uniform stiffness as illustrated in the following 
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If "yawing-type" motions (i.e. bending in the y-direction) were 

to be considered the rudder should be included as well; for our pur- 
pose (transversal bending of the structure), it has been neglected. 

The basic assumptions are summarized below: 

(1) The airplane structure is considered to be of the beam type. 

(2) The displacements of the beam components are governed by the 
classical one-dimensional equations (in particular, Bernoulli- 
Euler equation for the lateral bending) . 

(3) The geometrical and elastic characteristics are constant along 
the span for each component. This last assumption can be relieved 
using the approach of Chapter II. 

2 . Rigid body analysis (small amplitude motion about the initial con- 
figuration) 

Notations : 

R : Gallilean frame of reference 

o 

G: Center of gravity of the structure 

T (G/R ) : Acceleration of the center of gravity with respect to R 
o 


o 
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iJj: yaw angle in the rigid motion 

9: roll angle in the rigid motion 

<J> : pitch angle in the rigid motion 

(ip, 0, <j> are the classical Euler angles) 

ft ((S)/R ): angular momentum of the structure with respect to R 
o o 

Ij^ I 2 , I ^ : principal moments of inertia 

In absence of exterior forces Newton’s law yields. 

-*■ -»• 

(total mass) . r (G/R q ) = 0 (15) 

This relation decouples the translations in the three directions 
from the rotation motion abound the center of gravity. 

For zero exterior torque the angular momentum equation is simply: 

■+• 

d jt ( (S) /R ) = 0 (16) 

dt 0 

+ 

Because of the symmetry of the structure, x y z are prin- 


cipal axis of 
and therefore 

inertia with respectif moments of 
the linearized equations of motion 

intertia 1^, I 2> I 
of the rigid body are 


h 0< 

= 0 

(17) 



= 0 

(18) 



= 0 

(19) 


This set of equations shows that the roll, yaw and pitch motions 
are all decoupled from each other and that the corresponding natural 
frequencies are zero. 
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The same type of considerations apply in the elastic domain be- 
cause of the symmetry, as will be seen in paragraph 3. 

The study of the decoupled motion of "yaw-type” is not included 
in the present analysis, as mentioned before. 

3. Elastic modal behavior ; derivation of the frequency equation 

a) Statement of the elastic problem: 

degrees of freedom, derivation of the constraint equations 
As the motions investigated are of both roll and pitch type, 
each of the substructures has to be allowed to undergo transverse 
bending and twisting as well., Although it is sufficient to con- 
sider only rigid twisting for each component, the analysis will be 
developed for the general case in which the torsional stiffness 
of each conponent is finite. The previous case will then be a simpli- 
fication of the latter. 

We develop a Rayleigh Ritz analysis by defining: 

- for the fuselage : - L < x < L 

The bending displacement W (x, t) sought in the form: 

W (x, t) - 2 0 (t) * (x) C20) 

n-1 

The twisting angle. 0 (x, t) sought likewise: 

9 (x, t) = ? Pit) 9 (x) (21) 

n=l 

- for the wing : - < y < 

The bending displacement (y, t) : 

W (y, t) = ? q (1) (t) * (y) 

n=l n n 


C22) 
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The twisting angle 0 ^ (y, t) : 

0 (y, t) = I P C1) (t) 6 (y) (23) 

1 n=l n n 

- for the tail : - 1? 2 < y < 

The bending displacement W (y, t) : 

w 2 (y, t) = | =1 q n C2) (t) * n (y) (24) 

The twist angle © 2 (y, t) : 

© 2 V = I = i P n C2) Ct) 0 n (y) (25) 

The physical displacement of each component, respectively at y = + c, x 

x + C , x = x -x + C , is given by: 

11 2 12 

W phys ^ ^ r W ( x > + c '0 (x, t) (26) 

/'phys Cy> t ) = w ! (y> "t) + c i ©! t) (27) 

•V, p hys (y ’ ^ = W 2 (y ’ ^ +c 2 e 2 ty ’ ^ ( 28 ) 

Where c, c^ c 2 are the half chords of the respectif components. 

With suitable non- dimens ionalization x = L £ eq. (26) yields: 

" phys (;,t) . W «,t ) + £ 8 ({. t) (29) 

L L L 

and by similar changes of variables 

y = n and y = f 2 n for equations (27) and (28) re- 
spectively, they become: 
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and 


W 1 phys ^ 


W (n, t) c 

+ 1 

X' I- 


0j_ Ov t) 


(30) 


W 2 phys (n > = W 2 Cn > 


c 2 e 2 (n, t) 


(31) 


In agreement with the one-dimensional treatment of each substruc- 


ture, it can be assumed that c_, and c 2 are negligibly small. 


in fact small enough so that 

c_ 0 (£, t)<<l for any £ -1 < 5 < 1 

L “ 

0 ^ 0j (n, t) <<1 for any n -1 < n < 1 

c 2 0 2 (n, t) <<1 for any n -1 < n < 1 

J 2 

As a consequence, the corresponding second terms, will be neg- 
lected on the right hand side of equations (29), (30), (31). The 
physical displacements at the joints between the components must be 
equal, and this leads mathematically to the constraint conditions: 


W w 
phys 


(C,, t) = W 


1 phys 


(0, t) 


(32) 



or, by using (29), (30), (31) : 


W (q, t) - Wj (0, t) 


( 34 ) 


W U 2 , t) = W 2 (0, t) (35) 

The equality between the slopes in the x and y directions has also 
to be enforced as part of the constraint conditions and leads to the 
equations : 

0 U r t) = 1. 9W 1 (o, t) (36) 


0 U 2 , t) = 1^ 3W 2 (0, t) 

1’ 2 — 


(37) 


for the continuity of the slope in the x - direction, and similarly 


I 3W (C 1 , t) = - © x (0, t) (38) 

L H 


1 3W U 2 , t) = - 0 2 (0, t) (39) 

L 3£ 


for the continuity of the slope in the y - direction. 
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b) Derivation of the frequency equation 

The total kinetic energy of the structure is 

L L $ 

T = 1 ('ft 2 ( X , t) dm + 1 I © 2 (x, t) dx + 1 W 2 (y , t) dm, 
2 -L -L -I 1 1 


+ ' r % (7, t) dy 


+ 1 w 2 2 (y> t) dm 2 

~^2 


* J . 1 2 0 2 ^ dy ^ 

~ l 2 


C40) 


where I, I , X ? are the fuselage, wing and tail mass moments of 
inertia per unit length about the torsional axis. In particular 


1 = p £us U (y 2 + z 2 ) dy dz 

fuselage 

cross-section 


*1 * P wing (x2 + dx dz 

wing 

cross-section 


*2 ’ "tail '' . t" 2 * ** ** 

tail 

cross-section 


By substitution of the corresponding expansions for each of the 
bending and angular displacements and because of the orthonormal 



42 - 


property of the bending and twisting modes , equation (40) becomes : 



where 

M, Mj, M 2 are the half masses of the fuselage, wing and tail re- 
spectively (i. e. mass/unit length x half length) 

Likewise, the total potential energy is the sum of the potential 
energies of the components : 



L 2 2 ^ ? 

[ El / C9 W r dx + GJ / ,90. 1 dx + 
-L .2 -L 1 " 1 


(. EI)^ f (- 9 ^ yf dy + (GJ)^ / ^ 30 ^ dv + 

-ij 9y -fj 9y 

£ £ 2 

C ei ) 2 t 2 (»h2) 2 4y * CGJ) 2 / .DO, 2 ] (42) 

-i 2 9y 2 -< 2 V 

where 

^2 are the CC) rres ponding structural moments of inertias 
about the neutral axes given by 

I = // z 2 dy dz 

fuselage 
cross section 
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I = // z 2 ' dx dz 

wing 

cross-section 

I ? = ff z 2 dx dz 
tail 

cross-section 

E and G are the Young and shear modulus of elasticity (related by 

G = E ) and J is a constant which value depends upon the 

2 (1+v) 

shape of the beam's cross-section. 

In particular, useful values of J are tabulated in Ref. (39) . 

(GJ) is the classical twisting stiffness. 

By expansion of the bending and angular displacements and use of 
the corresponding beam equation, equation (42) leads to: 



The constraint conditions have been previously derived and they 
are given by equations (34), (35), (36), (37), (38), (39). Substitution 
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of the displacements in terms of the component modes yields: 


z (On 00 V ^ ) - q n (1) (t) 4 > n (0) ) = 0 

n=l 


( 44 ) 


E , ft) S 2 ) - q n (2) CO * n CO) J - 0 

n=l 


(45) 


Z (P n (t) 0 n ( 5 ) - q n C1) (t) d> ’ n CO) ) = 0 

n=l * 

1 

oo 

2 , (P „ "> V «2 > - q„ (2) ♦'„(») > * 0 

n=l a 

f 2 

OO 

2 (On ct))' n ( 5j) * P n fl) (t) 6 CO) ) = 0 

0=1 L 

OO 

* (SELiil *' n (5 2 ) + p n C2) (t) 0 n (O) ) = 0 


(46) 


(47) 


(48) 


(49) 


The Lagrangian of the total elastic system is then 


w th 

L=T-U+Z A x ( i = constr. equation) 
i=l 


(50) 


where T and U are given by equations (41) and (43) . 

As mentioned in Chapter II, the Lagrange Multipliers provide 
effective information, namely: 

*1 30(1 X 2 give the shear forces developed at the joint between the 
fuselage and wing and at the joint between the fuselage and tail 
respectively, in motions of plunging type. 

^3 ^4 S* ve the torque developed at the two joints respectively. 


O 
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in motions of rolling type. 

X^ and Xg provide similar information for motions of pitching type. 

By use of equations (41), (43), (44) - (49), equation (50) 
leads to the Lagrange’s equations 

. El k/ Qh ft ) - Xj - X 2 n C« 2 ) - 

L' 5 

^5 ♦'„(?!) - \ = ° GD 

L L 

V L P n (t) + 1GJ1 (n- 1) 2 t? P n (t) _ X 3 © n (Sp _ 

L 4 

X 4 0 n (C 2 )s '° (51 bis) 


M 1 q n (1)ct) + ^ EI >1 k n 4 +A i* n ( °) + A 3 ♦ , n (0) = 0 


(52) 


h £ i p n (1)(t) + CG^^n-DV P n UJ (t) - x 5 6 n (0) = 


2.2 D (1) 


(53) 


;■ (2) 


M 2 q n ^' (t) + (:EI) 2 k n 4 + A 2 V°) + A 4 V n f0) = 0 ( 54 ) 


d‘(2) 


y 2 + ^ 2 (n -^ 2 I P n C2) ^ - A 6 6 n (0) = 0 


(55) 


The harmonic time dependence in the normal analysis is taken in the 
following form: 
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v 


1 
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El 


GJ 


(GJ). 


= 


ML 

(El) x 

M £ ' 
1 1 


= 


1 L" 

CHI). 


M 1 V 


*3 = 


*1 

~WT{ 

M l 

1 l l 


(HI). 


(GJ). 


*4 = 


M t 3 
2 2 


(El) 


1 


M £ 3 

1 h 


*5 = 


*2 ^2 

(El), 


M 1 V 


p C 1 ) _ 2 , , 4 

n n 


C2) . ,2 


i>2 C n- 1) 2 I 


(3) _ ,2 


n 


= a) - k 


n 


(4) _ 2 


n 


2 2 


w - tpj (n-1) 7r 


n (5) 2 ,4 

V * “ - *4 V 


F n (6) = » 2 - * 5 (n-l) 2 , 2 


4 
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The generalized coordinated are then expressed in terms of the Lagrange 
multipliers as follows: 


On - - »! ( ♦„ (Cj) ♦ * 2 W * X 5 Yl *' n ( 5l ) * 


( 1 ) 


n 


X 6 ’ (“) 
F CD 


P n = - n "l C x 3 9 n «i) * x 4 e n U 2 ) ) (57) 

p (2) 

n 

% (1) _ x ! ♦„«» * X 3 ♦ '„(( )) (S8) 

F C3) 

n 

F n (1) ■ - Vl X 5 (59) 

F C4) 

n 

— (21 — _ 

q n = v 2 C X 2 ♦n (0) *_ X 4 y 2 ♦'n (tl) 5 (60) 

r ( 5 ) 
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p (2) 

n 


Y 1 Y 2 V 3 X 6 6 n W 

F C6) 


C61) 


Note that and are proportional to 9 (0) and hence they have 

zero value for all the odd modes; as a consequence the angular twist for 
the wing and tail has an even distribution along the span as expected 
from the geometrical symmetry. 

The eigenvalue equations are found by substitution of equations 
(56) - (61) into the constraint conditions (44) - (49) and form a set 
of six linear algebraic equations in the given by: 


n 1 


( h ( *1 V 


(ep 


*2 *n <«2> * X S n +„ («1> C5j) 


CD 


n 


* A 6 y 1 > + h ^ } f62 ) 

p (1) p (3) 

n n 


oo 


S' (u, C 

n=l 


X 1 V g 2 } + x 2 + 4> n Cg 2 ) + 

F CD 


WV * , n Cg 2 ) } + P 2 ^n^ 

F & F C5) 

n r n 


= 0 


( 63 ) 
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CO 


z 

n=l 


(y i v i 


X 3 9 n 2 «l1 * *4 W W ’ * X 3 ♦’,>> > 

F C2) p (3) 

n n 


00 


2 

n=l 


f Vl 


X 3 W V S 2 5 * 1 + 


( 2 ) 


n 


V2 2 X 4 +'n 2 t 0 ” 


(5) 


n 


l ml (P 1 Y 1 C A 1 W + _ X 2 + V l 

F CD 


2 


*'n C V * Vl ♦ W *V*2> 


CD 


n 


* V2 >5 e n 2 (°^ ) 

P 


Z Bml lV l y l x i Vfp * X 2 W_ »'n« ? ) *X,.Y, 

7"^ 


♦W *Vf£ * Vl 

F (i) 


* Ws X 6 9 n 2 f°) > 

p C6) 


0 (64) 


C65) 


(66) 

0 


0 (67) 
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or, in tensional form 

-*■ -> 

A(w) . A = 0 (68) 

where A is the six-dimensional vector of components A^ and 
A(u>) is a symmetric tensor of order six which elements in matrix 
representation are: 



CO 

E 


n=l 


( Uj yhtp * 

F CD F C3) 

n n 


) 


>12 - l ml H W ^n ce 2 ) 

U Cl) 


4 13- 4 14 -° 


4 15 = ^ n, Tl 

F CD 

n 

oo 

S 16 = E , "1 *1 W ♦'n«2 3 

n=l — — 

P Cl) 


22 


OO 


= £ <v, 

n=l 


4> n 2(S 2 ) + m 2 ^n 2(0) 5 

F (1) F (S) 
n n 
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A 23 ' A 24 " ° 


25 


£ Yj r n Ccp 4> n (S 2 ) 

n=l 


(i) 


n 


A 26 


Li yiY? - w tvv 


Cl) 


*33 - E , (Y 1 V 1 9 n 2 t^ + ♦'n^ > 

n=l 


(2) n (3) 


n 


n 


34 


= Z Y 1 V 1 W W 


n=l 


( 2 ) 


n 


35 


A 36 " 0 


44 


= Z =1 Cy 1 V 1 V CS 2 ) + ^2 Y 2 2 ♦'n 2 ^ } 


( 2 ) 


C5) 


n 


n 


^45 


*46 - 0 


i55 * n -1 *’n 2(t 1 } * Vl e n 2 (°^ ) 


( 1 ) 


(4) 


n 


n 
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4 S6 * £ , *W 

n=l 


S 66 = ^ 1 . ( “l lr l 2 ♦'n 2cs 2 ) * ^2 v 3 . 5 


Equation (68) has solutions only for the values of u for which 


det ( A(w) ) = 0 


- (69) 


and these values are precisely the non-dimensional resonant frequencies 
of the structure. 

Eq (69) can be written more explicitely. 


A 11 

A 12 

0 

0 

A 15 

A 16 

A 12 

A 22 

0 

0 

A 25 

A 26 


A 34 '° 


A. . 0 

44 


4 1S 4 

25 

0 

0 

A 55 

A 56 

4 16 4 

26 

0 

0 

A 56 

A 66 

or, by simple algebraic 

manipulation. 



A 11 


A 12 

A 15 

A, 

^33 A 44 - A 34) 





, 


A 12 


A 22 

A 25 

A, 


A 16 


A 26 

A 56 

A < 
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The decoupling between the rolling and pitching type of motion 
found where the rigid behavior was investigated is also valid in the 
elastic domain as shows equation (50) and therefore the two types of 
motions are to be studied separately. 

The natural frequencies of the motions of rolling type are solutions 
of the equation 

( ‘33 4 44 ' 4 34 1 ■ 0 ; (71) 


In such motions the fuselage undergoes only angular twisting and 
the wing and tail undergo anti- symmetric bending. 

The pitching type of motion has natural frequencies given by ; 


11 

A 12 

A 15 

A 16 


12 

A 22 

A 25. 

A 26 

= 0 

15 

A 25 

A 55 

A 56 


16 

A 26 

A 56 

A 66 



The fuselage bends only, while the wing and tail bend and twist as 
well . 


c) Accuracy of the method 

As part of the analysis relies on numerical computations, the 
question of accuracy is a consequence of precisely these computational 
errors . 

The sources of error are: 

1) The fact of considering a finite number of modes for each component 


✓ 
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in the Rayleigh-Ritz analysis 

2) the fact of solving with approximation the frequency equation, 
det (A (a) ) = o 

3) the fact that the various computer subroutines used in the numerical 
analysis have errors inherent to the way in which they have been 
written . 

This last type of error can be eliminated by using Double Precision 

in the computation or Extended Precision if needed, and hence will not 

be accounted for in the accuracy analysis . 

(N) 

Let det A.. (to) be the approximation of A. . (to) when N modes 
lj " lj 

per component are considered. Then 


A ij = A iT Cu3 + E ij Cw) (73) 

where (oj) is mainly the next-order term in the serie giving A_ (to) i.e. 
the term due to considering the N*= + 1 component modes. These errors 
propagate into the determinant in a way which can be determined by 
substituting equation (73) into equations (71) and (72) respectively. 

For the motions of the rolling type this substitution yields: 




(N) 


- A 


34 


(N)2, 


33 


M e + A M £ 
44 44 e 


33 - * 


higher order terms = 0 
which can be written in the form: 
Det M + E M = 0 


( 74 ) 
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where 

net OO =4 IM 4 #«■> - 4 

33 44 34 

is the determinant actually computed in the finite mode - approximation and 

E fo) . 4 ft.) , A c*3 . 2 , 0*3 

33 44 44 35 34 44 

is the truncation error. 

Hence, instead of solving the exact frequency equation, one solves 


in fact 


Det ^ (to) = 0 


Furthermore, the numerical solutions a> n do not satisfy equation (75) 
exactly, but satisfy in fact 

Det ^0 = e (76) 

where the values of e are known and inherent to the trial 
and error process. Therefore the question that naturally arises 
is: given the numerical errors listed above, how different is 

exact solution co n ^ from the solution ^ actually 


Let be an exact natural frequency satisfying equation (71) t 

w the approximated value of lj * 


Det C03) - A 33 (oj)A 44 (oi) - A 34 2 (m) = Det (to) + E ^ ^ (oj) 


A 1-term Taylor expansion of Det (to ^ e - 1 ) 


gives : 


Det (tu o ^ *) - Det (to o ) + (u> o ^ - co o ) 3Det 


9(0 I (Jj =6 
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fe] 

and since Det (tu v J ) = 0, 



(77) 


All terms in the right hand side of equation (77) are known: 
£., as mentioned before is known in the trial and error process of 
solving Det^(<i>) = o , 

E ' N) <“o> =■ 4 33 <N) t“o’ £ 44 ( “c) ' 2i 3 4 CN ’ 


E 34 ( V> 


(N) 


can be explicitely computed from the values of the A^^ ; and 

which are known. 


3Det 


(N) 3A (N) (N) f w . (N) f • 3A.,(N) 

_ I _ I , A 44 O q ) + A 33 (u 0 ) i 44^ 


3o) oj=oj 

o 0 


3(0 CjU = DJ 


-2 („ o ) | » g4 W 

3A 


U)= 0) 


is then explicitely known, and so is 
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3E 


(N) 


£ 44 ( V" 4 33 (N)( "o)| 


3w w=o) 3 10 oj=tij 

0 0 


3oo o)=co 


9A..CN) f -v A * (N) 3e„ 

+ i 44 | e 7 .k) + A .to . i . 33 


3(jo o)=ti) 


33 44 Co) 


3(j0 CO — 05 


0 3A- * (N) f , - A (N) r , 3e Ty1 

- 2 | 34 | £ 34 (^ 0 ) _ 2 A ?,i OO 1 34 


34 


3oj U)=U) 


3u) w=w 


Thus, equation (77) gives an analytic determination of the first order 
error for the natural frequencies in motions of the rolling type. 

It is interesting to note that equation (76) also shows that 
| w n ^ e ^ " w n I increases with n , i.e. the lower modes are obtained 
with higher accuracy for a given number of modes per component. This 
fact is expected, in a way similar to the convergence in a pure Rayleigh 
Ritz analysis. 

This statement is proven as follows: 

(N) 

The A_ J (id) can be written in the following form: 


rNn N c 
- £ Ji 

J n=l 




with D 1 - 1 >0 

n 


2 ij 
<0 -D J 


n 


(N) 


hence A.. (w) are of the order of _J_ for any i, j. and same 

J 2 
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conclusion is valid for the 


C 1 '* 
e. . (w) = N+l 

w 2 -D ^ 
N+l . 


The respective derivatives 


3A 


(N) 

ij 


3 to 


3 e 

and ij are then of the order 

3co 


of 1 . Hence Det ^ (w) and E^ (to) are of the order of 1 , and 


to 


3Det ^ and 3E^ N -* are of the order of 1 . 


3(0 


3(0 


(0 


From equation (77) it follows that [ gd ^ ■ « | is of the order 

of to , thus increasing with go . 

The initial truncation errors and the approximations of equation 
(76) propagate also in the value of the Lagrange multipliers as follows: 
The non-dimensional Lagrange multipliers corresponding to the 
antisymmetrical rolling motions, A^ and A^ are solutions of: 


A 33 Ca ° X 3 + A 34 Cu) A 4 = 0 

a 34 (uj) T 3 A 44 ( w ) J 4 = 0 

By letting A^ = 1 in equation (78) , 


A 


4 


A 


33^ 


(78) 

(79) 


( 80 ) 


and this solution substituded into equation (79) leads to: 
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A 33 (uj) A 44 ^ 


A34O) 


which, for each natural frequency computed, can be written: 


Det (u n ) 


X J = 1 * 2 


4 34 f “n’ 


or, at the first order, 


X * = 1 + 


+ E CN) Coi n) 


4 34 <"„> * 24 34<"„> E 34<“n’ 


It follows that the first order error in the value of X_ for the 

J 

resonant freq . co is : 

^ n 


Det ^ (to ) + E fN )(> ) 

1 JL _jv_ 

“ (N) 

>3 4 34 (“n> * 2 4 34'V e 34<V 


which can be written also by using a one-term binomial expansion 


Det^f^) E- CN) (^) 

\ C ° etfN) (“ n ) 


7 3 4 34 <«n> 


2 e 34 ( “n ) 
(N) 

A 34<V 


Equation (80) gives explicitely the first order error in the value 
of X ^ by similar algebraic manipulations: 
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From equations (58) and (60) it follows that E_ (1) is pro- 

q n 

portional to E_ and E_ (2) is proportional to E_ for each natural 

A_ q A. 

3 M n 4 

frequency (with obvious notations for E_ (1) and E_ (2)), and hence 

% ^n 

the initial errors are finally propagated into the mode shapes of the 
wing and tail as follows: 



Of course, E w d) (0) = E w ( 2 )f0) = 0 

n n 


From equation (57) E_ is explicitely computed: 



(E e (5 ) + e e re ) ) 

T n i T n 2 
A 3 A 4 


p C2) 

n 
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and thus the error in the angular twist of the fuselage for the 
natural frequency is given by: 


E 0 «) = - YjV E 

71 S=1 


(E_ (a) n ) 0 (5j) + E_ Cw n )0 p C? 2 ) ) 0 p (O (83) 
^4 

F C«( U ) 

p n 


The accuracy on the actual physical displacements can be computed 
from equations (29), (30), (31) with c_, c^ and c ^ = 0 (1) 

L *1 *2 

and by using equations (81), (82), (83). 

The error analysis for the motions of pitching type is far more 
complicated as shows equation (72) and it is faster to draw a conver- 
gence graph versus the number of modes used for a given natural fre- 
quency than to obtain analytical expressions for the errors by the 
procedure explained previously. 

Furthermore, an idea of the order of magnitude of errors of the 
type given by equation (76) is given at the first step by the process 
of trial and error itself as follows: 

The corresponding det (io) is plotted versus u and the regions 
where det (w) changes sign are considered. For a case such as the 
one illustrad below 

dchC^) + l fro, n's 


I 

X 
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where det (to) has steep variations in the neighborhood of the 
natural frequency, values such as and are a good approxi- 
mation of w despite the fact that neither det (w. ) nor det (uO 
are close to zero. 

The frequency solutions in such cases are quite accurate, however 
the values of the Lagrange multipliers (and hence the modal shapes) 
are given with less accuracy precisely because of the fact that 
det (wp or det (u> 2 ) are not small enough for this purpose. 

If the opposite case occurs in the neighborhood of a natural 
frequency : 

, . date point's. 



* to 


i.e., det (w) has a locally small slope, even though det (u> o ) and 

det foi ) are very close to zero, w and w* do not bound the resonant 
v o 0 0 

frequency very accurately. However, the mode shape and the informa- 
tion of structural type are given with better precision. 

Situations of the first type tend to occur for the lower natural 

C» 

frequencies and smoother variations of det (w) are seen for the higher 
frequencies and thus the corresponding conclusions can be drawn. 

d) Results and convergence study in a check case 


As in the previous chapter, the theoretical and numerical part 




64 


of the method are first checked with classical results known analyti- 
cally. 

The reference case is taken to be a cantiliver beam which can 
be considered a particular case of the airplane - like configuration 
in the following manner: 

Values of the parameters Physical significance 

^ ^1 = ^2 = ^ both the tail and the wing are placed 

at the middle of the fuselage 

2) Vj-*- 0 for rolling type of the fuselage is given an infinite 
motions inertia in twisting 

^l"* ® ^ or pitching type the fuselage is made infinitely heavy 

of motions 

’^2^ °° ^ or r °Hi n g type the fuselage is made infinetely rigid 

of motions in twisting 

00 for pitching type the fuselage is made infinitely rigid 

of motions in bending 

The set of parameters in ( 2 ) is such that the fuselage is con- 
strained to stay still, i.e. clamp the wing and the tail in both kinds 

of motions. In fact, mathematically one could have allowed and ij, 

1 v 2 

to be finite. 

For the sake of simplicity the two beams representing the wing and 
tail are also made identical by taking 

v 2 ~ y 2 = ^4 = 1 for tfle antisymmetric motions 
v 2 = v 3 = ^4 =y 2 = 1 for the symmetric motions. 
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and hence the frequency equation is : 

(84) 

i.e. the frequency equation in the case of a beam clamped in its 
middle, in an analysis using only antisymmetric free-free modes 
(see Appendix III). 

Of course, all solutions to of equation (84) are in fact double roots 
as we have two identical cantilever beams in vibration. Note that in 
the analysis developed in part b) the non-dimensionalization has been 
made using half-length of the wing, hence there is no adjusting factor 
to be computed for the non-dimensional frequencies. Figure 13 shows 
the rate of convergence of the first natural frequency towards that 
of a cantilever beam versus the number of antisymmetric modes used; 
the convergence is very fast for the first few modes. Furthermore 
the shapes of the first three modes have been obtained using ten anti- 
symmetric modes in the analysis and they are plotted in Figure 14, 

15, 16. The frequency bounds are given by the analyses developed in 
part c); the percentage errors are almost equal for the three natural 
frequencies (4.1, 4.2 and 4.5 respectively). The values of the slopes 
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at the origin is given as an indication of the modal shape accuracy. 
Vibrations of the pitching type 

With the limiting values for the parameters given previously, 


A 


11 



z 6 2 CO) 
1 “ 



n 


12 


= A 


IS 


* A, . = A__ = A„ = A 


16 


25 


26 


56 


0 


55 



Z 

n=l 


e n 2 C 0 ) 

2 , ,.2 2 

(l> - C n "l) Tf 


The frequency equation becomes in this case: 
1 

DO OO 

C 2 2 ( 0 ) ) 2 c 2 e 2 co) ) 2 = o 

n=1 ~1 4 n=1 -1 2~2 

« - V w 2 - (n-l)V 

4 

- - — ... . 

The first factor corresponds to the bending modal behavior of 
a beam clamped in its middle, the second factor corresponds to the 
cantiliver beam in twisting (see Appendix III) . 

As expected, the bending and twisting motions are decoupled 
from each other. The convergence rates towards the first clamped- 
free natural frequency respectively for bending and twisting are 
both shown in Figure 17 so that they can be analysed comparatively. 

The percentage error in the first bending frequency computed 




- 67 - 

using symmetrical modes only decreases very fast to zero in fact much 
faster than in the case in which only odd modes are used. This result 
is to be expected as the natural frequencies of the even modes are 
lower than the ones of the odd modes . 

By comparison, there is slower convergence in the torsional 
modes and Figured which represents a more extensive graph of the 
twisting convergence shows that in order to get 0.9% error in the 
non-dimensional frequency (precision attained with only 3 symmetric 
modes in bending), 21 torsional modes have to be included in the 
analysis . 

This fact is surprising as there is no physical reason for 
which it could be expected, however, it can be justified on mathe- 
matical basis. 

Consider the finite N-termes frequency equations: 


PCco) 


E 

n=l 



(0) 



- 0 


(85) 


and 


0 » 


E 

n=l 


e n 2 C0) 

2 ( 1 1 2 2 
0 ) - (n-1) tt 


= 0 


( 86 ) 


and the first natural frequency corresponding solutions and Qj . 

The problem is that of determining how much will u)j and change 
by considering one more term respectively in the analysis, or equiva- 
lently, how much will P(w) and Q(w) respectively be modified by use 
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of the additional term. 

By separating the rigid body contributions in bending and 
we can write P(co) and QC<o) as follows: 
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The N = +1 terms in each of the equations (87) and (88) 
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_! > _2 

2,4 2 2 2 

o) - k N+2 w - (N+l) 7T 

By use of this result together with inequality (89) it follows: 


♦ N ’ 2 (0 > 

2 . 4 

“ ' k N,2 


1 

~~2 , M ,.2 2 
U> - (N + l) TT 


(90) 


This inequality shows that an additional bending mode in the 
analysis provides more contribution (in fact much more contribution) 
to P (w) than an additional twisting mode in Q(u>) . 

Thus the value of will be more significantly modified by an 
additional mode than the value of 

This explains the relatively slow convergence rate of the tor- 
sional frequency by comparison with that of the bending frequency. 

Furthermore, this result suggests that for the next investigations 
more modes should be used in twisting than in bending (in fact quite 
a large number of modes) . This fact does not present any noticable 
inconvenience however, as it corresponds just to performing more 
additions; from the point of view of computer work, the additional 
time and region necessary are negligible. 

4 . Experimental work 
a) Introduction 

The idea of a completely unconstrained structure presents no con- 
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ceptual or theoretical difficulties, however, practical problems 
arise when such a structure is to be simulated in the laboratory. 

For experimental purposes the specimen has to be supported, and 
thus one of the first problems which is .to be solved in the design 
of the experiment is that of finding a suitable type of supports 
which characteristics will not affect the resonant frequencies of 
the structure. In particular, for the present analysis care has 
to be taken that neither the natural frequencies of motions of 
rolling type, nor the frequencies in the pitching motions will be 
changed by the type of supports used. 

The next subparagraphs describe the experimental model, the 
type of supports which simulate its free motion and the experimental 
setting. 

b) Experimental model 

The beam-like structure on which the experiments have been per- 
formed is shown in Figure 19. The components are aluminum alloys 
beams of rectangular cross-section and following characteristics: 


length 

width 

thickness 

E 


beams (1) and (2) 

10 in 

0.5 in 

0.1 in 

10 7 lb/in 2 


v 0.33 

p 0.1/385 lb/in 3 /in/sec 2 


beam (3) 

10 in 
0.25 in 
0.05 in 
10 7 lb/in 2 
0.33 

0.1/385 lb/in 2 / in/sec 2 


The values for E, v, p are taken from Ref (40). 
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Beams (1) and (2) are mounted perpendicular to each other on 
their middle, beam (3) is mounted parallel to beam (2} at 3 in. dis- 
tance from i t . 

Small steel coupons are mounted at :one of the tips of each com- 
ponent so that the structure can be excited at various tip locations 
by a magnetic transducer, 
c) Choice of the optimum supports 

The experimental model with the two identical supports is shown 
in Figure 20. The elastic supports at both tips of beam (1) have been 
perfectly alligned so that they do not affect the rolling of this 
beam. The constraining effects on the bending motions of the beam 
have been minimized by a proper choice of the spring constant, proceed- 
ing as follows: the idea is to choose the springs soft enough so that 
the frequency of the rigid body translation of the model resting on 
the two springs is much lower than the lowest free-free bending fre- 
quency of beam (1) . 

This lowest frequency is found to be of the order of 200 cps . 

and with the choice of a spring constant equal to 0.5 lb/in, 

the frequency of the corresponding rigid body translation is of 

9.5 cns. The value of the spring constant is related to that of 

the frequency of the rigid body translation by 

2 

k = Ma> , where 
2 

M is the mass of the specimen 

k is the spring constant 

u> is the dimensional frequency (rad/sec) and the factor 1 

2 
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corresponds to the fact that the two springs are identical, 
d) Experimental setting 

The set up of the experiment is shown in Figure 21 and more 
detailed explanation is given in the sketch which follows: 
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The oscilloscope has a double beam feature and a polaroid 
camera connected to it, so that the input and output signals can 
be compared and picture information recorded. 

Two types of resonances were found as shown in Figure 22, and 
labeled type a and B, 

The lower part of the Figure 22 shows a resonance of type a: 
the output signal (lower beam) is seen to be at -90° phase with the 
input signal. 

The upper part of the Figure 22 shows a resonance of type 3: 
the output signal on the lower beam now precedes the input signal 
by 90°. 


e) Frequencies measurements 

Three sets of measurements have been performed with the follow- 
ing positions for the driver and pick up: 

. set of measurements (1) 

Both the driver and pick up are at the tips of beam (1) ; the 
driver is positioned farthest away from beam (3), the pick up is 
placed at the tip closest to beam (3) . With this arrangement the 
following resonant frequencies were detected in increasing value: 


62 cps 
112 cps 
361.1 cps 
480 cps 
575 cps 


weak signal type a 
weak signal type a 
strong signal type a 
strong signal type a 
weak signal type a 
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. set of measurements _ (2) 

The arrangement is as shown in Figure 21; the driver and pick 
up are at each of the tips of beam (2) . The measured resonant 
frequencies are: 


32 . 5 cps 

weak signal 

type 8 

62 

cps 

weak signal 

type a 

110 

cps 

weak signal 

type a 

264 

cps 

weak signal 

type 8 

489 

cps 

strong signal 

type 8 

575 

cps 

weak signal 

type a 

578.7 cps 

strong signal 

type 8 


. set of measurements _( 3) : 

The arrangement is similar to the previous case, however, the 
driver and pick up are positioned now at the tips of beam (3) . 
Because of the larger flexibility of this beam, the resonant frequen- 
cies were easier to detect and in fact this set of measurements is 
considered to be the most satisfactory of the three. The following 
resonant frequencies were detected: 


.32.5 cps 

strong signal 

type 8 

62 

cps 

strong signal 

type a 

112 

-cps 

weak signal 

type a 

264 

cps 

strong signal 

type 8 

361 . 1 cps 

strong signal 

type a 

480 

cps 

strong signal 

type a 

489 

cps 

weak signal 

type 8 
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57S cps strong signal type a 

578.7 cps strong signal type 6 

As it can be seen, many natural frequencies were found in more 
than one set of measurements; very often their value agree axactly 
from one measurement to the other, for the few frequencies for 
which this is not the case the data corresponding to the strongest 
signal on the screen of the oscilloscope have been retained. 

Tabulated below is a summary of the measured resonant fre- 


quencies , 

the way in which they have been detected and 

their type 

natural 

driver' s 


pick up's 




frequency 

position 


position 



type 

32.5 

tip of beam 

(3) 

other tip 

of beam 

(3) 

B 

62 

it 


n 



a 

112 

tip of any beam 

other tip 

of the 


a 




same beam 




264.3 

tip of beam 

(3) 

other tip 

of beam 

(3) 

3 

361.1 

M 


IT 



a 

480 

tip of beam 

CD 

other tip 

of the 


a 


or (3) 


same beam 




489 

tip of beam 

(2) 

other tip 

of beam 

(2) 

B 

575 

tip of beam 

(3) 

other tip 

of beam 

(3) 

‘ a 

578.7 

tip of beam 

(2) 

other tip 

of the 


3 


or (3) 


same beam 





The resonant frequencies of this beam-like model have also 
been investigated theoretically, and the next paragraph deals with 
the numerical results which have been obtained and the conclusions 
which can be drawn from the comparison with the experimental work. 
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5 . Numerical analysis for the experimental model: theoretical results, 

comparison with the experimental results 

We remind the reader at this point that although the numerical 
results are presented with the frequencies in increasing order for 
easier comparison with experimental data, the frequencies of rolling 
type and pitching type of motions have been obtained by independent 
computations as discussed in paragraph 3. This factor will be taken 
into account however, in the quantitative analysis. 

The lowest computed natural frequency is = 36.7 cps and the 
physical deformation of the structure is shown in Figure 23. In this 
mode beam (1) undergoes only angular twisting; beams (2) and (3) 
bend antisymmetrically, predominantly in the rigid mode. The value of 
this resonant frequency can be related to that of 32.5 cps. measured 
experimentally; as it can be seen in Figure 23 this resonance could 
not be detected in the first set of measurements. A comparison bet- 
ween the amplitudes of the rigid body rotations of beams (2) and (3) 
shows that reason for which a stronger response has been obtained where 
the driver and pick up were both placed at the tips of beam (3) . 

The value of the second calculated natural frequency is im ^=69 cps; 
as shown in Figure 24, the modal behavior is of the pitching type: 
beam (1) rotates as a rigid body, beam (2) undergoes weak bending in 
the first symmetric mode, and higher amplitude bending of the same 
shape (but opposite sign) is observed for beam (3) . 

The corresponding second resonance measured experimentally 
occured at a driving frequency of 62 cps (with the strongest response 
in the third* set of measurements); it is interesting to notice that 



this resonance has been found experimentally to be of different type 
than the previous, specifically of the type a. It can be anticipated 
on these basis that resonance of the type a will correspond to motions 
with symmetric behavior of beams (2) and (3) while type & resonances 
will be characteristic of the "rolling motions". In the third mode 
(for which the computed value of the frequency is = 125.6 cps) the 
deformation of each component is similar in shape and amplitude (see 
Figure 25); this resonance has been detected in all three sets of mea- 
surements, and found to occur at a frequency of 112 cps. Note that 
indeed it corresponds to a type a deformation. 

The next higher natural frequency has a theoretical value of 295.5 
cps; Figure 26 shows the deformation of each component and it can be 
seen that in this second rolling mode both beams (2) and (3) do show 
an elastic behavior. The corresponding measured value of the frequen- 
cy is of 264.3 with strongest response picked up on beam (3); indeed 
it is a resonance of type 6. 

The next four natural frequencies obtained theoretically and the 
corresponding modes are shown in Figures 27, 28, 29, 30, 31, respectively; 
similar analysis show qualitative checking between theoretical and 
experimental results. As the comparison shows good qualitative agree- 
ment, it is of interest to further investigate in a quantitative 
manner the agreement (or discrepancies) between theoretical and ex- 
perimental results. For the time being the quantitative analysis will 
be done separately for the modes of rolling type and pitching type, as 
the theoretical results have been obtained by independent computations. 
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The results are given in tabulated form for clarity, and shown below: 

Modes of rolling type: 


01 

n 

theoretical 

(cps) 

(l) 

n 

experimental 

(cps) 

percentage error 
normalized by the experi- 
mental data (%) 

36.7 

32.5 

11.5 

295.5 

264.3 

11.8 

546.6 

489 

11.6 

648.2 

578.7 

12.0 


At first view the percentage error seems to be too high for 
being satisfactory; further analysis shows, however, the interesting 
fact that this error is constant for the four frequencies. This remark 


suggests that possibly part of the error comes from the differ- 
ence between the exact and assumed characteristics of the experi- 
mental model. If this is true then better agreement should be 
provided by comparing the two sets of frequencies in a non- 
dimensional manner. This can be achieved by comparing the ratios 


Cm /wj theoretical and (u AO experimental as shown in the table: 
Tl 1 n i 

^ ^ ((j ) Percentage error normalized 

,, “ T , by experimental data 

1 theoretical 1 experimental 

1 1 


8.05 


7.93 


1.6 


14.89 


14.67 


1.5 


17.66 


17.36 


1.7 



80 - 


Column 3 illustrates that the previous assumption was correct, but 
at least two other possible facts have to be taken into account: 
one is that the components were not welded in the structure but 
glued together and hence the elastic properties of the glue play 
a role at the joints, the second fact is that there is a stiffening 
effect due to the steel coupons, especially on beam (3). It is 
interesting to notice that, probably because of the facts mentioned 
previously, all computed values of the frequencies are higher than 
the measured ones. But, of course, this is also inherent to the 
Rayleigh Ritz procedure per se. 

A simular analysis is now developed for the modes with symmetric 
behavior. 


“n theoretical 
Ccpsj 


Modes of pitching type : 


“n experimental 
(cps) 


Percentage error normalized 
by the experimental data 
(%) 


69 

62 

11.2 

125.6 

112 

12.2 

404.8 

361.1 

12.1 

536.2 

480 

11.7 

643.3 

575 

11.9 


Consideration of the percentage error leads to conclusions similar 


to that of the previous case and thus we develop again a non-dimensional 
manner of comparing the results, as shown below: 
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to 

: n ) 

W 1 theoretical 

to 

t-io 

1 experimental 

Percentage error normalized 
by the experimental data 

1 

1 

- 

1.82 

1.79 

1.7 

5.87 

5.78 

0.9 

7.77 

7.68 

1.2 

9.32 

9.20 

1.32 


Column 3 shows satisfactory agreement between the non-dimensional 
frequencies obtained theoretically and experimentally. Furthermore 
it is of interest to notice that even for different types of modes 
the percentage errors have almost the same value; this suggests that 
the analysis performed for the two classes of frequencies applies in 
fact globally which is an indication of succesfull agreement between 
theory and experiments . 

The global analysis is presented below in tabulated form: 


< jO 

n 

theoretal 

Ceps) 

03 

n 

experimental 

(cps) 

percentage 

error 

(%) 

to 

to," ^ ( 

theoretical 

_IL. ) 

“l 

experimental 

percentage 

error 

36.7 

32.5 

11.5 

1 

1 


69 

62 

11.2 

1.88 

1.86 

l.l 

125.6 

112 

12.1 

3.42 

3.36 

1.8 

295.5 

264.3 

11.8 

8.05 

7.93 

1.6 

404.8 

361.1 

12.1 

11.03 

10.83 

1.9 

536.2 

480 

11.7 

14.61 

14.40 

1.4 

546.6 

489 

11.6 

14.89 ' 

14.67 

1.5 

643.3 

575 

11.9 

17.53 

17.25 

1.6 

648.2 

578.7 

12.0 

17.66 

17.36 

1.7 
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For all of the nine resonant frequencies there is less than 
2% disagreement between theory and measurements when compared in a 
non-dimensional manner. 

It is important to mention that because of the choice of very 
soft springs in supporting the model, the frequency of the rigid 
body translation of the structure resting on springs could be very 
easily filtered out in the experiment and was not accounted for in 
the previous discussion. 

Even if this had not been the case, a valid comparison between 
the two types of data should still be possible, and for this purpose 
the author has developed in Appendix V a theoretical analysis allow- 
ing for the effects of the two springs . 

6. A modal analysis of the Lockheed model L-100 Hercules 
a) Presentation of the model 

Following the production of early C-130 military Hercules air- 
craft, Lockheed - Georgia had decided to offer a commercial version 
of this heavy transport, and the model L-100 is one of the six models 
which have been offered for that purpose. 

Customers of this aircraft have included Alaska Airlines, Delta 
Air Lines, Interior Airways and Pacific Western Airlines. 

A picture of the airplane is shown in the upper part of Figure 
32 and a general arrangement drawing is shown below. 

More detailed information can be found in References (32) and 




q 


(33). 
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b) Numerical results 

The theoretical analysis developed in paragraph 3 has been applied 
with the following values of the parameters representing the model. 


Parameter 


K 

2 


V 


1 


u 


Y 

Y 
v 


2 

1 

2 

1 


v 

V 

V 

* 


2 

3 

1 

2 

3 

4 

5 


Physical significance Value for 

L-100 

point along the fuselage connected with 0 

the wing 

point along the fuselage connected with 1 

the wing 

ratio (mass wing/mass fuselage) - engines not 
taken into account 1 

ratio (mass wing/mass tail) 38.97 

ratio (wing span/fuselage span) 1.36 

ratio (wing span/tail span) 2.52 

mass to inertia ratio between wing and 91.4 

fuselage 

ratio between mass and inertia of the wing 166.84 

mass to inertia ratio between wing and tail 2669.44 
frequency factor of the fuselage in bending 65.47 

frequency factor of the fuselage in torsion 81.05 

frequency factor of the wing in torsion 245.67 

frequency factor of the tail in bending 39.98 

frequency factor of the tail in twisting 1560 


The first twelve natural frequencies of the aircraft have been 
obtained and their values are listed below. A total of 53 modes 
(38 bending modes and 15 twisting modes) has been considered for the 
motions of the rolling type; 125 modes in total (50 twisting modes 
and 75 bending modes) have been used in the analysis of the modes 
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of pitching type in order to compensate for the eventual loss in 
accuracy in the matrix manipulations especially when they are 
nearly singular. 


non-dimensional 

frequency 

type of 

frequency u> 

ft (cps) 

modal behavior 

4.20 

18 

pitching 

8.30 

36 

rolling 

12.03 

53 

rolling 

16.78 

73 

rolling 

18.75 

82 

pitching 

21.65 

95 

pitching 

23.5 

103 

pitching 

24.3 

106 

rolling 

25.14 

110 

pitching 

40.09 

176 

rolling 

40.6 

179 

pitching 

45.17 

199 

rolling 

52.55 

231 

rolling 


The deformation of the structure in the first mode is shown 
in the Figure 33; this mode appears to be mainly a wing mode going 
only small rigid body translations. In fact the natural frequency 
of this mode is close to the first free-free bending frequency 
of the wing alone, more precisely 

W 1 beam cantilivered ^1 < W 1 free free beam 
in the middle 
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The figure shows that the wing translates with the fuselage but 
mainly bends in the first even mode; the ratio of the corresponding 
generalized coordinates (bending to translation) which is an indication 
of the relative contribution of the initial modes, is of 1.6 in absolute 
value. Structural information obtained from the values of the Lagrange 
Multipliers is given below (up to multiplication). 

1 lb 

(taken as reference) 
0.019 lb 

0.24 lb x in 

0.011 lb "x in 

As expected, the shear force developed at the point where fuse- 
lage and wing join, is predominant in this mode. 

The structural deformation of the aircraft in the second mode 
is shown in Figure 34. This mode is characterised by an antisymmetric 
behavior of the wing and the tail and only angular deformation of the 
fuselage. This angular deformation is composed of rigid twisting, 

j . 


Shear force developed at the point of attachment 

between : 

the wing and the fuselage 

Shear force developed at the point of attachment 

between : 

the tail and the fuselage 

Torque developed at the point of attachment 
between the : 

wing and the fuselage 

Torque developed at the point of attachment 
between the : 

tail and the fuselage 
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with contributions from the second and first elastic torsional modes 
by order of importance. For purpose of qualitative information: 




3.9 


The wing rotates with the fuselage and bends in the first anti- 
symmetric mode; 


| q ^ bending | = 1.3 

q ^ rotation 


The tail shows practically no elastic behavior. 

The value of the torques per unit length developed at the points 
of connection between fuselage and tail and between fuselage and wing 
respectively is given by the ratio 


= 0.15 


The next higher mode presents again antisymmetric behavior of the wing 
and tail as shows Figure 35, in fact the gross-deformation of the 
structure is very similar to that of the previous mode except for the 
fact that in this mode the angular twisting of the fuselage has same 
direction at the points of connection with the wing and the tail. 

Detailed analysis shows that in this mode the behavior of the 
fuselage is predominantly elastic, in fact mostly in the first elastic 
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and very little 

. 1 ) 

contribution from the 4th elastic mode 

f 2 

( r — - 32 .i ) 

p 

5 

As previously, the wing bends predominantly in the first antisymmetric 
mode, but undergoes less rotation than previously and shows also 
influences due to the second free-free antisymmetric mode. 


a ^bending | = 4.92 

| n ^bending j = 

a ^ rotation 

a ^bending 2nd mode 


as it can be seen in the lower part of Figure 35 the tail acts mainly 

rigid, however, slight elastic behavior can be graphically detected 

(quantitatively, | a ^rotation | = 18.3 ) 

(21 

— 'bending 
a 

The ratio of the torques per unit length developed at the wing and 
tail is now equal to 0.66 as compared to 0.15 in the previous case. 

The modal behavior for the fourth natural frequency is shown in 
Figure 36 and similar discussion can be developed for this mode. 

For a more complex analysis taking into account the masses and 
inertias of the engines (represented by rigid structures attached to 
the wing), the author refers the readers to the theoretical develop- 
ment of Appendix VI . 

The present model does not provide information about the chord- 
wise bending occuring in the wing and the tail because of the 


mode with some rigid contribution P 
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one -dimensional assumptions for these components. A more sophisticated 
analysis taking these additional effects into account, is provided in 
the Appendix IV where the wing and the tail are represented by plate 
elements with uniform properties. 
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CHAPTER IV 

MODAL BEHAVIOR OF A HIGH ASPECT RATIO AIRPLANE 
WITH OBLIQUE WING .AND TAIL. 

1 . Statement of the problem and basic assumptions 

The present chapter may be considered as a generalization of 
Chapter III. This analysis has not only a pure academic interest 
in showing the generality of the method, but also a practical inter- 
est. Since the oblique wing configuration was proposed, there has 
been an increasing concern over the aeroelastic problems arising 
from such an uncommon geometry. As most of the flutter analysis 
are based on the structural normal modes, the purpose of the present 
chapter is to provide precisely these normal modes and frequencies. 

It is of interest to note that the assumption of a high aspect ratio 
wing is particularly suitable for such a type of configuration as 
it can be seen in Reference (34) and Figure 37; the additional assumptions 
under which the problem is treated are similar to those of Chapter III 
and are reminded below: 

1) the airplane structure is considered to be of beam type 

2) the displacements of the beam components are governed by the 
Bemouilli-Euler equation 

3) the geometrical and elastic characteristics are constant along 
the span for each component. 

The geometry of the structure is illustrated in the top view 


given below: 




-L 


The angle a, angle between the oblique components and the fuselage, is 
the main new parameter of the problem. 


2 . Rigid body dynamics in the linearized theory 
Notations : 


R : 
o 


G: 


3 (G,S): 
M : 

fi(S/R ) : 
o 


Gallilean frame of reference 

Center of gravity of the structure 

Classical Euler angles as defined in Chapter III 

Inertia tensor of the structure at the center of gravity 

and with respect to its own system of axis 

Total mass of the structure 

Vector defining the rotations between the frame of the 
structure and R . 


o 

f (G/R : Acceleration vector as defined in Chapter III. 

o 

H ((S)/R ): Angular momentum as defined in Chapter III. 
o 

In absence of exterior forces Newton's law yields: 


M * r(G/R Q ) = ($ 


( 91 ) 


This relation decouples the translations in the three directions 
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from the rotation motion around the center of gravity. 

For zero exterior torque the angular momentum equation is simply: 

H ( S/R ) = 0 (92) 

dt o 

which can be also written as: 

^ (J(G,S) . 3(S/R 0 ) )= o (93) 

Due to the absence of symetry, the tensor 3 (G,S) has a non-diagonal 
representation and can be written in the following matrix form: 

r IjOO Jj(a) 9 

2(G,S) = Jj(a) 1 2 (a) 0 

0 0 I_ 

l 3 

where (tt/2) = 0. 

Hence the linearized equations of the rigid body motion are: 


(I x (a) 

cos a - J^a) sina)0 

+ (Ij(ct) 

sina + J^(a) cosa)t() 

= 0 (94) 

(Jj(a) 

cosa - 12 (a) sina)0 

+ (J a (a) 

sina + IjCa) cosa) 

= O' (95) 

V - 

0 



(96) 


This set of equations shows that the yawing motion decouples from 
roll and pitch, but the later two remain coupled unless a = rr/2 

This fact introduces a fundamental difference from the types of 
structures investigated in the previous chapter: indeed, if previously 

the modal analysis could be done from the onset separately for symmetric 
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and antisymmetric modes, this is no longer the case for the present 
problem as it can be seen from the rigid body analysis. 

3. Elastic modal behavior, derivation of the frequency equation 

a) Statement of the elastic problem: degrees of freedom, derivation 
of the constraint equations. 

The motions investigated are of roll-pitch type; the bending vi- 
brations in the y-direction (yaw-type of motion) are not considered 
in the present chapter. 

Each one of the components is allowed to bend and twist rigidly 
and elastically and the corresponding modal displacements are sought 
in the following forms: 

- for the fuselage: - L <_ x L 

bending displacement: 

W (x,t) = S 0(t)6 (x) (97) 

n=l ™ n 

twisting angular displacement: 

0(x,t) = l P (t) 0 (x) (98) 

n=l n n 

- for the oblique wing: - 1^ <_ x^ <_ 1^ 

bending displacement: 

GO flV 

w 1 (x 1 ,t) - I (t) ^(Xj) (99) 

n=l 

twisting angular displacement: 

0 lf x l-O ■ 1 . p " ) Ct) e/xp 


( 100 ) 
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- for the oblique tail: - ^ — x i — ^2 

bending displacement: 

w 2 ( Xl ,t) = I q^ 2) (t) * n (xp C101) 

n=l 

twisting angular displacement: 

g 2 (x 1>t ) = z ^ p^ 2) (t) e n (x 1 ) C102) 

As explained in the previous chapter, the physical displacement of 
each component is rigorously the sum of the bending displacement 
and an additional term due to the twisting of that element; however 
in agreement with the one-dimensional assumption it is reasonable to 
neglect the contribution due to twisting (based on the fact that the 
ratio of the chord to the length of the element is considered small). 
The constraint conditions insuring the equality of the physical dis- 
placements at the points of connection between the components can be 
written then . 


W(5j,t) = w^ (0, t) 

(103) 

W(S 2 ,t) = w 2 (0,t) 

(104) 


after having suitably non-dimensionali zed W , and W 2 so that their 
domain of definition be [-1, 1], 

Another set of constraint conditions is derived from enforcing 
the equality of the slopes of the components in the x and y directions 
and leads to the equations 

©(£, » t) = e.(0, t) cosct + 1_ 1 (o.t) sina 


(105) 
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1 3 W (5 , t) = - 0 ( 0 , t) sina + _1 _1_ (0,t) cosa (106) 

L 3C 1,3? 

1 1 

for the wing and fuselage, and 


0(£,» t) = 0(0, t) cosa + 1 2 (0,t) sina (107) 

1 2 

w 

1 3W. (£ ,t) = - 0 (0,t) sina +1 2 (0,t) cosa (108) 

L 3C 2 1 2 3 ?1 

for the tail and fuselage'. 

The slope condition of Chapter III can be found by setting 
a = tt/2. 


b) Derivation of the frequency equation 

The total kinetic energy of the structure is 

L • X 1 h 

T = 1, ( J W 2 ( x ,t) dm + 1 w, (x ,t) dm + { w 2 (x ,t) dm- + 

2 -L -1 1 1 1 1 ~ l 2 1 1 

L .2 1 .2 l 2 .2 

! l I0 (x,t) dx 1^0 (X x ,t) dx : +_{ r 2 0 (Xj , t) d Xj ) (109) 

where T, t^, I are the respective mass moments of inertia per unit 
length about the torsional axis of the fuselage, wing and tail. Their 
integral expressions are given in the paragraph 3 section b) of the 
previous chapter. 

Substitution of equations (97) - (102) into equation (109) yields 


T 


. 05 

1 1 

2 n=l 


.2 

(MQ^t) + Mj 



2 

(t) + m 2 



2 

(t) + 
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• 2 .( l ) 2 • ( 2) 2 

IP n Ct) + I l P n Ct) + *2 P n C0) (ll0) 


by virtue of the orthonormality of the free-free component modes. 

The total potential energy of the structure in bending and torsion 
is : 


L 

U = 1 (El / 

2 ' L 


(3JV 

3x 2 


(Xjt))' 


dx + El, 


[h (x. t)) 2 dx + 

L 2 1 

3X. 


1, a 2 w 9 1 

EI 2 -[ ( —1 ( x i' t5 ) dx i + GJ _l 
2 3x 2 


( 3£ (x,t)) i ’dx + 
3x 


CGJ) { *( ^1 (x t) ) 2 dx + (GJ), J ) 2 ( ^2 (x ,t)) 2 dx ) (111) 

1 3Xj 1 1 2 3x 1 1 

where I, 1^, I are the structural moments of inertia about the neutral 
axis of the fuselage, wing and tail respectively, E and G are the Young 
and shear modulus of elasticity, and (GJ) , (GJ) ^ and (GJ)^ are the 
respectives torsional stiffnesses of each of the components. 

Equation (111) can be rewritten by use of (97) - (102) and the 
free-free beam equations in bending and torsion in the following way: 

U = 1 f (El k 4 Q 2 (t) + (EI) 1 k 4 q a) (t) + 

T t,-i “ 7 n — 7~ n M n v 

2 n=i 3 1 3 

L l 1 

CEI) 2 k 4 ) (t) + GJ (n-1) 2 -n 2 P 2 (t) + 

3 n n ~ 4 n 

2 
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< GJ >1 (n-1) 2 tt 2 P (1) (t) + CGJ) 2 Qi- 1) 2 * 2 P,? 2:) Ct) ) (112) 
ll 4 1 2 4 

Using again the Rayleigh Ritz expansions given in the equations 
(97) - (102) for substitution into the constraint conditions (103) - 
(108), the following modal expressions are obtained: 


5 ~ (t)* n (0) ) = 0 

n=l 


(113) 


? - q,S 23 ct) * n (0) ) = 0 

n= 1 


(114) 


( 1 ) 


5 , W>W - P n 11 ' t ) V°> CM “ l - t — n 

n=i 


- V^j,' (0) s ino) = o (115) 


(1) 


? (V!l <f>’_ (5,) + pf 1 ^(t)6 (O)sina - \ (t) <f>' (O)cosa) = 0 (116) 
* ni n n n 


n=l L 


,( 2 ) 


( 2 ) 

- 


E ' n ; (t) 0^ (COcosa - (0): sina) = 0 (117) 

n=l 1 ^ 

00 0 (t) (2) a ^ 2 ^ ft) 

E ( V J t'ttJ * P n (t)0_ (0)sina - q n ^ V_(0)cosa) = 0 (118) 
n=l L n 1 n n — n 


These equations weighted by their respective Lagrange multipliers 
contribute to the formulation of the Lagrangian for the entire elastic 
structure: 

flfct) = T (t) - U(t) + X x (t) x eq. (113) + * 2 (t) x eq.(114) + 


X 3 (t) x eq. (115) + * 4 (t) x eq.(116) + A g (t) x eq.(117) + 



As mentioned in Chapter III the values of the Lagrange multipliers 
give interesting structural information: (A^ = 12 relat ^ ve t0 

displacement conditions give a comparative value of the shear force 
developed at the points of connection between the fuselage and wing 
and the fuselage and tail; (A^)^ - 3 6 give a relative value of 

the bending moments about the x and y axis developed at the points 
of junction between the components. 

Moreover the values of the Lagrange multipliers indicate which 
type of behavior is predominant in a particular mode before even having 
actually determined the mode shape. This statement is a little of an 
anticipation at this point but its proof is clear from equations (125) - 
(130). 

The equations of motion, known as Lagrange's equations for the 
system are: 

♦ MQ^(t) + EI_ k^ QJt) - V tH n C? l } ' X 2 (t)<!> n (5 2 ) " 

L 3 

- V^i<t' n (s 2 ) - o dip) 

L L 

. ILP n (t) + GJ (n-l)V P n (t)(t) - X 3 (t)e n (e i ) - 
L 4 ' 

8„(5 2 )-» 


( 120 ) 
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. M^(t) + (EI) 1 k* q^^Ct) + X 1 (t)cf, n CO) + 

h* 

X 3 (t) <j)' (0] sina + X 4 W <f> r (0)cosa = 0 (121 

1 n 1 n 
l l X 1 

• I.l.P^V) + CGJ) 1 ■ (n-nVpf^Ct) + A,(t)e (0) cos'o - 

l 1 n — : * -r n 3 n 

H 4 

A 4 (t)6 n (0)sina = 0 (122) 

• • (2) fEIl ^ (2) ^ i>'> 

• M 2 % Ct) + 2 k n % W + x 2 (t) *n C0) + — ^ti (0) 

1 2 3 h 

sina + X 6^ <(>' (0) cosa = 0 (123) 

X 2 

• I 2 1 2^i| C )(t) + (GJ)2 ( n -l) 2 ^ 2 p n C2:, Ct) + X 5 (t) 6 (0) 

*2 4 1 2 


cosa - X^(t) 6^(0) sina = 0 (124) 

together with the constraint equations (113) - (118). 

The harmonic behavior in the modal analysis is taken in a similar 
form to that of the previous Chapter and is reminded here: 
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M, 1, M L . ML 

v, = 1 1 ; v„ = 1 , v, = 1 

l T i j o i 

1 12 


*1 = 


SI 

ML 3 

(EI) ] 

V? 


GJ 


2 = 


IL 


(El). 


M 1 ‘ 
11 




(GJ). 


(El). 


M 1 
11 


(El). 


4 = 


M 2 1 2 ' 

(El). 


5 = 


(GJ) 2 

15 

(EI) , 


M 1 

11 


M iY 


2 , , 4 

p (1) = o) - 


n 


1 n 


P (2) 2 , r ..2 2 
F ^ } = (i) - (n-1) it 

Tl 2 -r 

4 


c (3) 2 - 4 

F ' = tu - k 

n n 


r (4) 2 , n 2 2 

F n v = a> - <J* 3 (n-lj £ 

4 


F < S > » „ 2 - Kk 4 
n T 4 n 
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v ( 6 ) 2 
F = (1) 
n 


- ip- (n-l)‘ 


— 

4 


In terms of these parameters and of the non-dimensional Lagrange 
multipliers the generalized coordinates are expressed as follows: 



C*lW + X 2*n«2 )+ Vl*'n (E l )+ Vl*Vy) 

CD 

F 


(125) 


P 

n 


Y l v l 


(hW + x sW) 

( 2 ) 


F 

n 


(126) 



A <j> (0)+A (J)' (0)sina+A ' (0)cosa 
x ti 5 n 4 n 


(3) 

F 

n 


(127) 



v„ ( A 3 C0Sa " X 4 Sina )0 (0) 
(4) " 


n 


(128) 



( A 2^n^ +A 5 Y 2^'n^ sinot+A 6 Y 2 < * ,, n^ COSa ) 

(5) 


( 129 ) 
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P 


n 


( 2 ) 


y 2 V 3 


,A._cosa - X.sina 
{ 5 O 

( 6 ) 

F 

n 


•V 0 ’ 


(130) 


It can be noticed that only even modes contribute to the twisting 
of the wing and tail as it can be seen from equations (128) and (130). 

It is also interesting to point out by the use of two examples 
how the modal behavior of the structure can be predicted from the 
value of the Lagrange multipliers directly. 

- Example 1: 

Assume that for a particular mode the values of the Lagrange multi- 
pliers are such that 

T. 

j — | <<1 for any i = 2, . . .6 
A 1 

Then equations (94) - (99) can be written at zeroth order approximation. 


Qh = - Vl 




CD 


n 


P 

n 


= 0 


5 
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n 


for any n . 

By substituting these equations into equations (97) - (102) 
we obtain the expressions of the modal displacements at the lowest 
order: 


W(5) 


2 

n-1 F (D 


(131) 


WjCO 


2 +n (0) +n (O 

n=l F (3) 
n 


(132) 


e(0 - Bj(c) = w 2 (?) = © 2 (c) = 0 (133) 

The modal displacements being defined up to a multiplicative 
constant, division by X^ is justified. Equation (133) shows that 
the tail does not move almost at all in this mode, nor is there almost 
any twisting in the fuselage or wing. Equation (132) shows that the 
wing bends symmetrically about the fuselage; equation (131) gives the 
mode shape of the fuselage. Such a mode is in fact a plunging mode, 
totally independent on the angle a and could possibly look as illus- 
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trated below: 

I 2 

t 



Of course, the analysis has only the value of an example and no claim 
is made as to the possible existence of such a mode. 

- Example 2: 

For such a hypothetical mode we assume that the values of A^ and 
A^ are much larger than the values of the other Lagrange multipliers. 
Then equations (125) - (130) give the following zeroth order approxi- 
mation: 

On = 0 (134) 

? n " - Vl 7 sW 035J 

F 

n 


q n i:) = 

7 ^ 


(136) 
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p_(‘> = v 2 X 3 cosa V 0) 


(137) 


n 


(4) 


n 


7 C2) . 


l n 


■ p 2 y 2 


X <J> T (O)cosa 
o n 


(138) 


(5) 


n 


(2) A sinaO (0) 

= - Y 2 v 3 6 n 


(139) 


( 6 ) 


n 


It is easily seen that in such a mode the fuselage twists only; an even 

more refined analysis can be done by comparing X^ and A & and hence the 

generalized coordinates for given values of Y2> v 2* v 3’ u 2" 

It can be noticed that for a = ir_ such a mode cannot be possibly 

2 

encountered. 

Indeed, equations (134) - (139) reduce to: 


o -^ (2) =0 


ti n 


P = - y v 

n 


- w°> 

r (3) 
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T 2 V 3 


Vn<°> 
F t 6 ) 


(140) 


and equation (140) introduces coupling between rolling and pitching 
modes. In fact the values of and X^ are obtained by totaly indi- 
pendent analysis. 

For the general case, substitution of equations (125) - (130) 
into the constraint conditions given in equations (113) - (118) 
leads to the set of linear algebraic equations in terms of (X^) i = 1, 
...6 from which the frequency equation will result. 

This system of equations can be written in tensori/al form 

_ -V -* 

A . X = 0 (141) 

with self explanatory notations for X and with the second order 
symmetric tensor A defined by its matrix representation as follows: 


11 


!_, (p, *„ 2 <V * 


( 1 ) 


(3) 


n 


n 


A = Z u ^n^l^n^2^ 

n=l F ^ 

n 


A 13 = A 15 = A 23 = A 25 =A 36 = A 45 = 0 


(142) 


(143) 


(144) 
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4 14 = £ , >Vl 
n=l 




a) 


n 


(145) 


A = Z y y. 

It) -« J. i. "" " ” -11 /'ll 

n=1 F (!) 

n 


(146) 


A 22 = ” (u, 'n'V + y = !> b 

11 n-1 F (1) p (5) 

n n 


(147) 


* 24 ’* !W!W 

n=l F (1) 

n 


(148) 


4„ - ? p Y *n fC 2 ) ' ,, 'n t5 2 :i 

26 n-1 1 1 7TI) 

n 


(149) 


. ® , 6 2 (£..) + 0 2 (0) cos a + (J> f 2 (0)sina) 

A,, = 2 (Y^, n '•^1-' n y n (150) 

n=l p (2) 1 F (4) F (3) 

n n n 


2 . z 

K r 0 (0) <£ * (O).sina cosa 

A,, = z (-Y,v- n + v n 1 ') 


n=l 


(4) - (3) 


n 


n 


(151) 


4„ = f Y,v ] W W 

n= l r (2) 


n 


( 152 ) 
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, s , 2 2 (E.) + 9 2 (0)sinS * v 2 [0)coso 

A 44 = 1 , C “l Y l V, 1 V 2 _2L 1 " 


n=l 


Cl) 


) C153) 


(4) n C3) 


n 


n 


n 


A,. = " u,Y, 2 *'n tE l , *'n (5 2 ) 

46 n=l 1 1 TJT) 

n 


(154) 


55 


Z 

n=l 


c Vi 


n 


CO 

n 


Y 1 Y 2 V 3 


e n 2 (0)c° s a 

jTW 

n 


U 2 Y 2 


4>' n 2 (0)sina 


(155) 


(5) 


n 


56 


00 

Z 


n=l 


t-’iVj 


9 (0) + .. 2, 


n 


( 6 ) 


2 2 


n 


|i , ri (0)^sina cosa 

jTW 

n 


(156) 


66 


= Z 


n=l 


(u i T i 


2 <(.• 


n 


(5 2 ) 


CD 


Y 1 Y 2 V 3 


0 n 2 (O)sina 

_ ( 6 ) 


; (5) 


2 m 2' 


(157) 


n 


n 


It can be pointed out that the elements are in fact mathe- 
matically Green's functions or in physical language influence co- 
efficients and can be easily understood. For instance A represents 
the bending displacement of the structure at the connection point 
between the fuselage and wing due to a unit force applied at the point 
connecting fuselage and tail. 

It is interesting to note that a Green's function approach to 
eigenvalue problems *as suggested by OtakarOanek in 1971, how- 
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ever no correlations of the work with a variational approach to the 
problem is mentioned, nor a practical way of actually solving the 
eigenvalue problems given. 

The eigenvalue system (141) leads to the frequency equation 


det & (oj) = 0 


(158) 


necessary and sufficient condition for existence of solutions X. 

Although equation (158) can be reduced to a polynomial equation 
in cu , in order to obtain a reasonably accurate solution for the 
first few modes, the degree of the polynome has to be high and hence 
it is impractical and time consuming to obtain the roots by iterative 
methods. For this reason.it is preferable to solve equation (158) by 
a trial and error method, which besides its simplicity has the ad- 
vantage of giving information about the slope of the determinant in 
the vicinity of the natural frequencies. 


c) Remarks on the sensitivity of the natural frequencies with 

respect to small changes in a from the symmetric configuration . 
The ultimate purpose of the modal analysis for the oblique air- 
plane is to determine the functional dependence of the natural fre- 
quencies upon a and also the changes in the modal shapes due to 
variations in the oblicity angle. It is by now obvious that these 
two questions cannot be solved entirely by analytic procedure, in 
other words it is impossible to express the natural frequencies 

ui (a) in the form 
n ^ 


w n = f n 
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where / (a) are known functions of a. 

However, particular aspects of the problem can be investigated 
analytically; for instance it is of interest to determine the order 
of magnitude of the change in a given natural frequency when the 
structure is given a small oblicity angle from its symmetric con- 
figuration . 

Based on physical evidence it is consistent to assume that 

to^Ca) are continuous functions of a ; then, for small values of 

1 ct -7T | it is possible to approximate to (a) as follows: 

2 ' n 


a> n (oO = to n O/2) + (a - ir/2) du n 

da 


a= 7i/2 


(159) 


provided that du n 
da 

Using the relation 


exists and has a finite value. 

a = tt/2 


dw 


n 


da 


a = ir/2 


3 2 et(A) 
9a 

£et (A) 

9(o 


a= tt/ 2, (o=io. 


n 


(160) 


tt/2, w=co 


n 


and the expressions of A_ given in equations (142) - (157) it can be 


seen that ^n 
da 


exists and it is finite except for situations in 


which 


a = tt 
2 


5 D et (A) 
9(jJ 


a = tt/2, oj = U3 


is zero 


n 



Ill 


This type of situation occurs only if is a double root of 
the frequency equation for the symmetric configuration (for instance 
this would be the case if the structure were composed of two identical 
beams connected at their middle), generally the natural frequencies 
are only simple roots of the frequency equation. 

Substitution of equation (160) into equation (159) yields: 



(161) 


3. Tiet(A) can be computed as follows: 
9a 


6 

9 Det (A) = I Det (A ) 

3a 1=1 1 

where is the matrix defined as follows: 


(162) 


'V 15 - A ij 

if j ^ 1 for 

all i 

<Vn ■ !ii! 

for 

all i 


3a 


This results from the fact that the function determinant is 
linear with respect to row or column. 

By using equation (162) it is found that 

3 Det (A) is identically zero in w. 

3a a * tt_ 

2 
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Ps a consequence, if the structure is given only a small 
oblicity angle B from the symmetric configuration there is no 
change in the natural frequencies at the order B . (B =a -tt/ 

A more refined analysis can be done by looking at the second 
order problem: 


2 

w (a) = u 0/2) + 1 (a - tt/ 2) 2 ** n 

n n 2 ,2 


da 


(163) 


a - tv/2 


Because 9 D et (A) 
9a 


= 0 for all u) 


n 


a = it/2 


,2 

d u). 


n 


da 


a = 7r/2 


9^et (A) 
9a 2 

9 0 et Q) 
3co 


a = tt/2,w=(o 


n 


a = tt/ 2, oo=co 


n 


and hence d 


n 


da 


exists under the same conditions as 


a = tt 
2 


dw 


n 

da 


a= ti_ 
2 


By substitution into equation (163) it follows that 


9 ^iet(A) 


w n (a) - to n (ir/2) = 1_ (a- tt/2) 


2 9a 


a= 7r/2,a)=aj 


9 D et(&) 


9to 


a = 7r/2,cu=a). 


n < 


(164) 
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3 2 Det (A) is computed by using again the 

a 2 
3a 


linearity property: 


2 6 6 

3 Det (A) = £ l Det (A ) 

~ 2 ~ k=l 1=1 

3a 

where 

\l= <Vl 

with the previous notations. 

As A^ and A^^ are independent of a for all i. The summation 
above extends in fact only to the subscripts 3, 4, 5, 6. 

A practical application of the analysis is given in the example 
below. 

This example assumes a pivoting wing version of the Lockheed 
L-100 (for which the structural frequencies and modes have been 
investigated in Chapter III) and by use of equation (132) predicts 
its first four natural frequencies for an 80° configuration. The 
results are tabulated below: 


3 Det (A) 
3a 2 

I 

|a = r/2,(i)=aj n 

3 S 

to 

o 

o 

CO 

o 

o 

3 2 Det(A) ! 
3a) 

a = 7r/2,co=a) n 

0.00016 


4.20 

4.200024 

0.00206 


8.30 

8.29969 

0.0058 


12.03 

12.03087 

0 . 8666 


16.78 

16.65 
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The results show that for a 10° deviation from the symmetric 
configuration only the fourth natural frequency changes by a 
significant amount; moreover it is seen that the first frequency 
remains virtually unchanged while the second frequency is slightly 
decreased and the third frequency increased. 

4. Numerical analysis: a pivoting - wing version of the Lockheed 

model L-100 Hercules 

The present section provides a more extensive study of the 
previous example; the modal behavior of such an airplane is in- 
vestigated for ranges of pivoting angles going from 90° (case of 
the straight wing] up to 30° (reasonable limiting value for practical 
purposes) . The value of the first four frequencies as functions of 
the angle a are plotted in Figure 38; for the straight wing case the 
type of modes is indicated (rolling or pitching) , 

As predicted by the perturbation theory, it is seen that no 
significant changes occur up to the fourth natural frequency where 
a decreases slightly from the 90° value; moreover, the trend of the 
variation is correctly predicted by the perturbation theory. 

It is seen that the first natural frequency is insensitive to 
values of a; the second natural frequency decreases significantly 
only when a varies from 60° to 30°; the third frequency can be con- 
sidered to increase as early as a = 70°, and the fourth frequency is 
seen to decrease even more rapidly. 

The fact that the first natural frequency remains constant 
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through changes in a is not in fact surprising; in the analysis developed 
in Chapter III it has been pointed out that the corresponding mode is 
practically a "wing mode"; all other components show rigid behavioT. 

The values of the ratios X^/X^ are tabulated below and give more 
detailed information about the first mode shape. 


a 

W 

Yh 

T Al 

x /X 
s' 1 

Yh 

90° 

1— ■ ■ 

-0.02 

- 


- 

- 

cr\ 

o 

o 

-0.02 

0.3 

0.2 

-0.26 

0.26 

45° 

-0.02 

-0.03 

-0.03 

0.02 

0.03 

04 

o 

o 

1 

o 

o 

to 

-0.017 

-0.03 

-0.014 

0.03 


The results show that for any angle of rotation the non-dimensional 
Lagrange multipliers X^to T. are small compared to X^; as a consequence 
both the fuselage and the wing will keep a symmetric behavior through 
changes in o f, practically the behavior of the straight configuration. 
These conclusions follow directly from equations (125) and (127) . More 
specific results can be drawn from the table above for each configuration 
in particular, it can be seen that when a = 60° the value of X ? /X^ is neg 
ligible compared to any X^/T^i ^ 2; in this case equation (129) shows 
that the tail will undergo mainly antisymmetric vibrations and since 
is less than the first corresponding free free natural frequency of the 
tail, its modal behavior will be a rigid rotation. 

Same argument applies to justify a rigid behavior of the fuselage. 
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As a decreases further, the values of the ratios A. A, i ^ 1 
become comparable and hence the antisymmetric characteristics alternate 
in the tail’s mode. 

A more complete picture of the first mode is given in Figure 39. 
The behavior of the airplane in the second mode can be analysed in 
similar fashion, given the values of the ratios i = 2, ....6. 


a 

V A i 

Vh 

Vh 

Vh 

Vh 

0 

0 






60° 

-3.56 

. 54.60 

36.16 

-6.32 

-3.04 

45° 

-4.0 

33 . 33 

! 

39 . S3 

-2.71 

-2.11 

30° 

: -4.94 

1 

22.36 

48.92 

-1.12 

-1.30 


A first important remark is that the values of T^/T^i = 3,5,6, 
decrease as the wing and tail pivot from the straight configuration; 
this shows that the second mode which originates from a mode of the 
rolling type becomes progressively coupled with pitching. 

It is also interesting to notice that ~ 4 /X 1 , T 4 /X 2 and X^Tg. are 

large compared to unity for the different values of angle o tabulated; 
as a result, the fuselage will mainly undergo antisymmetric vibrations, 
in fact rigid rotations in this range of frequencies. 

Because A^A^ and A 4 A^ are large compared to unity, another general 
feature of this mode is the antisymmetric behavior of the wing. 

Finally, it can be seen that the plunging mode becomes increasingly 
important in the tail as it rotates from the straight position to 30°. 
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The evolution of this mode with the pivoting angle is shown in 
Figure 40; the displacements have been normalized so that the larger 
generalized coordinate for the wing be equal to 1. 

Figure 41 shows the physical displacement of the components in 
the third mode with same normalization as in Figure 40; this mode is 
similar in nature to the second natural mode. 

The analysis can be pursued for infinitely many modes, but very 
seldom more than the first 10 or so modes are needed in a flutter in- 
vestigation; they are very easily provided by the present analysis. 
As an example, the first 10 natural frequencies are given for the 
straight wing configuration (Chapter III) and for a 45° oblique wing 
and tail: 

Straight wing oblique wing percentage difference: 


0) 

n 

mode type 

w n 

I“n90° - “n45° 




“n90° 

4.20 

pitching 

4.20 

0% 

8.30 

rolling 

7.87 

5.2% 

12.03 

rolling 

12.5 

3.9% 

16.78 

rolling 

15.2 

9.0% 

18.75 

pitching 

20.51 

8.9% 

21.65 

pitching 

21.83 

8 . 3% 

23.5 

pitching 

23.82 

14.0% 

24.3 

rolling 

37.54 

52% 

25.14 

pitching 

40.41 

61% 

40.09 

rolling 

48.50 

21% 

40.6 

pitching 

53.88 

26% 

i 
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The results show that the higher frequencies are more subject 
to change when the wing and tail rotate; moreover the coupling between 
a particular bending mode and a torsional mode has the effect of making 
the higher of the two frequencies higher and this trend is shown in the 
numerical results. 

5. Discussion 

The validity of the numerical results for such types of problems 
is generally subject to the assumptions made in their modeling (in the 
mathematical sense) and also very strongly dependent upon the methods 
used in handling them. 

If the basic assumptions for the problem are suggested by geometric 
or elastic considerations, it is relatively easy to make a choice giving 
satisfactory results; a more detailed analysis has to be made concerning 
the method by which these results are obtained. 

For instance, numerical results obtained by a perturbation method 
can be subject to doubt if it is noticed a posteriori that the small 
parameter was in fact of the order of unity. The present discussion 
has the purpose of precisely pointing out the situations in which the 
results given by the previous analysis have to be considered very cri- 
tically. It has been mentioned that the important parameter in the 
problem is the pivoting angle a defined between 0 and 180 degrees. The 
case where a = 90° is a limiting case of uncoupled motion, case in which 
both the elements of the eigenvalue matrix A . and A vanish because 
of the factored term sina cosa. If this result is similarly true for 
a = 0 or 180 degrees , care should be taken to avoid drawing conclusions 
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from it. Indeed, the latest situations correspond to a structure 
made of superimposed beams for which point contraint conditions are 
unsuitable. The constraint conditions should either apply along the 
line of contact in the one-dimensional theory, or else in this type 
of vibrations the beams would perpetually hit each other. 

The last remark would be true also when a is close to 0 or 180 
degrees; such a type of constraint condition (assuming that its formu- 
lation made sense) would be non-holonomic and thus the use of the 
Lagrange multipliers would be invalid. 

It is with this consideration in mind that the previous analysis 
has not been pursued for oblicity angles of less than 30°; these types 
of situations are anyway very difficult to realize in practice. 

The aforementioned limitation appears thus to be of theoretical 
rather than practical importance; but it is necessary to have it in 
mind when the method is used^for a critical analysis of the results. 
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CHAPTER V 

THE RAYLEIGH RITZ COMPONENT MODES METHOD USING 
LAGRANGE MULTIPLIERS APPLIED TO PLATE PROBLEMS. 

1. Statement of the problem and basic assumptions 

The problem of the lateral bending of plates is considerably more 
complex than the one-dimensional problem of the beam theory. While the 
eigenvalue problems leading to the natural frequencies and modes of 
vibrations of uniform beams with classical boundary conditions can 
be solved analytically, this situation is completely different in 
the two-dimensional theory of thin plates. Analytical solutions are 
known rigourously only in very few cases (eg. the case of circular 
plates by use of Bessel functions, or the case of rectangular plates 
simply supported along at least two opposite edges) . 

For the other types of boundary conditions numerical methods 
must be considered; for design purposes they should be general enough 
so that they still apply when the plate is constrained in a slightly 
different way. Unfortunately, most of the methods used so far have 
been ad-hoc methods . 

The purpose of the present Chapter is to show that the Rayleigh- 
Ritz Component Modes Method, through its use of the Lagrange Multi- 
pliers, applies to a very large class of problems and gives accurate 
results for less effort than the corresponding ad-hoc methods. The 
method requiers mainly basic intuition, (usually suggested by the 
physics of the problem)and is very flexible in its use. 




121 


The bending vibrations which are under consideration correspond 
to the classical small-deflexion theory of plates based on the follow- 
ing assumptions: 

(1) Points which lie on a normal to the mid-plane of the undeflected 
plate lie on a normal to the mid-plane of the deflected plate. 

This assumption corresponds to the dual assumption in beam theory 
that ’’plane sections remain plane” and ’’deflections due to shear 
may be neglected”. 

(2) The stresses normal to the mid-plane of the plate, arising from 
the applied loading, are negligible in comparison with the stresses 
in the plane of the plate. 

(3) The slope of the deflected plate in any direction is small so that 
its square may be neglected in comparison with unity, and the de- 
flections are small in comparison with the thickness of the plate. 

C4) The mid plane of the plate is a neutral plane, i.e. any mid plane 
stresses arising from its deflection may be ignored. 

Based upon these assumptions, the following probl ernes are considered 

2. Lateral vibrations of cantilivered plates 

As mentioned in Chapter I, this problem has references in the 
literature and thus the results obtained by the present method can be 
checked for accuracy with numerical results already existing, 
a) Principle of the Method 

The structure under consideration is shown in Figure 42; it con- 
sists of a rectangular plate of length L and width 21, clamped along 
the edge x = 0. 
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Its normal modes are sought in a Rayleigh Ritz expansion in terms 
of the modes of a completely free plate,with the clamped condition 
at x = 0 reinforced by means of Lagrange multipliers. 

Although the normal modes on an unconstrained plate are themselves 
given by approximate numerical methods, the foregoing development has 
been adopted in view of its generality as being valid for any type of 
support conditions . 

These modes themselves are found in a Rayleigh Ritz analysis in 
the form: 


Oo oo 

(x,y!> = E E 

k=l m=l 


( n ) 

\m 


(165) 


where X. and Y. are the natural modes of a free-free beam and the 

l l 

coefficients A^ 11 ^ are related to the natural frequencies in the 
following eigenvalue problem: 


K 2 A.t n > = A.< n) Ck 4 * s 4 k. 4 ) 
n 13 lj 1 j J 


2 °° 

+ S Z 


(n) 


1=1 r=l 


lr 


(v (to., uj . + to.. to. ) + 2fl-v)a.,ot. ) 
v v ll rj li K J ll 


where 

k^ : natural frequency parameter of a free-free beam 
s : aspect ratio of the plate 
v : Poisson ratio 


to. . 




jx!'(q x^o d£ 

non dimensional length of beam 


(166) 
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“ij -J X 'i'« X V« d5 

non dimensional length of beam 

The normal modes of the cantilivered plate are thus sought in the 

form: 


Wtx.y,') ■ I £ ' V« \m n) 

n=l k=l m=l 

(167) 

b) Problem of the constraint conditions 


According to equation 

(167), W(x,y,t) must 

satisfy the following 

boundary conditions: 



W(x = 0, y,t) = 0 

■ 1 1 y 1 1 

(168) 

9W(x - 0,y,t) = 0 

>\ » 

- 1 <_ y <_ 1 

(169) 


3x 

The two boundary conditions above are functional equations. 

They can be reduced to one single boundary condition by con- 
ceptually considering an unconstrained plate of double aspect ratio 
in symmetric or antisymmetric vibrations about the y axis, as shown 
in the lower part of Figure 42. 

In the first case the only condition which has to be enforced is 

W(0, y,t) * 0 for - 1 <_ y <_ 1 (170) 

in the second case the boundary condition to satisfy is 

9W (0, y,t) = 0 for - 1 <_ y <_ 1 (171) 

3x 

Because the analysis of Chapter III shows that a more rapid convergence 


>■ 


<7 
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is obtained when using symmetric free-free beam modes, the first 
possibility has been chosen. Then the summation in k in equation 
(167] extends only over the even beam modes. 

The additional term in the modified variational integral is: 


1 

J X (y) W (0,y,t) dy - 0 (172) 

-1 

Direct manipulation of this integral is very cumbersome in the 
Lagrangian method and therefore the integral has to be approximated. 
Three basic alternatives are considered in the present analysis; the 
theoretical development is pursued for each one of them and the results 
compared 


* First alternative 

The following approximation is taken for the left hand side of 
equation (172) : 

1 R 

1 X (y) W(0,y,t) dy = lim I X(y)W(0 ,y,t) «(y-y r ) (173) 

-1 R -+■ » r— 1 


where 



< y < 

V 2 


< 


R-l 


1 


* Second alternative 

It is a step function approximation of the integral : 
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1 R-l 

; A(y)W(0,y,t)dy = lim l A (y ) W (0 , y , t) Cy r+ 1 - y ) 074) 

-1 R -*• 00 r=l 

where the definition of the y^ is the same as in the Dirac approxi- 
mation. 

* Third alternative 

This alternative consists in computing the integral by a trapezoidal 
rule : 

1 , R-l 

5 A(y)W(0,y,t) dy = lim 2 (A(y )W(0,y ,t)+A(y )W(0,y . ,t))(y -y 

-1 R-oo r= i_^ L Ell E1 1 ... L 

with same definition of the y . 

r 

In each one of the cases the following notation is adopted: 

A (y ) - A . 
v/ r r 

c) Theoretical development; frequency equation 

The kinetic energy of the unconstrained plate which mass density 
per unit area is p , is given by: 

1 1 

T = 1 p L 1 1 1 W 2 (s,n,t) d£ dn (176) 

2 -1 -1 


Since 


W(C,n,t) 


= E 

n=l 


QnCt^Ce.n) 


(175) 
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and by virtue of the orthonormality of the plate inodes it follows: 

(177) 

If D is the bending stiffness of this plate of thickness h 

(D = Eh ) , the potential energy can be written as : 

12(l-v 2 ) 
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* in the first alternative 

The modal expression of the constraint equation is: 


lim Z Z X 0 (t)* (0,n ) = 0 
R -*■ » r=l n=l 


so that the Lagrangian of the elastic structure is: 


of = Z ( _1 pLl 0 2 (t) + 1 D1 K 2 0 2 (t) + 

n=l 2 2 .3 n 

L 


E A r 


The subsequent Lagrange's equations are thus: 


pLlQ^t) + D1 K \Ct) - Z X r*n (0 *V = 0 n s ^ 
L 3 r=1 


^ Q_Ct) * (0,n ) = 0 
n=l 71 


r = 1 


With the assumed harmonic time dependence of the following form: 


- On 


i - 2 (i ) 1/2 

e L 1 


x r (t) 


X D1 e 

r 7 


1 72 B) 1/2 - 
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The generalized coordinates are solved in terms of the Lagrange 
multipliers from the first set of equations: 


On " - 



2 

U) 


V°’V 

- 


(.182) 


and their substitution into the second set of equations gives the 
eigenvalue equation: 


(to) 


A = 0 


where 


( 1 ) 


A 

TV 


<30 


I 


n=l 


$ 


n 


(0,n r ) * n (0,n v ) 



(183) 


r = 1, . . .» ; v = 1, . . .» 

In the numerical analysis both the eigenvalue matrix and its series 
elements will be given a finite form, i.e: 


n = 1,2— N 
r,v = 1,2 R 


N finite 
R finite 


The frequency equation is 
det (A^((o) ) = 0 


( 184 ) 
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and its solutions are the natural frequencies of the cantilivered plate 
initially considered 

* in the second alternative 


The constraint equation can be written in modal form as: 

R -i M 

lira E z 1\ o (t)* (0,nJCn^, -nj (185) 

R-x” r=l n=l r n n r r+i r 

and the new expressions of the Lagrange’s equations are: 

. pLl^Ct) + D1 K n 2 Q n Ct) - Z i iVn^’V^r+rV = 0 Cl86) 


n = 1, °° 


V 0 ’ n r ] C Vl ■ V ■ 0 


( 187 ) 


r = 1 , 


Same harmonic time dependency is assumed in this case, and 


. 1 


.2 »■ - 


( D ) 1/2 u t 


x r Ct) = X r D_ e 
L 3 


The generalized coordinates are given as functions of the Lagrange 
multipliers by: 


On = 


2 *r V°’V<Vl - V 


r=l 


2 Y 

£jl) - K 


n 


2 


( 188 ) 
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and the elements of the eigenvalue matrix are in the finite form: 


( 2 ) 

A 

rv 


N 

E 

n=l 


V 0 - n r i V 0 - n v K Vr' 1 r ) ( vr\i 


2 

OJ 


2 

K 

n 


( 189 ) 


r,v = 1, — R-l 


It can be further noticed that A 


( 2 ) 

rv 


can be written in fact as : 


A 


( 2 ) = 

rv 


(n 


r+1 


V 


Cn 


v+1 


- n ) A, 

v 


CD 

rv 


and hence 


( 190 ) 


This last equation shows that the difference between the two alternatives 
is that for the same finite number of plate modes and same accuracy 
expected in the frequency, it is necessary to take one less point 
in the subdivision of the y axis. The practical meaning is that the 
order of the matrixes handled in the numerical computations decreases 
by one in the second case. Although this fact is an improvement, it 
is quite negligible compared to the improvement brought by the procedure 
described in part a) which reduces by half the order of the eigenvalue 


matrix. 
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* in the third alternative 


The constraint equation can be written in the double serie form 


as : 


R-l 


liL E , E , V t,1( VrV t> r.lV 0 'Vl ) * > rV 0 'V ,) = 0 (191) 

r= l v ' 


and the subsequent Lagrange's equations are: 


0 L1 (^(t) * m Vt) - 1 J C n r4l - n ) 

r * E : — 


(A . $ CO , n ) - X <t> (0,n )) 
v r+1 n v r+l' r n v r J 1 


(192) 


n = 1 , . . . °° 


• with R “ 


l =1 V t} V 0 ^ Cn 2 “ V = 0 


(193) 


00 

J =1 °-n Ct) V 0 'V ( vrvp = 0 


r = 2,... R-l (194) 


l =1 V°’V (n R ■ VP • 0 


( 195 ) 


With same non-dimensionalizat ion as previously. 
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R-l 


- E * Vn(0, n rl) 

= • r ’" 


2 V 

oj - K 


n 


C196) 


and the eigenvalue matrix 




of order R is defined by its elements: 
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A rR^ = 1 ^n^’V *n^ 0 * n R^ n R“ n R-l^ n r.+ l ’ n r-P 

11 = 1 — ~ 


2 „ 
co - K 


n 


r = 2 , • - - R-l 


A RR } ? (°’V fn R n R-l } 

n=l r 

co 2 - K 
n 

The eigenvalue determinant can be written after factorisation: 

R-l R- 1 

det (A^ (w) ) = rx C^ r+1 - ^ r ) 2 f n r+ i _n r _i^ 2 detCA^hu)) (197) 
r=l r-2 

and hence this approach is seen to be completely equivalent to the 
first approach considered. This fact is somewhat surprising because 
generally the trapezoidal rule provides better approximations of 
integrals than does the Dirac functions method. 

It can be noticed that in the process, the present method gives 
also the natural frequencies of point-clamped plates; the number 
of points and their location along the edge is completely arbitrary. 
Furthermore, for each natural frequency the distribution of the 
shear forces and/or bending moment is given along the constrained 
edge at the points of constraint. 



d) Natural frequencies of a square cantilivered plate 

The first three natural frequencies of a square cantilivered 
plate are computed by the present method using the normal modes and 
frequencies of an unconstrained plate of aspect ratio 2. The 
analysis of Chapter III shows that five symmetric modes in the x 
direction insure good convergence (they correspond to less than 
0.1% error in the beam case); the number of free-free beam modes in 
the y direction is chosen as a function of the number of points at 
which the constraints are applied along the y axis and as a function 
of the number of natural frequencies sought. For the present analysis 
five free-free modes have been taken in the y direction. The changes 
in the first three natural frequencies with the number of constraint 
points are shown in Figure 43, both for an approximation by Dirac functions 
and by characteristic functions. As mentioned earlier, the second 
alternative gives faster convergence and the results will be dis- 
cussed for this approach. The asymptotic values corresponding to the 
square cantilivered plate have been taken from Reference (35) . 

The lower part of Figure 43 shows that as the number of constraint 
points increases, the first natural frequency increases monotonically 
towards that of the cantilivered nlate. The results are further 
analyzed in Figure 44 where the percentage difference with the canti- 
livered case is plotted versus the number of constraint points. 

It can be seen that with as little as three constraint points 
(one in the middle and two at each extremity) the plate can be con- 
sidered clamped at the edge x = 0 in the first mode, with an error of 
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only 1%. This result is to be expected because cant iii vere d 
= 3.494 is very close to the first natural frequency of a cantilivered 
beam 3.516 and because in this mode the plate undergoes only 

translation in the y direction in a "beam-fibers” type of behavior. 

The second natural frequency is seen to converge as well towards 
that of a cantilivered plate when the number of constraint points 
increase, as shows the middle portion of Figure 43. In this case the 
variation is not monotonic per se, and this fact will be commented 
upon later. 

The percentage error from the case of the cantilivered plate 
is plotted in Figure 45; it can be seen that with the characteristic 
function approach^a square plate clamped at the points = -1, 

= -0.5 , rij = 0 can be considered as completely clamped along the 
edge in the second mode, with an error of- only 0.7% in the natural 
frequency. It is interesting to notice that as far as convergence 
is concerned, there is no improvement in constraining the plate at 
n l = _1, n 2 = ° Up0n having only the point Pj = -1 clamped. This 
result comes from the fact that the modal behavior of a cantilivered 
plate in the second mode in the y direction is a combination of 
plunging and rotation. By imposing zero displacement at Pj = -1 
we constrain Y q and Y ^ , as for the point 5 0 it is a redeen- 
dant point for this particular type of constraint because Y^fO) = 0. 

The method is seen to converge remarkably fast also in the 
third mode as shows the upper part of Figure 43; this result is 
even more remarkable considering that only 5 modes are used along 
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the y axis for as many as 3 points of constraint for instance. For 
this mode also, the percentage error with the member of constraint 
points is plotted in Figure 46 and the corresponding mode shape of a 
squared cantilivered plate is given for reference. 

Higher natural frequencies could be obtained provided that more 
modes are considered in the y- direction in the analysis giving the 
normal modes and frequencies for the unconstrained plate. 

For cantilivered plates with higher aspect ratio, same conclusions 
are expected to be valid with an even smaller number of constraint 
points. Also, the method is expected to converge even faster for 
edge conditions physically closer to the free condition (as simply 
supported, for instance). 


3. Vibrations of a rectangular plate with an internal support 

As mentioned in Chapter I, an analytic solution is proposed to 

; 

the problem of symmetric-symmetric vibrations in Reference (36) 
for the following configuration: 

y ‘ 

i 



The mixed boundary value problem is formulated in terms of dual 
series equations; they are first partially solved by means of the 
Bessel function and then cast into the form of Abel's integral eouation 
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for which the solution is given explicitely. 

Although the method is mathematically interesting, it only 
applies to the particular foregoing problem and would fail to give 
results if the line of support were to be displaced in the plate, 
a] Principle of the method. 

As in the previous paragraph, the mode shapes of the structure 
are sought in terms of mode shapes of a less constrained plate. The 
physics of the problem suggest these initial modes to be the modes of 
a simply supported plate along all four edges; mathematically, they 
offer the advantage of being known by their analytic expression. 

Hence the modal solutions are sought in the following serie form: 


W(x,y,t) = 22 E 0 (t) sin mirx sin n-rry 

m=l m=l L 1 


(198) 


0 < x < L ; O^y^l 

For motions symmetric about x = L/2, y = 1/2, the summations in m 
and n extend to the odd subscripts only. 

The boundary conditions at the edges of the plate are already 
satisfied by the initial choice of the spectral expansion; the only 
additional constraint to satisfy is that along the internal line of 
support , i . e . : 


- c < 

1 


x < c. 


W(x, i , t) = 0 
2 


for 
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when only doubly symmetric modes are sought, this condition further 
simplifies to: 


W(x, _1 , t) * 0 for _L - c. < x <_ _L 
2 2 1 2 

The above equation requires the function f(x,t) = WCx^.t) restricted 

2 

to L_ - c < x <_ _L to be identically zero in this interval. 

2 2 

Assuming the function to be continuous in x in the interval of 

definition, if we choose a dense set (x ) „ L - c„ = x,< x„< 

r r=l,R 2 * 1 12 

< x R = L it is sufficient to require f(x R ,t) =0 r = 1,....R 
In the limit, equation (199) is equivalent to the set of equations: 


W(x 1_, t) = 0 r * 1 , R (200) 

2 

when R -*■ <». 

The eigenvalue frequency equation will be derived based on this 
form of the constraint conditon. 

b) Lagrangian formulation of the problem; frequency equation. 

The total kinetic energy of the structure is : 


11 . 

T = i pL1 0 r o W (5>n,t) dn 


( 201 ) 


where p is the mass density per unit area of the plate. 
By use of equation (198) it follows: 
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( 202 ) 


Likewise, D being the bending stiffness of the plate and s its aspect 
ratio, 
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or, after substitution of equation (198) 


U=J_in! l Q 2 (t) ( (m7r) 2 +s 2 (n7r) 2 ) 2 
2 l 3 m=l, 3 n=l , 3 ^ 


In equivalent form, 


u -ijn e 

2 ^3 m=l,3 


E 

n=l , 3 


<(*5 K 


2 

mn 


C203) 


where 

K = (mir) 2 + s 2 (nn) 2 
mn v J 


By substitution of (198) into equations (200), the set of constraint 
conditions is obtained in modal form as follows: 
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E 2 (-i)^' 1 ^ 2 o (t) sin mTtC = 0 r=l,R 

m=l ,3 n=l,3 ™ n T 


where E, = — r- 
r 


( 204 ) 


l 


If X^ is taken as notation for X(£p, the Lagrangian of the system 
can be written as 


JE, = 


E 

m=l , 3 


E 

n=I ,3 


1 _ pLl 
2 



1 Dl_ K 

2 L 3 


urn 



(t) + 


C-D^-O/Z Q^ct) I J r sin n,»5 r ) 

r=l 


and the subsequent equations of motion are: 


pL1 + ^i c 

r. 


R 

- E 


r=l 


, ,.(n-l)/2. 

(-1) X Sinun7TC 

r r 


n = 1, 


(205) 


t E 


00 

E 


m=lj3 n=l,3 


(-l)*- n 1 ^ 2 Q jrin (t) sin imr? r = 0 r=l,..R 


( 206 ) 


In a harmonic time dependent motion of the form: 


^mn^ ^n 


i ~ ( D ) 1/2 w t 
k P 
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Vt) = 


D1 

.3 


~2 C 

L 


M2 


ait 


the respective amplitudes of the generalized coordinates are given in 
function of the Lagrange Multipliers by : 


'mn 


R (n-l)/2 

S (-1) 

r=l 


A sin nurE 
r r 


2 v 
u) - K 


mn 


and the eigenvalue equations are: 


(207) 



(208) 


where 

^,(“0 = E sin mirg sim m-rr£„ ? 1 

m=l , 3 r V n=l , 3 2 

*) z _ k 

mn 


(209) 


The natural frequencies of the structure are solutions of the equation 


det ( A (to) ) = 0 


c) Numerical results. 

The analysis of Reference (36) treats the case of a simply 
supported square plate and gives the variation of the fundamental 
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natural frequency with the length of the internal support. For 
purpose of comparison, the method is applied in the same case. Six 
symmetric modes are considered in the y direction and five more 
symmetric modes than constraint points in the x direction. 

The results are summarized in the table which follows: 


2 c i 

L 

* 

W Ref (36) 

* 

ClJ Present Method 

Percentage 

disagreement 

# constr. 
points used 

0 

2.665 

2.668 

0.9% 


Z = 0,5 

0.1 

2.898 

2.966 

2.6% 

Z 

: (0.45,0.5) 

0.2 

3.164 

3.160 

0.13% 

Z 

: (0.4, 0.5) 

0.3 

3.353 

3.360 

0.21% 

Z 

: (0.35,0.5) 

0.4 

3.451 

3.456 

0.15% 

Z 

: (0.3, 0.4, 0.5) 

0.5 

3.492 

3.480 

0 . 36% 

Z 

: (0.25,0.325,0.5) 

0.6 

3.507 

3.509 

0.06% 

Z 

:(0. 2, 0.3, 0.4, 0.5) 

1.0 

3.512 

3.514 

0.06% 

z 

: (0.1, 0.2,0. 3,0.4, 


It can be 

seen that except for 

2C 1 

the value — j-~- 

- 

there is less 


than 1% disagreement between the fundamental frequencies obtained by 
the two methods. Moreover, as shows the fifth column, the number of 
points approximating the line of support is never greater than 5 and 
hence the corresponding sizes of matrices are being kept small. 

The functional dependence of the fundamental frequency upon the 
length of the line of support is plotted in Figure 47 and indeed the 
results are in good agreement considering the scale used. 
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d) Conclusion . 

The numerical results tabulated in column four make any comment 
on the convergence of the method unnecessary. 

It has to be pointed out that such a good convergence is due par- 
tially to the particular choice of the initial modes which verify the 
boundary conditions at the edge and which are known in analytic form 
and hence do not have an intrinsec error. 

The problem can be solved alternatively by using unconstrained 
plate modes, but this choise has for effect in increasing the number 
of constraint conditions required. More precisely, the analysis of 
Reference (44) shows that for the first six natural frequencies a simply 
supported square plate can be approximated to a good degree of accu- 
racy by a point supported plate with four points of support along each 
edge, and thus the last alternative would include 12 additional pone- 
tual constraint equations. 

If the line of support were to be moved in the plate, there the 
method would still apply with the following modifications. 

(1) the sumation in m and n in equation (198) should be extended 
to both even and odd subscripts 

(2) the points of constraint should lie along the entire line of 
support and no longer be restrained only to its half. 

The first modification does not bring any significant increase 
in the computation time because it corresponds to operations of adding; 
the second modification would increase the size of the eigenvalue 
matrices by a factor of two (as their initial size was already small, 
such a modification does not represent any real inconvenience) . 
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CHAPTER VI 

CONCLUSIONS AND RECOMMENDED FUTURE WORK . 

1 . Conclusions 

The Rayleigh-Ritz Component Modes method using Lagrange Multipliers 
has been shown to predict with good accuracy the modal behavior of 
linear elastic structures of various degrees of complexity; it has been 
seen that the analysis offers computational advantages through its 
convenient matrix formulation and clear physical interpretation. 

Good convergence has been obtained in several cases where other 
analytical results were available. Moreover, for structures with non- 
uniform properties, because the convergence is based upon several 
independent parameters, there is more than one procedure for obtaining 
good accuracy in the solutions. The discussion presented in Chapter II 
suggests a systematic procedure economical from the point of view of 
computation time. 

The method may be used either as a Rayleigh-Ritz type method or a 
Finite-Element type method. In the Rayleigh-Ritz option one employs a fixed 
number of components and systematically increases the number of modes per 
component; conversely in the Finite-Element option one fixes the number of 
modes per component and systematically increases the number of components 
(or elements) . Thus, potentically, it is more flexible in its use and has 
better convergence characteristics than either of the conventional methods. 

The main factors upon which the success in the application of 
the method depends are the following: 

1) an adequate mathematical model of the components of the structure 
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(within the framework of linear elasticity and small deflection 
theory) , 

2) an accurate knowledge of the natural modes and frequencies of 
vibration of the unconstrained components, 

3) a suitable formulation of the constraint conditions describing the 
connections between components and, of course, a suitable choice of 
the degrees of freedom for each component. 

As a consequence of 1), the present theory has the limitations in- 
herent in linear models. The second point mentioned above is particularly 
critical when one or several components of the overall structure are 
modeled by shells, for instance, because of the difficulty of ob- 
taining a good modal representation for them. This fact will be commented 
upon later. 

Finally, as is pointed out in Chapter IV, care has to be taken in 
checking that the constraint equations can be expressed in holonomic 
form so that the use of the Lagrange Multipliers is valid. The problem 
of two structures striking each other during vibration for instance, can- 
not be treated by the present method. For the problems for which the 
use of the Lagrange Multipliers is justified, the present method is 
convenient also for computing the internal forces developed at the 
points of constraint. The analysis of Appendix VII shows that the 
values of the Lagrange Multipliers are much closer to the actual values 
of the shear forces and bending moments at those points, than the 
values directly given by differentiation of the modal shape. 
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2 . Recommended future work 

The modal behavior of arbitrarily constrained plates with non- 
uniform properties can be investigated by extending the analysis 
developed in Chapters II and V; it is expected that the more sensitive 
parameter for the convergence of such problems is the number of substruc- 
tures (or uniform plate elements) considered, rather than the num- 
ber of unconstrained plate modes used for each component. (The latter 
depends upon the number of constraint points taken along each edge of 

the plate elements and their respective bending stiffnesses). 

Vibrations of structures with shell components can be also treated by 
the present method; however, the problem is then of greater complexity. 
First of all, it is very complicated to formulate and find solutions 
of the characteristic free vibration equations; while the form of the 
classical fourth order equations of motion for plates is universally 
approved, there are numerous shell theories which have been derived 
for a given shell configuration. Hence, first the theory adopted for 
obtaining the unconstrained frequencies and mode shapes should be chosen 
by checking the different assumptions which give rise to that theory. 
Additional complexity enters into the problem because of the question 
of the boundary conditions. The classical bending theory of plates 
requires only two conditions to be specified along an edge, while a 
corresponding shell theory requires four specified conditions. How- 
ever, there are already available results, primarily for circular cylin- 
drical shells, as can be seen in Reference (47). 

Another problem which is of interest in structural dynamics is that 
of the determination of the modal damping of a given complex structure. 
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The analysis presented in Appendix II takes into consideration only the 
structural damping per se. However, sometimes the damning occurring in 
structural joints is the major contributor to the total damping of the 
structure, and hence an open area for further research is that of mo- 
deling the joint damping so that it can be incorporated as a constraint 
condition in the present method. In Reference (48), some experimental 
attempts are presented for modeling the various types of joint damping 

occurring in the orbiting Skylab cluster. 

In Chapter IV, an analysis for determining the structural modes 
of an airplane with oblique configuration is developed; the results 
can be directly used for a flutter or gust response analysis. As 
a first step, a simple aerodynamic theory could be considered (such as 
piston theory, or aerodynamic strip theory for instance); then the 
problem can be pursued with more complex aerodynamic operators. 
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APPENDIX I 

MODAL BEHAVIOR OF NON UNIFORM CONTINUOUS BEAMS 
BY A PERTURBATION - GALE RK IN METHOD. 

1 . Introduction 

The purpose of this Appendix is to present an alternative method 
to the Rayleigh Ritz component modes method presented in Chapter II, 
for continuous beams with variable stiffness along the span. 

The method previously presented applies with great simplicity 
to the treatment of stepped or almost - stepped beams, where the number 
of components is determined by the number of ’’steps" of the beam. This 
number is then finite and the convergence problem identifies with the 
convergence of the Rayleigh-Ritz method. Because of the simplicity 
and efficiency of the analysis, hardly any other alternative method 
could be considered for such problems. 

For this reason, only the case of beams with continuously 
varying properties has been investigated by an alternative method. 

2. Theoretical development 

The variations in the bending stiffness of an elastic beam can 
be caused by the change in the elastic properties of the material 
(in particular the Young modulus) and/or variations in the geometry 
along the span. If changes in E occur, we shall assume them con- 
tinuous and smooth; as for changes in the geometry, only slowly vary- 
ing properties along the span are meaningful 1 within the one-dimen- 
sional theory of continuous beams. The previous considerations 
suggest a "perturbation - type" approach. 
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Unlike the method presented in Chapter II, this method does not 
adopt a variational point of view and the starting point is the differ- 
ential equation and the associated boundary conditions. 

a) Notations 

EI(x) : bending stiffness at station x 

El : mean value of EI(x) 

o 

p : mass density 

A(x) : cross sectional area at station x 

A : mean value of A(x) 

o 

y(x,t) : transverse displacement of the bent beam 

w(x) : normal mode shape of the bent beam 

6 : Kroenecker’s symbol 

np 

b) Development of the problem 

The differential equation governing the transversal vibration of 
a non-uniform beam is: 


Kj (El (x) 3 2 y (x,t)) + pA(x) 5 2 y (x,t) =0 0 <_ x <_ .L 

o 2 a+ 2 

3x 9t 


( 210 ) 


and can be written in non-dimensional form as follows: 


1 9 2 (Ei (O afx + pA(e) sfv a,t) = 0 


( 211 ) 


4 2 

L 2C 


n' 


9t 


The normal modes are found by a solution of the form: 


iJlt 


y(C,t) = w(C) e 
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which when substituted into equation (211) leads to: 

1 C (EI(5) 4 CO J - p*ft)AcO - o 

l dr dc 

We introduce the following definitions: 

f(0 = El CO - ei o 

(EI(0 - EI ) 
v v o'max 

and 

g(0 = ACO - A 0 

(ACC) - A ) 

K J o-'max 

The above defined functions satisfy 

I f CO 111 and I gCO !l I vce[-l, 1] 
and the quantity 

c = larger of ( I fU£ - , 1 ^ ~ A oL, ) 

EI A 

o o 

is small compared to unity. 

Using (213), (214), (215), EI(c) and A(0 can be written: 
EI(C) = EI q ( 1 + N(e) f (O ) 

A (?) = A q (1 + M(e) g(0 ) 


( 212 ) 


(213) 


(214) 


(215) 


(216) 


(217) 
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with obvious definitions of N(e) and M£e) from £215). 

By substituting £216) and £217) into (.212) this equation becomes: 


El 

-r- ( 1 + N(e)f(C) 
L 


N(e) ™ 

dg 4 L 4 dg 


d 3 W(?) 


dV 


El 

+ —J N Ce) 

L 4 


iiiD 


pA o £1+M£e)g£0) « 2 W(?) = 0 


£218) 


So far no specific boundary conditions have been assigned, in fact the 
method will be developed for any arbitrary type of boundary conditions, 
noted B.C. 

Solutions are sought in the form: 


W (O = ♦ (O + Z E Ce) * in (?) 
“ i=l 


£219) 


fl = a) + Z P. (e) 0. 
n n T . l m 
1=1 


C220) 


with 


i • E. 1 £e) . . F (e) 

lim l+l ^ J = lim l+l J = 0 


e-+0 


E.£e) 


e -+0 


F i £ e ) 


Substitution into £218) yields the perturbation equation: 


(1+N£e)f £0 ) £ d + E (e) d %l (g} + 

* a * a 


dt 


d, 4 
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V«J d *-2 [ » * ....) ♦'fL> NCC) ( 

dS 4 L 4 dc 


d V« „ , Al (5) r r , A,2<« . . 

' ,T3 _ + E l (e;) ~3 + E 2 (e) , 3 ••• } 

d? d£ dB, 


EI d 2 fCO dV CO d 2 S> .(5) 


dC dt 


dr 


E,(e) d $ n 2 C0 + . ..) - pA (l+M(e)gCO ) 
2 

dr 


(w 2 + 2F-(e)cu n , + F 2 Ce)Q 2 + 2F„(e)w fi _ +...) 
v n l v J n nl 1 1 ' nl 2 V n n2 


(♦ n CO + E iC £ )*niCO + E 2 Ce ^n2 (0 + •••) = 0 (221) 

From the boundary conditions B.C. we determine boundary conditions for 

the functions $ . (5) . 

m 

For example, in the case of a cantilivered beam clamped at the end 
5=0 the boundary conditions are: 

W(0) = dW (0) = 0 
d 

EI Cl) d\ (1) = d_ (EI (5) d^W (5) ) (1) = 0 
d ? 2 d? dC 
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from which we determine the boundary conditions on (?) using 

(219) 


* ni (0) = CO) = 0 
d? 


( 222 ) 


dV. (1) d 3 4> . (1) = 0 

m v J - m v 


(223) 


dr 


dr 


The problem is now welj defined at each order and solved as follows 
* Zeroth order solution 

It is as expected, <f> n (?) satisfying the zeroth order problem.: 


El 
o 

L 4 

B.C. 


d <t> n (€) 


- pA ui (?) = 0 
on n v 


(224) 


* First order solution (solution at order E (e)) 

For almost all the cases of interest the slow variations in El 
(?) and A (?) have comparable lowest order terms, which can be written 
using Landau's notation: 

N(e) =0 (M(e) ) 

(i.e. lim (N(e)/M(e) ) is finite) 
e -* 0 

Solutions are sought such that: 


Ej(e) = 0(N(e) ) : Fj CO = 0(M(e) ) 


(225) 
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and hence the first order problem is: 


El d 4 $ . (?) pA u> 2 $ ..CO 
o nl l ' -on nl v ' 

4 4 

L 4 d ? 4 


(f co J \g), 

de 


+ PA 0 ( 2 <V\iiV 5) + g( ° w n V C) } 


(226) 


We denote by the elements of the spectral decomposition of 


$ . , i.e. : 
ni 


*nl CO " S , ni 


(227) 


S CP) . J 
nd “ o 


* nl CO 4 (?) d? 


(228) 


and we apply a Galerkin procedure to equation (226) . 
It follows, using B.C. that: 


, 2 2 . (p) 

PA (oi - w ) s , = 

p n' ni 


T o r d (f(?) d \ ( 0 )a (?)d? 

T 0 75 ~75 p 


+ pA 0 C2 “iAii 6 np % o 8 CO+ n CO+ COd?) 


for n = p, ft is determined by: 


( 229 ) 


2 
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where is non-dimensional frequency parameter of the reference beam. 



The solutions at the next order are determined by a similar 
type of procedure. 


3. Beams with linear bending stiffness and mass distribution along the span. 
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A beam of such type is illustrated above. The width distribution along 
the span is given by b CO = b (1 + e?) where e = o(l) 

The analysis of the previous paragraph applies with 

fco = g(e) = e (23n 

life) = M(e) = e (232) 

and solutions are formed with 


Ef ( e) = (e) = e 1 


(233) 


for various types of boundary conditions. 

It can be noted that relations (231) and (232) hold not only for 
the particular type of beams shown in the figure, but are more general 
(there are multiple ways in which EI(x) and the mass distribution /unit 
length, m(x) can be made to vary linearly along the span). In fact 
formula (231) and (232) could have been taken as starting points. 


Free-free beams. 

The first order collections to the natural frequencies 
are found to be zero; the corresponding corrections to the mode shapes 
are : 


*nl<« 


E 

P =1 

p different 
parity than 


16a a k k A (?) 
n P B P P 


( 234 ) 


n 


where the values of the are given in Reference (21) . 
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On the basis of this result it can be concluded that the natural 
frequencies of such beams can be considered within a reasonable degree 
of accuracy as being the same as the natural frequencies of the 
uniform beams of reference, the degree of accuracy being determined 
in its major part by the value of (the next higher term in the 
serie (220},. It is also significant to notice from formula (234) 
that the corrections in the mode shapes at the first order are the 
corrections due to the fact that there is no longer symmetry about the 
middle of the beam, and thus although the boundary conditions are 
symmetrical, the modes are no longer even or odd functions. 

Finally, the ^ aTe determined by solving the second order 
problem and found to be: 



Clamped - free beams. 

The natural frequencies of such cantilivered beams are given at 
the first order by: 


£1 = ii) ( 1 - 2e an ) 

n n — sr 3 


( 236 ) 


and the associated modes: 
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w (O = <f> (O - £ 1 

" n P=1 

P^n 


8a a k 

n p n p 



n+p 


( C-i K 



>V° 


Equation (236) shows that at first order a> n bounds ^ from above if 
e>o and from below if £>o and that for a given £ the first order 
correction decreases in absolute value with n^ the order of the natural 
frequency sought . 


Clamped - clamped beams 

This type of boundary conditions leads to similar conclusions 
about the natural modes and frequencies as in the free -free case; such 
a result is expected because the free-free and the clamped- clamped 
modes in the constant case are related by: 


(♦_CO ) 

free-free 


2 

clamped -clamped 


Beams simply supported at both ends. 

The previous conclusions apply again and the natural frequencies 
up to the second order are found to be: 


ft = U) 

n n 


1+e [ 


00 
1- E 


16 


(2mr) ' 


p=l 

p/parity 
than n 


5 4, 2 2, 2 

n 7r (n -p ) 


, 2 

Op 

2 


i ) 


n 


The serie in the right hand side of the equation is convergent, its 

3 

general term being of the order of 1/p . 

The modes shapes are determined at first order by: 
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The other possible types of boundary conditions present no 
difficulty and the modal characteristics can be found from paragraph 
2 in similar fashion. 

It is to be noted that for the foregoing class of beams, the 
natural frequencies are closer to the natural frequencies of the 
beams of reference when symmetrical boundary conditions apply that for 
the other types of boundary conditions. 


4, Tapered truncated cantiliver beams with elliptical cross section, 
This class of beams is a reasonable approximation of high 


aspect ratio airplane wings and therefore it is of interest to 
investigate the modal behavior of such structures by the method 
developed in paragraph 2. The type of configuration considered is 
the following: ) Q 





16 n 


MO = h 0 - e £ ^ b(0 = b 0 " ae ^ 


o<e<<l ; 0 1 5 1 1 

The corresponding mathematical model is given by the differential 
equation (212) aid the boundary conditions B.C. are: 


y(o) = 


d^ 

dC 


(o) . 0 


(237) 


El Cl) Zjl CD 

d5 2 


d_(EI(0 Z_x. CO ) 


dS 


dO 


= 0 


(238) 


5 = 1 


The cross sectional area and structural moment of inertia are: 


A(0 = * MO MO 


(239) 


1(0 = n b C0 MO 

64 


(240) 


and consequently, the method of paragrpah 2 
functions : 

ae£ 


o"h 

A(0 - A = - irb ec ( 1 + o - 

o o — ~b 

o 0 


) 


applies for the particular 


EI(£) - El = En eil h 2 b ( 3 - ah o ) 
o — oo 

64 b 


t rL u f i ah >. 3 __ 3i « 3 ah 4 4 

3e£h b o ) - e £ b (1+ o^ae £ j 

0 0 0 f* * 
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The first order solutions for the natural frequencies are 


found to be 


Id ( 1+ T- 
n h 


2a f a . ah 
n ( 4 + o )) ) 


( 241 ) 


The first order correction in the right hand side increases linearly 

with ah o the taper ratio between the chord and the thickness; more- 
_ 

over thi§ correction is always negative for the natural frequencies 

for which > 2 Z - 4 + taper ratio. 

This formula is significant for design purposes because it indi- 
cates in which manner to adjust the taper ratio such as to increase 
the value of a resonant frequency which is within an undesirable range. 

This formula shows that the effect of the taper is to increase 
the fundamental resonant frequency regardless of the taper ratio, but 
the value of the second natural frequency can be increased by choosing 


configurations for which a > l.S b£. 

h 

o 

The first order solutions for the normal modes are found to be: 






It can be easily checked (by looking at the order of magnitude 
of the general term) that the serie on the right hand side of the 
above equation is convergent. 

5. Conclusion and discussion on the limitations of the method. 

Like every pertubation method, this analysis is based i n large 
part on the assumption that z is a small auantity. The term ’'small" 
however has no intrinsec meaning; for perturbation problems z is 
considered small if the specific perturbation method gives consistent 
results. Hence the important problem when dealing with such methods 
is to know what error in the results is associated with a given z . 

If this error is acceptable for the specific engineering purpose for 
which the problem is solved, then e can be assumed small; if on the 
contrary, a sufficiently large number of terms in the perturbation serie 
does not give a satisfactory result, then the value of z does not 
justify the validity of that perturbation method. 

As an illustration of the above, the method is applied to tapered 
cantiliver beams for which the theoretical results are given in the 
paragraph 3 of this Appendix. The accuracy in the first three 
natural frequency is based on the comparison with analytical results 
of Reference (22) . 

One- term perturbation solutions are compared to the exact solu- 
tions for different values of z , and the percentage errors are 
plotted versus e for each of the three frequencies in Figures 48, 

49, 50. 

Figure 48 shows that the one-term solution for the fundamental 
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frequency is really accurate only for e small compared to unity. 

For e = 0.5 the percentage error in the frequency is 6.03% and in- 
creases relatively fast to 9.21% for e = 0.6. 

If a more accurate solution is needed for the fundamental 
frequency of cantilivered beams which such taper, then it is 
necessary to proceed to a second order analysis, or to higher order 
if needed. 

However, as it can be seen in Figures 49 and 50 when the second 
and third natural frequencies are sought, even a one-term solution is 
highly satisfactory by this method. 

For values of e as high as 0.7 there is only 5% error for the 
second natural frequency and only 2% error for the third natural 
frequency. 

This perturbation analysis does not present loss in accuracy 
for higher modes and therefore seems particularly suitable when the 
determination of the higher harmonics is desired. 
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APPENDIX II 


ANALYSIS OF MODAL DAMPING BY COMPONENT MODES METHOD 


USING LAGRANGE MULTIPLIERS 


1. Introduction 

In a rencent note Hallquist and Synder (Reference (42)) formally 
analysed a linear, damped vibrating system with arbitrary support con- 
ditions using Lagrange multipliers. This method is dependent upon 
determining the coefficients of the Rayleigh dissipation function 
for damping. Most commonly, in fact, what can be determined experiment- 
ally are the damping coefficients relative to the classical undamped 
normal modes of a structure. In their recent paper (Reference (43)) 

Kana and Huzar developed an empirical method for predicting the modal 
damping of a Space Shuttle model by means of damping measurements 
performed on the individual substructures. As an alternative, presented 
here is a theoretical analysis which extends that of References (3) 
and (4) to the determination of the modal behavior of a damped linear 
elastic structure of arbitrary complexity by a component mode method. 

2. Theoretical Development 

We conceptually disassemble the structure into N simple components 


for which we know for each mode: 

a) the generalized masses M^ n ^ j = 1,2, 


; n = 1, N 


b) the damping coefficients Z- j = 1,2, .... ; n = 1, ....N 

c) the natural undamped frequency j*l,2,... ; n = 1, N 

For the total system: 
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Kinetic energy: 


CO 

J_ 2 Z 
2 n=l j=l 


M. 

1 


(n) ' 2 (n) 

4 j 


Potential energy: 


II 


N 
7 . 

=1 j = 



M. ^ 
1 




( 243 ) 


(244) 


Damping dissipation function: 


D 


± 2 f 2 Z . (n) (o. Cn) M. Cn ^ q. Cn)2 

2 n=l j = l J .1 1 


(245) 


Here we assume that the damped and undamped modes of components (though 
not necessarily of overall structure) are indistinguishable. 
Interconnecting conditions between the components : 


f 

r 


N 

Z 

n=l 





r=l , 


.R 


(246) 


Lagrangian : 


L = T + D 


R 

U + Z \ f 
— i r r 


where T, D and f are given respectively by equations (243), (244), 
(245) and (246). 
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Lagrange’s equations are thus 


M 


Cn) _ 


, (n) [ qj' n) * 2 C . tn V (n) q/ 1 ” ♦ qj (n) 1 ' | =1 Vr)"" ’ 0 (247) 


Assume time dependent motion of the form, 
q_. (n) Ct) = q^ n ) e(-8J2 + iSl /l^ 2 )t 

X (t)= X e(-8fl + in /H3 2 ')t 
r r 


The Eigenvalue equation is then obtained from (246) and (247) , as 



- 0 


where 


pq 


P Cn) R Cn) 
2 2 B Pi e pi 


N 
E 

n=l j=l 


M.™ { [n 2 (28 2 -l)+2?.^io.^ftB+w.^ 2 ] + 

1 1 3 3 1 


2 in /r-6 2 '(c^ n ^u. ^ - ne) } 


(248) 


For the case where the damping coefficients of the components 
are small for each mode, a first order approximation of equation (248) 


leads to 
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A 

pq 


N 

00 

E E 
n=l i=l 


B [u. 2 -« 2 ) 2 -2iS2 CC - ] 

pj qj Lj : _j 3 i 

m .™ Cm .^ 2 - a 2 ) 2 
3 3. 


(249) 


Thus to first order, for small damping, the natural frequencies, 
ft, are the same as for the undamped system. 


3. Example : 

The utility of the method is now illustrated by a simple example 
which is that of a structure composed of two identical beams connected 
at their middle. In particular, for doubly symmetric vibrations, the 
single constraint condition reduces to the equality of the displacements 
at the mid-point. The natural frequencies and damping coefficients in 
each mode are given by (from equation (248) for N = 1) 


o» 2 2 2 

z (0) (to. V) - 0 

j=0,2,4 3 — 

(w . 2 - ft 2 ) 2 

3 


2 V(°) ft (e; . w.-ft8) - 0 

j=0,2,4 ] 3 _ l 

(U). 2 - ft 2 ) 


(250) 


(251) 


where represents an unconstrained beam mode. 

The solutions ft. = u> ; 6. » £. corresponding to A = 0 are phy- 

J J J J 

sically valid and represent modes in which the beams vibrate symmetrically 
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together, without any constraint force required, 
Other solutions for X ^ 0 are given by: 


2 ip . 

3 = 0 , 2, 4 -1 


CO) 


= 0 


W. 


> 


■ s * (o) a ( c,w. - w ) = o 

1 j=0,2 ,4 -2 2-J 


I 


L 


2 2 

Go. - cr) 


( 252 ) 


(253) 


One can show that there is one root of the frequency equation (252) 
in between each two consecutive unconstrained natural frequencies <o . 

A two-term analysis for the fundamental bending mode gives: 

0 « 0.758 w 2 

8 = 0.758 z 2 

In this mode, the overall structure is more flexible than its 
individual components and contains less damping. 


4 . Discussion and Conclusions: 

For any arbitrary structure the natural frequencies and modal 
damping can be predicted from the knowledge of the natural frequencies 
and damping coefficients of each of the components. All parameters 
appearing in the analysis can be measured experimentally and the 
technique is simple and offers physical clarity. Further study needs 
to be given to the damping which may arise at the connection of the 
components per se. 
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APPENDIX III 

FREE VIBRATIONS OF A BEAM CLAMPED IN ITS MIDDLE BY A 
RAYLEIGH RITZ METHOD USING LAGRANGE MULTIPLIERS. 

1 . Introduction ; statement of the problem. 

The present analysis is mainly intended as background for the studies 
of the limiting cases considered in Chapters III and IV. The problem 
considered is that of a beam of length 2L, free at both ends and clamped 
in the middle, as illustrated below: 



M: 

half mass of the beam 


El: 

bending stiffness 


^(C): 

unconstrained beam modes 

!<_ £ <i 

k. : 

corresponding non-dimensional 

frequency parameters 


W (5 , t) :bending displacement of the beam 

The kinetic and potential energies of the beam are: 

1 *2 

T = 1 M { w (£,t) d£ 

2 


( 253 ) 
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u = i ei { r , 2 (c,t) d? 

2 L 3 ' 

2 . Modal analysis using even modes of an unconstrained beam. 
The modal displacements are sought on the form: 

W(S,t) = I <^(t) 4 n CO 
n=l,3 

and thus (253) and (254.) become 

T = 1_ M * q 2 (t) ; U = i HI t k* q 2 (t) 

2 n=l,3 n 2 l 3 n=l , 3 

L iimiiib PI U I MU I ■ - JL "• ' ~ 

by virtue of the orthonormality of the <|>. (£).' Because of 
of even beam modes, the only condition for the beam to be 
the middle is : 

W(0,t) = 0 

or, by use of (255) 

? q (t) A CO) = 0 
n=l,3 

The corresponding Lagrangian is : 

EI \ 4 q„ 2 Ct) ♦»q„(t)» n C0) ) 

21 

and the subsequent Lagrange's equations are: 


JC = ? (M q 2 (t) - 
n=l , 3 2 


( 254 ) 


(255) 


(256) 

the choice 
clamped in 


(257) 
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• M q (t) + El k 4 q (t) - (0) =0 n=l,...« 

L 3 


• £ _, , %<*> V°> ■ 0 

n- 1 j j 


Harmonic solutions are sought in the form: 


1 C El ) 
■ On e ML 3 


•h 


(Ot 


X(t) = X El 


i ( El ) 
ML 3 


wt 


and thus : 


On = - * V 0) 


2 v 4 
co - k 

n 


(258) 


1 , , x *n™ - 0 
n=l,3 

2 v 4 

to - k 

n 


( 259 ) 


Equation (259) leads to the frequency equation: 
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(260) 

and equations (258) determine the generalized coordinates for each 
natural frequency. The mode shapes are then given by equation (255) . 
3. Modal analysis using antisymmetric modes of an unconstrained beam. 
Only odd free- free beam modes are used in the modal serie : 

W(C.t) = Z qft) * (?) -1<£<1 (261) 

n=2, 4 

and equations (253) and (254) become: 

(262) 
(263) 

The zero-displacement condition at the middle of the beam is 
already satisfied by the coice of the antisymmetric modes in the 
analysis; however the beam has to be constrained so that: 

1 W* (o,t) = 0 
L 
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or using (261) 


op 

1 z 


L n=2,4 


q^Ct) ** n (0) = 0 


(264) 


If y is the Lagrange Multiplier corresponding to the above con 
straint, the equations of motion are: 


• M a (t) + El k q ft) - y <j>» (0) * 0 n=l,...» 

i* » n n n 

L L 


* 2 q ft) (0) = 0 
n=2 ,4 — — n 


Harmonic time dependent input is assumed as follows: 

i ( El ) 1/2 tot 
ML 3 

q n (t) = q n e 

i ( El ) 1/2 tot 


u(t) = y El e 
L 2 

and hence the Lagrange's equations are equivalent to the system: 


y <j>' n (0) 



Z 

n=2 ,4 


y 


V n 2 (0) = 0 



n 


(265) 


( 266 ) 
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From the eigenvalue equation (266) we obtain the frequency 
equation : 


00 

£ 

+ ’ 2 
y n 

(0) = 0 

n=2 ,4 

! 

CNJ 

3 

k 4 



n 


and equations (265) and (261) determine the corresponding modal 
displacements . 

4. Discussion 

Although a more elaborate analysis is presented in Chapter III, 
a two-modes approximation of the frequency equation can give some 
insight on the rate of convergence of the fundamental frequency 
according to the type of modes used. 

A two-terms expansion of equation (260) gives: 

_ 1 _ ♦ *3 2 W * 0 

2qj 2 (o 2 -k^ 4 

and yields the following solution for the fundamental frequency: 

to = 3 . 553 

with the values of <^(0) and given in Reference (35). 

Hence, a two-symmetric-modes analysis is seen to produce an error of 
about 1% in the fundamental natural frequency. 

If two antisymmetric modes are used in equation (267), the fun- 
damental natural frequency is given by 
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w - K 2 (1 + 2 cj>' 2 COD 

3 4 

and with the values of and ^'^(O) given in Reference (35), 
yields the numerical value : 


u = 4.709 

In the analysis using antisymmetric modes, a two-modes approximation 
gives a 5.4% error in the fundamental frequency and thus we can 
anticipate better convergence using symmetric modes. 
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APPENDIX IV 

MODAL ANALYSIS OF A FLEXIBLE AIRPLANE WITH LOW ASPECT RATIO WING AND TAIL 
1. Introduction 

The present analysis extends the study of Chapter III to the 
more general case of a flexible airplane with a wing and tail of 
finite or low aspect ratios. 


Such a configuration is illustrated below. 

Vr 



In this analysis the fuselage has a one-dimensional beam repre- 
sentation with freedom in bending and torsion, while the wing and 
tail are represented by uniform plates undergoing bending in the x 
and y - directions . 

Each one of the components is given a linear treatment, geo- 


metrically and elastically. 
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The rigid body dynamics of this structure are similar in concept 
to those investigated in Chapter III. 

Because of the geometrical symmetry, pitching and rolling motions 
can be considered independently of each other. Yaw-type of motions are 
not in the purpose of this analysis. 

2. Modal analysis for motions of the rolling type . 

Such motions are characterized by pure twisting of the fuselage, 
and antisymmetric behavior of the wing and tail about the x axis. 

The approach is based upon the use of free-free twisting modes 
for the fuselage and unconstrained plate modes with antisymmetry 
about the x axis, for the wing and tail. More specifically, the 
latter modes are of the form: 


* (x,y) = Z Z A.™ X (x)Y (y) (268) 

n k=l m=2,4 Km k m 

where 

X. and Y are free-free beam modes, and the "weighting" coefficients 
k in 

are found together with the natural frequencies by solving 

the eigenvalue problem: 



aP*) + s^K,, 4 ) + s 2 2 A> J [v(oi. (0 ) 

km v k m .. „ 0 . lr L v Hcl rm lk mr J 

1=1 r=2,4 


00 00 


(n) 


2 (l-v)a, -a ] = K 2 fS 71 '* 
v J kl mr J n km 


( 269 ) 


In the above equation: 

k. : natural frequency parameter of free-free beam. 
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s : ratio length/width plate 

v : Poisson ratio of plate 

1 

un = ' X” (0 X,(0d5 

■*-j 1 - 


CL. . - 

11 


! Y » 


•1 


X', COX', COdC = a.. 


1 


Ji 


The angular displacement of the fuselage 0(x,t) is sought given 
by the following seriet 


0(x.t) = l P„Ct) 0_ (x) 


n=l 


n ' n 


- L < x < L 


(270) 


0 n (x): free-free modes in twisting 

and the bending displacements of the wing and tail W^(x,y,t) and 
(x 3 y,t) are sought in the forms: 


W 1 (x,y,t) = Z (t) (x,y) - L x £ x <_ Lj (271) 

n=l 


- lj <y < 1. 


n=l 


W 2 (x,y.t) = f ^ q n < ' 2 ' ) Ct) ® n ^ 2 ^(x,y) - L., <_ x <_ L„ 


2 - - 2 




(272) 


If the wing and tail have different ratios length/width, then: 


♦/^(x.y) J* 'f n ( ' 2 ' ) (x f y) 


The total kinetic energy of the structure is : 
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T = 1_ [ ' I G 2 (x,t) dx 

2 -L -L, -1 


U 1, . 2 

+ 1 1 Wj (x,y,t) dm 1 + 


1 


1 2 1 2 w 2 (x,y,t) dm ] 

■ L 2 -»2 


( 273 ) 



where I is the mass moment of inertia per unit length of the fuselage 
about the x-axis. By sybstitution of (270), (271), (272) and use of 
the orthonormality of the component modes, equation (273) becomes: 


(274) 


In the above equation and are one fourth of the masses of the 

wing and the tail respectively, i.e. relative to half-length and half- 
width of these components . 

The potential energy of the airplane is: 


U = 1 f ,l GJ(30 (x,t)) 2 to ♦ I °i 'V* 

2 - L 9X 1=1 - L i -h 9x 2 9y 2 


2 2 2 2 

. r , 31. * 3%. 3 w . , • i 
+ 2 (1-v^ ) [ ( i ) - i i ] dxdyj 


3x3y 


2 2 
3x 9y* 


(275) 


In a first step equations (270), (271), (272) are substituted into 
(275), then use is made of the classical beam equation in twisting, and 
of equation (269) applied respectively to the wing modes and the tail 
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modes . 

After lengthy algebraic manipulations, the potential energy is cast into 


the form: 


» ? o ? 2 W 2 Ci) 2 

U = I 2 ( GJ (n-1) 1 i P n (t) + 2 D.L. K n (t) ) 

2 n=l L 4 i = 1 <? 

1. 


C276) 


Finally, the continuity of the displacements and slopes between the com- 
ponents has to be enforced. The slopes in the x direction and the 
displacements are already zero because of the choice of antisymmetric 
modes only, in the y-direction for the plates; as for the slope con- 
ditions in the y direction they are: 


aw 

0(x,t) = 1_ (x,0,t) - Lj <_ x <_ 

3y 

0(x,t) = 5 ^2 (x,0,t) - L_ <_ x <_ L- 

3y 

The variational formulation of these conditions is: 


1 

' * C1) (x,t) [0(x,t) 

1 


3 ^1 (x,0,t) ] dx 

3y 


(277) 


(2) r ; 


3W, 


-L„ 


(x,t) [0(x,t) - 2 (x,0,t) ] dx 

9y 


(278) 






181 


where and A^ are the corresponding Lagrange multipliers. In 

the present case A^(x,t) and X^(x,t) are the generalized bending 
moments per unit length, developed along the widths of the wing and 
tail respectively and necessary to keep the continuity of slopes. 

The integrals (277] and (278) are considered in a Dirac 
approximation; the alternative procedures discussed in Chapter V are 
prooved to be equally suitable but give rise to more elaborate com- 
putations . 

For this purpose, the following notations are introduced: 


(x 



r=l , R 


partition of the interval [-l^, 1^] 




r 




(x 


( 2 ) 


r 


) 


r=l,R 


partition of [~1 ? , 1 0 ] 




A^^ (t) , A^ (t) : discrete values of A^ and A^ at the corre- 

r r 

sponding points. 

Then, the integrals (277) and (278) are approximated by: 


R 

Z 

r=l 


X<« 


I © (5' r ,t) 


aw. 


9n 


(5 r (1) ,0,t) ] 


(279) 
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E A ^2) [ 0 (5”>t) - 1_ ^2_ (C^i^t) ] (280) 

- 1 r " X 1 3n 1 

which can be written in modal form by use of (270), (271), (272) as 
follows : 


R 


( 1 ) 


E X'- ■’ E (P (t) 9 ( V ) 
, r , n l 7 n k r 

r=l n=l 


- q n ( 1 ) (t) s* n ll X ( Vo)) = 0 (281) 


'1 3n 


£ 

E X ^ E CP (t)0 (£’ ') 
, r , ^ n v J n r 

r=l n=l 


_ q n (2) (t0 9\ (2) (€ r (2) ,0)) = 0 (282) 

1 2 3P 


The subsequent Lagrange's equations are thus: 


ILP (t) + GJ (n-l)Vp (t) 
n l 4 n 


e (x (g’) 

, r n r 
r=l 


x (e”))=o 

r n v ^r 


• * Vi q/'ht) ♦ E . 0 ,i . 1,2 


r=l 1 


i 9n 


together with the continuity conditions (281) and (282) 

The following hamomic behavior is assumed in the analysis: 


p n (t) 


^i/ ^i=l,2 J 


n 


^i*l p 2 J 


i 1 (D^pp 172 at 


1 


x e 
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i JL CD, !p,) 1/2 cot 




T (i) D I. 

^ x 11 x e 


where is the mass density of the wing per unit area. 

The following non-dimensional parameters are defined: 


1 


M,l, 


M, 


y.= 1 ; y = l ; u = 1 1 ; = 1 


'2 — ’ M l 

L 2 


M, 


GJ 


1C 


2*2 


*1 = 


l 1 ! 


I 1 ! 


- ( 1 ) 2 ..2 2 

F J = to - (n-1) 7T tjj 

" 4 1 


F (2) - „ 2 - K C1)2 
n n 


„ ( 3 ) 2 ( 2)2 . 

F = (0 - K ip„ 

n n 2 


Then the generalized coordinates are expressed in function of A ^ 

r 

( 2 ) 

and A^ from the equations of motion as follows: 
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P 

n 


R 

Z 

Vi Z± 


(i) 


( 2 ) 


e te' ) + x 

r r 


e_ (C') 

n r 


a) 


n 


(283) 



R 

E 

r=l 


_(i) 

X 


(i) 

94 > 

n 


p (i*l) 
n 


(i) 

U r . 0) 


i=l ,2 


(284) 


Substitution of the above equations into the constraint conditions lead 
to the eigenvalue equation and then to the frequency equation 


det (4 (u) ) = 0 


(285) 


with the following representation of the eigenvalue tensor: 
. 1 < r < R ; l<v<R 


V M ! , <Vl 9 n«'r e n«V * 

n=1 — r® — 


(1) 

94 > 

n 

9n 



0) 


( 1 ) 


. . 9$ 

■ n 

9n 



0) ) 


F 


12> 


n 
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.l<r<R ; R + l <v<2R 


A 

rv 


Cw) 


00 

z 


n=l 


Y l u l 


0 Cc" )0 (C* 

n r n v 


F CD 


.R + l<r<2R ; R+l<v<2R 


3$ 


( 2 ) 


A (u>) 
rv 1 


z 

n=l 


( Vl 


e C€") 

n r 


e fc") 

n v 


Y 2 y 2 


n 


sn 


Cl) 


C£^ 2) ,o) 

C C2) 


3$ 

r 

9n 


C2) 


(E 


n 


3. Modal analysis for motions of the pitching type 

Such motions are characterized by pure bending of the fuselage and 
symmetric behavior of the wing and tail about the x axis. The procedure 
is similar to that of paragraph 2, except for the facts that free-free 
beam modes in bending are used for the fuselage, and that the uncon- 
strained plate modes are now based on even A fv) . 

m 

Hence, equations Cl) and (2) are still valid formally, but the sub- 
scripts m and r are 1,3,5, 

The bending displacements of the three components are sought 
given by the series: 


W(x,t) = Z q^ft) <f> n (x) 
n=l 


- L < x < L 


C.286) 
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W.(x,y,t) = ? q/^Ct) $ n (l) (x,y) -Lil x l L i 

n=1 ( 287 ) 

- 1. < y < 1. 
l — — i 

i = 1,2 

and thus the kinetic and potential energies of the structure are: 


-T 


II (Mq 2 (t) ♦ I M q_ Cl) (t) ) 
2 n=l 71 i=l 1 n 


( 288 ) 


U 


1 

2 


? , <Ei k „ 4 % 2 (o * 

n =l l 3 


2 

Z 

i=l 


D.L. 

l l 

l . 3 

l 


(i)' 


(i) 


Ct) ) 


( 289 ) 


where M is half of the fuselage mass, and El its bending stiffness. 

The above equations are derived in a procedure similar to that 
of the paragraph 2 and therefore the details of the computations are 
omitted . 

Because symmetric behavior of the wing and tail has been assumed 
in equation (1), the continuity of slopes in the y-directions between 
components is automatically satisfied; the continuity of the dis- 
placements and slope in the x-direction at the connections between 
components is expressed by: 


W (x,t) = W^x.O.t) 


( 290 ) 
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9W 


(x>t) 


9x 


aw. r 

1 0,0, t) 

9x 


C291) 


for - L. < x < L. and i = 1,2 
i — — l 

By following a similar procedure to that of paragraph 2, the 
following expressions are obtained for equations (29)) and (291) in the 
variational form: 


R 

Z 

r=l 


X ™ 


2 (q. 

n=l 


n 


(£') - q < - 1 ' ) (t) 
r n n v J 


4 ^ (s ^ ,0)) = 0 
n r * JJ 


(292) 


£ >«>[ fa ft) V (5- ) . % a) ft) 80 (I) C5 (1 ),0)) . 0 (2931 

- r in n v r 

n=1 T~ 4 « ■ 


1 I 1 ^n^n^P “ q n C2)c ° $ n (2) O r C2) ,0) ) = 0 (294) 


R 

Z 

r-1 


? ( q n ft) 4>' (C ') 

- n=l T n ' r 


^ ( 2 ) ct) ™ n (2) a r ( 2 ) ,o)) 

K 9S 


= 0 (295) 


Also, the harmonic analysis in the Lagrange's equations assumes similar 
harmonic time dependence. 

The following non-dimensional parameters are defined here: 
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Vi! 


= ~ ’ ^ ; v 2 = ^ 


*1 = 


El 

ML' 


^2 = 


D, 


Y 2 1 2 


Vl 


Vl 


F ™ - * 2 - k \ 

n n 1 


C2) , „ 2 . K (11 2 


F ' ' = to 
n n 


F < 3 > . „ 2 - K P) 2 *- 
n n 2 


The equations of motion are first solved for the generalized coor- 
dinates : 


% 


R Cl) _(2) _C3) _(4) 

Z A 4> (OM $'(£') + * <f> (E M )+M <1>’(5") 

. r 1 r n^r J r r n K r J 1 r T n r 

Uj r=l 


CD 


n 


R 


-a, 

a ' = r=l 

'n 

„ C2) 


-C2) 3$ CD 

r 

3C 
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-(2) X r (3 ’ t „ C2) «r C2) - 0 > *. 5 2 X r 4) ^ ^ ^ 

V ‘^ 2 ?! ^ 


(3) 


n 


and then the eigenvalue equations are found by substitution into 
(292) - (295) . The eigenvalue matrix is of order 4Rx4R, symmetric, 
and its elements are: 

, l<r<R ; l<v<R 


A 

rv 


(w) = 


Z 

n=l 


(V^vcp^nte 1 ; . * n o) te r cl) ,0) * n (1) « v (1) ,0) ) 
p Cl) F CD 

n n 


. 1 < r < R ; R + 1 < v < 2 R 


3 <f n a) rr C1 > OH 

'A (to) = Z ( n,M’ C5’ )<f> (V ) — r—- U r ,lJ n C V 015 

rv v ' , 1 1 T n^ r J n v + 3£ 

n=l —2 


( 1 ) 


( 2 ) 


n 


. 1 r < R ; 2R+1 < v < 3R 


A ((d) = £ Pl ^n^’r^n^V 
w n=l 1 


Cl) 


n 


. 1 < r < R ; 3R+1 < v < 4R 
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A 

rv 


(to) 


£ . Vi ♦’„ t5 'PV 5 V 

n=l 


Cl) 


n 


R+l < r < 2R ; R+l < v < 2R 


‘rv'” 5 


d) m .* Cl),, (2) ... 

3« < E C l \o) fn. «r *°» 

2 , C^, 2 *' n cc r )V n C5' v ) tTr n 

n=l — 2 


(1) 


C2) 


R+l < r < 2R ; 2R+1 < v < 3R 


A C«) = E y 6 ‘•’n^'r^'n^ V 

rv v J .11 

n=1 F (i) 

' n 


. R+l < r < 2R ; 3R+1 <_ v <_ 4R 


A 

rv 


M 


z 

n=l 


Vl 


V n V>» r ^T n ».S, V J 


F 

n 


CD 


. 2R+1 < r < 3R ; 2R+1 < v < 3R 


A 

rv 


Cw) 


E 

n=l 


(y . (f> (£" )<{• (5") 

v 1 n v r n v 


v 2 V C2) cc r (2) . 0 )»n t2) a v t2 !o)) 



F C3) 
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. 2R+1 < r < 3R 

• A rv Cw) f Cy 

p CD F C3) 

n n 

. 3R+1 _< r <_ 4R ; 3R+1 £ v <_ 4R 

2 2 3t n (2) C m 9 *n t2) f 5 v t2) ' 0) 

4 rv (“) = ^ . V J^2 JT C r- 0 2 E 1 

F CD F C3) 

n n 

Non-trivi’al solutions for the Lagrange multipliers exist if: 

det ( A(w) ) = 0 

i.e., for the natural frequencies of the structure. 

4. Discussion of the convergence of the method. 

The main sources of error for the natural frequencies are the use 
of a finite number of modes per component and the choice of a finite 
number of points at which the constraints are applied. 

An obvious way of solving the problem would be to give to both 
parameters very large values form the beginning, but then the computation 
may be too time consuming for the accuracy desired. A better approach 
is to aim at optimizing the values of the two parameters based on pre- 


; 3R+1 < v < 4R 

C2) °» 

“ 2 { 2 — ’ 
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vious information acquired in the dissertation. 

As a first remark, it can be noticed that the bending or torsional 
free- free modes of the fuselage are analytically known functions and 
are used only in elementary operations such as summation. Therefore, 
a large number of these modes can be used without any loss in accuracy 
or significant time increase. 

The problem is much more comnlicated, however, for the plate modes. 

The unconstrained modes are obtained by numerical techniques applied 
to equation (?69) , and therefore neither the A^ ^ nor the are 

rigourously known. As a consequence, it is desireable to avoid using a 
larger number of modes than is strictly necessary, in order to minimize 
the initial errors propagating to the final result. 

The minimum number of modes necessary is, of course, determined 
by the number of natural frequencies desired and by the number of con- 
straint points considered. For instance, if five constraint points are 
considered along the width of the wing and the first three natural fre- 
quencies are desired, at least ten degrees of freedom in the x-direction 
should be used in the analysis, for a reasonably good accuracy. Obviously, 
a suitable number of modes in the y-direction has to be considered as 
well . 

The purpose of this discussion is to suggest a systematic approach 
for the optimum choice of plate modes versus points of constraint, for 
the wing and the tail. The procedure is based upon convergence in a 
limiting case, i.e. the case in which the fuselage behaves r’gidly and 
therefore the constraints at the connections are most critical. 
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The analysis of the Chapter V shows that by using five modes 
X^C*) and five symmetric modes Y^(y), an unconstrained plate with 
aspect ratio 2 can be considered clamped in the middle of the first 
three natural frequencies with as little as seven constraint points. 
The fundamental natural frequency is then predicted with less than 
0.1% error, similar precision is obtained in the second natural fre- 
quency, and an error of only 0.4% is obtained in the third natural 
frequency. If the plate has higher aspect ratio, an even smaller num- 
ber of points should suffice for comparable accuracy. Therefore, 
setting R=7 in the analysis of the pitching modes should be very satis 
factory for at least the first three natural frequencies. The order 
of the eigenvalue matrix would then be 28, which is not unreasonable 
for computations of such type. 

The analysis of Chapter III indicates that for motions of the 
rolling type, a number of twenty antisymmetric modes should be used 
in the y-direction in order to obtain same accuracy in the natural 
frequencies. This disadvantage is compensated by the fact that in thi 
analysis the eigenvalue matrix would be of order 14 only and hence the 
total computation time should not be higher. Finally, it is also suit 
able to check the unconstrained plate modes for orthonormality at the 
beginning of the numerical analysis in order to avoid possible errors 
comming from an incorrect numerical assumption. 
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APPENDIX V 

EFFECT OF THE ELASTIC SUPPORTS ON THE NATURAL FREQUENCIES 
OF THE EXPERIMENTAL MODEL CONSIDERED IN CHAPTER III 

1 . Introduction 

As discussed in Chapter III, the effect of the supporting springs 
on the resonant frequencies of the rolling type can be completely- 
eliminated by a perfect allignment of these supports along beam (1) . 

In the present study this is supposed to be the case and only the 
corrections to the frequencies of the pitching type are considered. 

The springs are supposed to be massless and perfectly elastic. 

The analysis is developed by a Component Modes Method using La- 
grange Multipliers. The notations of Chapter III are used in the 
present analysis; some new notations are introduced below: 
k , k 2 : elastic constants of the two springs (the experiment was per- 
formed with k^=k 9 but generally this might not be the case) . 
Z^(t), Z ? (t): longitudinal rigid body translations due to the elongations 
of the springs. 

2 . Theoretical development 

For the total system in which the supports are included, the 
kinetic and potential energies are: 
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(297) 


oo oo f i 

£ <L(tH (£.) - £ O j (t)i (0) = 0 i=l >2 C298) 

n-1 71 n 1 n=l ,3 n n 


£ (Q-Ct)*' (5.) + P n (i) (t) 0 (0) ) = 0 i-1,2 (299) 

n=l 

L 


( £ 0 n (t)4> n (-l) ) - Z 2 (t) = 0 (300) 

n=l 


( ? (^(0^(1) ) - Z 2 (t) = 0 (301) 

n=l 

Equations (298) and (299) are found in the analysis 
developed in Chapter III, equations (300) and (301) express the con- 
tinuity of the displacements between the beam (1) and the supports 
placed at the tips. 

The Lagrange's equations of the system, are: 
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. M Q^t) . Ell k n 4 ^ft) - i 1 +„te 1 ) - V n C5 2 ) - ^ ♦W 
L 3 L 

- n U 2 ) " X 5 V" 13 ” X 6 K a) = 0 C 302 ) 

L 

• M i % (i)w + — 3 L k n 4 % Ci) ^ + X i^ n (0) = ° C303] 

1 3 
1 

for i - 1,2 and n=l,3,5,... 

• Vi '’n (i)(t) * — 1 O- 1 ) 2 ! p n (l) t t ) - X 2.i 9 n (0) = 0 t304) 

1 . 4 
1 

for i = 1,2 

. k.Z.(t) + X 4+ . = 0 C305 

for i = 1,2 

. equations (298] - (301] 

The physical interpretations of the Lagrange Multipliers A^, 

Aj X 4 has been given in Chapter III; A,. and A g are the internal 
forces developed at the tips of beam (1] when connected to its supports 
(this physical interpretation is particularly obvious from equations 
(305)]. 

The equations (302) - (305) are solved for the generalized coordi- 
nates, when harmonic time dependent motion is assumed as in Chapter III. 
In particular: 
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The resulting tensonial eigenvalue equation is 


— -> ->■ 
AO) . X » 0 


(305) 


where 

A : ( ^ 1 ' ^2 ’ 



and AO) is a sixth order symmetric tensor with the following matrix 
representation : 
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For given values of kj and k^, 
supported structure are solutions of 


the natural frequencies of the 


det (AO) ) = 0 

3. Remark 

The foregoing analysis can be easily extended to the case in 
which the structure is resting on two beam-like elastic foundations 
perpendicular to beam (1). 
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In such a case, ( k., Z. ) . , - are substituted in the analysis 
by the free- free frequency parameters and combinations of normal modes 
and the additional generalized masses terms appear in the kinetic energy. 



- 202 - 
APPENDIX VI 


EFFECT OF THE MASSES AND INERTIAS OF THE ENGINES ON THE NATURAL 
FREQUENCIES OF THE LOCKHEED L-100 MODEL CONSIDERED IN CHAPTER III. 

1 . Introduction 

In the analysis of the L-100 model considered in Chapter III the 
effects of the engines were ignored. This omission was intended only 
for the purpose of computational simplicity and the present analysis 
indicates the modifications to be made in the previous theoretical 
study in order to take into account the masses and inertias of the 
engines. 

The Lockheed L-100 Hercules is equipped with four 4,050 eshp Allison 
T 56-A-7A turboprop engines whose characteristics are given in Ref. (32) . 

A schematic representation of the structure is illustrated below: 



Because the geometrical symmetry of the structure is preserved, 
motions of the rolling and pitching type can be studied independently. 
2. Effect of the engines on the rolling natural frequencies 


Additional notations to those of Chapter III are now introduced: 
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n 1 , ri 2 : location of the engines along the span of the wing 

m_ : mass of each engine 

C 

Z 1 (t),Z 2 (t) : plunging displacement of the engines 

I : moment of inertia of each engine about its axis of symmetry 

a 1 (t),a 9 (t) : angular rotation of the engines about the axis previously 

defined. 

The presence of the engines accounts for an additional terra in the 
expression of the kinetic energy given by equation f41) ; 

T engines ’ * V > < 3 °« 

while the expression for potential energy given by equation (43) remains 
unmodified. 

Additional constraint conditions expressing the continuity of dis- 
placement and slope at the connection points along the wing span are: 

( Z q^(t)* (rij) ) - Z 1 (t) = 0 (307) 

n=2,4 

( E q^^ftH (nj ) - Z (t) s' 0 (308) 

n=2,4 * z 


( £ J (t) * t n (n 1 ))-a 1 (t) - 0 

n=2 , 4 " 


00 Q ^ ftl 

( 2 % + f n (nO) - ct (t) = 0 

n=2,4 * 1 
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(309) 


( 310 ) 
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and ( A. ). , A are the corresponding Lagrange Multipliers, or 
generalized forces developed at the points of constraint. The modified 
Lagrange’s equations are: 

. I L P n (t) + GJ (n-l) 2 1 2 P n (t)- X 3 0 n (? 1 ) - * 4 9 n (£ 2 ) = 0 (3U:i 

Li 4 


q^ 1 -* (t) + ^^1 q^ 1 -* (t) + ^ n (0) - A^ n (Hj) - 


A + (n ? ) - A 3 p W " A 4 ° 

2 n 2 ' — — 


(3121 


q „ ( 2 ) (t)^ kn V 2 ’ w -o 

l 2 *2 


(313) 


2 Z. (t) + A = 0 
1 1 
E 


(314) 


2m E Z 2 (t) + A 2 = 0 


(315) 


2 I x ctjft) + A 3 = 0 


(316) 


2 l x « 2 (t) . A 4 = 0 


(317) 


We assume harmonic motion as defined in Chapter III and introduce 
the following new non-dimensional parameters: 


* M iV 

= 1 1 
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p „ and \ (2) 


remain as defined in Chapter III, but: 
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Substitution of (318), (319), (320), and q^ ^ J into the constraint 
equation leads to the eigenvalue matrix and hence to the frequency 
equation for the new configuration. 

The natural frequencies obtained in Chapter III are found as 
particular solutions of this frequency equation, in the case where 
n l = n 2 - 0, by factorising out the non-identical ly zero determinant. 

3. Effect of the engines on the pitching natural frequencies 

The problem is slightly more complex in this case because additional 
degrees of freedom have to be allowed for the engines, namely rotations 
about the y-axis. 


Let 
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y : mass moment of inertia of each engine about the y-axis 

: angular rotations of the engines about the y-axis. 

The presence of the engines corresponds to an additional term in 
the kinetic energy: 

T = (Z x 2 Ct) + Z 2 2 Ct)) + Ct) 2 *a 2 2 (t)) + 

IyCe^Ct) + 6 2 2 (t)) (321) 

but the potential energy of the overall structure remains unchanged. 
The new additional constraint conditions are : 


C £ V " z i Ct) = ° 
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n-1 1 n 1 


e 2 (t) = o 


(5 S Iil*' n (n 2 )) 
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(323) 

(324) 

(325) 

(326) 

(327) 


From this point on same procedure as previously applies, and 
therefore the details will be omitted. 
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4. Conclusion 

The expected changes in the natural frequencies due to the presence 
of the engines can be discussed on a qualitative basis even before 
performing the numerical computations. 

Because mass and inertia has been added to the structure, the 
new natural frequencies are expected to be lower than those computed 
in Chapter III (ef. Ref. (45), p.83) 

In particular, the first natural mode of the structure was shown 
to be essentially a wing mode; consequently, its corresponding natural 
frequency is expected to change by a larger amount than the natural fre- 
quencies corresponding to modes in which the elastic behavior of the 
fuselage or tail is predominant. 
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APPENDIX VII 

EVALUATION OF INTERNAL FORCES (SHEAR AND BENDING MOMENTS) 

ILLUSTRATED FOR A CANTI LIVE RED BEAM. 

1 . Introduction 

The Rayleigh-Ritz Component Modes Method using Lagrange Multipliers 
is characterised by two main factors. 

(1) The use, for each component, of its unconstrained modes for 
determining the deformation in each overall structural mode. 

(2) The identification of the Lagrange Multipliers as the shear forces 
and bending moments at the points of constraint. 

The computational advantages introduced by (1) are particularly 
clear at the stage in which the constraint conditions have to be written. 
In the one-dimensional case of beam structures for instance, the bound- 
ary conditions satisfied by the free-free modes are such that the con- 
tinuity of the second and third derivatives of the bending displace- 
ment is automatically satisfied. This results from the fact that 
W’(x) and W* * * (x) (with the notations used in Chapter II), have zero 
value at the points of constraint, as given by the Rayleigh-Ritz 
expansion. Obviously, the bending moments and shear forces which are 
proportional to Vf' * and W * 1 should not be zero on physical grounds, 
and consequently the problem which arises is that of the convergence 
of W M (x) and W ' 1 (x) as given by the method. This problem is also 
connected to (2) above. 

For the sake of clarity, the author has chosen to illustrate 
the convergence analysis on a conceptually simple case, which is that 
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of a cantiliver beam. In other words, the problem is studied for one 
constraint point. 

2. Convergence of W* 1 and W*'* illustrated for a cantiliver beam. 

ji 

Consider a beam of span L, 0 x <_ L, clamped at the end x = 0, 
as shown in the figure below. 

% — ,,-rr -J * X 

o L 


The treatment of this problem is based upon seeking mode shapes 
as given by the expansion: 

W(x) » 2 q n 4> (x) 

n=l 

subject to the constraints: 

Z q n * (0) = 0 

n=l 

12 ^ ^ n (0) = 0 

L n=l 

By proceeding as stated in Chapters II - IV, the Lagrange multi- 
pliers and generalized coordinates are determined, and consequently 
the functions 

W*'Cx)=2 q„ +" n 00 


(3 


(•3^*0 
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and 


W"'(x) - ? q n <t>”’ (x) 
n=l 

are also determined. 

The variations of W ' (x) along the beam span in the fundamental 
mode is represented in Figure 51; the different curves correspond 
to the number of modes considered in the analysis which is indicated 
along the curve. The exact variation of W* ' (x) can also be seen, 
as given by Reference (21). 

Two fundamental remarks can be made: 

(a) As the number of modes considered in the analysis increases, W* 1 (x) 
to the known exact result except in the neighborhood of x = 0. 

(b) However, the exact value of W' * (0) given by Reference (21) 

W' ’ (0) = 0.8), is only 8.75% higher than the value of the corresponding 
Lagrange Multiplier, obtained in a 40-modes analysis. 

This fact is not surprising, as the value of Xj/X^ (the Lagrange 
multipliers associated with equations £3*^2 and (.3303 is precisely 
that of the bending moment at x=0 (up to multiplication) , i.e. 

Hence, as shown in Figure 51, the method offers a good representation 
of W* ■ (x) as the number of modes increases; moreover, because the 
ratio X^/X^is computed in the process, it is possible to have an 
insight on the convergence of W* ' without even having to actually 
evaluate W' • (x) along the span. For fairly complex structures, this 
fact represents a significant saving in computation time and effort. 
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Figure 52 shows however , poor results for the convergence of 
W' * ' (as expected, the rate of convergence becomes slower as the order 
of the derivative of W increases), but the value of the Lagrange mul- 
tiplier predicts much better the value of W' 1 ' (0) than the direct 
result as it can be seen in the figure. This emphasizes again the 
convenience of using the Lagrange multipliers in the method. 

3. Generalization and conclusions 

The convergence results obtained in the one -constraint point- 
analysis can be readily generalized for the cases in which the iso- 
perimetric (constraint) conditions have to be applied at more than 
one point. A behavior of the type observed in Figure 51 is expected 
to occur in the neighborhood of each of the constraint points, with 
the values of W’ ' and W ’ 1 at the constraint points given by the 
values of the Lagrange multipliers associated with the particular 
constraint equation. 


The following picture may help illustrate this statement. 
This picture is drawn for the model of the rotor blade considered 
in Chapter II, in the fundamental mode. 



* 



The type of behavior shown above can be explained from a mathematical 
point of view by the uniform convergence of the serie giving W(x) al- 
most everywhere, i.e. except at the discrete number of constraint 
points. As a consequence, for example. 
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APPENDIX VIII 

EFFECTS OF FINITE CROSS-SECTIONAL DIMENSIONS ON BEAM INERTIAS AND CON- 
STRAINT CONDITIONS 

For structures of the type shown below: 



where, in contrast to the assumptions of Chapters III and IV, 
the points A and B, C and D do not coincide, it may be necessary to 
include additional (rigid) body degrees of freedom in the analysis 
as shown in the Figure. The corresponding additional term in the 
kinetic energy is: 
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and the following constraint conditions must also be included: 


u (B) = u(A) 

- d 

©(A) 

1 
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u .(D) = u(C) 

- d 

0(C) 

2 

2 


v (A) - V (B) 

+ d 

0 (B) 

1 

1 

1 

v (C) = v (D) 

+ d 

0 (D) 


1 2 2 

It could be argued that a similar treatment would be necessary also for the 
experimental model considered in Chapter III or for the L-100 model of 
the same Chapter. Although this is strictly true for the kinetic 
energy, such an analysis would be inconsistent with point constraint 
conditions . 

Indeed, if we consider the experimental model recalled by the 


sketch which follows: 
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we see that h<<2c h << 2c , h << 2 c . 

1 12 2 

For this structure, if it is desired to include the thickness effect, 

then care has to be taken to express the contraint conditions consistently. 

Namely, the constraint conditions between beams (1) and (2) should hold 

on the entire area F F F F and similarly the constraints between beams 

12 3 4 

m and (3) should hold at all points of the rectangle G G G G . Of 

12 3 4 

course, this induces a severe penalty on the numerical computations \ the 
latter problem can be avoided by reducing the analysis to point contraint 
conditions as indicated in Chapter III. 

Similarly, for the L-100 model recalled below: 
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where c + h and c are of the same order of magnitude, the analysis 
11 

including the cross-sectional dimensions should be applied consistent 

with constraint conditions formulated along the entire line H H . 

1 2 

Among the other numerous ways in which the aforementioned structures 
can be modeled, very interesting alternative possibilities are: 

plate representation for all components of the experimental model 
- shell representation for the fuselage of the L-100 and plate re- 
presentation for the lifting surfaces. 
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FIGURE 1 - SIMPLY SUPPORTED BEAM OF VARIABLE GEOMETRY 
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FIGURE S. - UNIFORM REAM 

CONVERGENCE OF THE ’’FINITE ELEMENT” TYPE 











FIGURE 8. - SIMPLY SUPPORTED UNIFORM BEAM 


BEAM DEFLECTION IN THE SECOND MODE 
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FIGURE 9- - STEPPED BEAM SIMPLY SUPPORTED AT BOTH ENDS 
VARIATION IN THE FUNDAMENTAL FREQUENCY WITH THE THICKNESS RATIO 











o mode shape predicted by cml method 
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FIGURE II - NON UNIFORM, CANTILEVERED HELICOPTER BLADE 
SECOND MODE OV24.6cps 
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• MODE SHAPE PREDICTED BY CML METHOD 

MODE SHAPE PREDICTED BY IVBV METHOD 



x (inches) 


FIGURE 12- NON UNIFORM, CANTILEVERED HELICOPTER BLADE 
THIRD MODE ta),«6lcps 
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FIGURE 13 • CHECK CASE OF A CANTILEVERED BEAM, 
ANTISYMMETRIC VIBRATIONS. CONVER- 
GENCE IN THE FUNDAMENTAL FREQUENCY 
VERSUS NB. OF MODES USED 
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FIGURE 14 -LIMITING CASE ^,=0 3.49 <0U,< 3.66 

FIRST BENDING MODE 
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FIGURE 19. - EXPERIMENTAL MODEL 
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FIGURE 21. - EXPERIMENTAL APPARATUS 
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FIGURE 23. - THEORETICAL RESULTS FOR THE EXPERIMENTAL MOOF-L 
u = 36.7 cps PHYSICAL DISPLACEMENTS OF THE COMPONENTS 
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FIGURE 25. - THEORETICAL RESULTS FOR THE EXPERIMENTAL MODEL 
w = 125.6 cps PHYSICAL DISPLACEMENTS OF THE COMPONENTS 
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FIGURE 26. - THEORETICAL RESULTS FOR THE EXPERIMENTAL MODEL 
w = 295.5 cps PHYSICAL DISPLACEMENTS OF THE COMPONENTS 
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FIGURE 27. - THEORETICAL RESULTS FOR THE EXPERIMENTAL MODEL 
w = 404.8 cps PHYSICAL DISPLACEMENTS OF THE COMPONENTS 
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FIGURE 28. - THEORETICAL RESULTS FOR THE EXPERIMENTAL MODEL 


oj = 536.2 cps PHYSICAL DISPLACEMENTS OF THE COMPONENTS 
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FIGURE 30. - THEORETICAL RESULTS FOR THE EXPERIMENTAL MODEL 
u = 643.3 cps PHYSICAL DISPLACEMENTS OF THE COMPONENTS 
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W IS NOT DETERMINED ACCURATELY AT THIS STEP IN THE COMPUTATION 
FIGURE 31. - THEORETICAL RESULTS FOR THE EXPERIMENTAL MODEL 


a) = 648.2 cps PHYSICAL DISPLACEMENTS OF THE COMPONENTS 



FIGURE 32. - LOCKHEED L-100 (CIVIL HERCULES) 
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FIGURE 34. - LOCKHEED L-100 HERCULES 


PHYSICAL DEFORMATIONS OF THE COMPONENTS IN THE SECOND MODE 
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FIGURE 36. - LOCKHEED L- 100 HERCULES 

PHYSICAL DEFORMATIONS OF THE COMPONENTS IN THE FOURTH WOE 
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FIGURE 37. - CONCEPT OF A PIVOTING-WING SST FOR COMMERICAL OPERATIONS 
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S) shows the model in flight with the wing at 45 deg. 
Burnett L. Gadeberg of Ames controlled the 
model. 

Variations of wing angle up to and beyond 45 
deg produced no apparent changes in stability and 
only a slight change in lateral trim — requiring a 1- 
or 2-deg offset of the ailerons. Elevator and aileron 
effectiveness remained normal and we observed no 
change in longitudinal trim. 

Ordinary maneuvers such as loops and rolls were 
performed without difficulty at wing angles of 45 
deg. Coupling between longitudinal and lateral 
motions did not appear in aileron rolls, but was 
quite apparent in the response to elevator control. 
Thus loops performed with the wing at 45 deg 
appeared to take the form of a 45-deg helix, 
indicating that the rotation produced by the 
elevator tends to align with the long axis of the 
wing. With the left tip forward, use of the elevator 
in a left turn tended to steepen the bank angle. 
Analysis shpws that in this case a certain amount of 
aileron deflection must be employed with the 
elevator to prevent banking toward the forward tip. 

Varying the sweep by turning the wing as a whole 
has several practical advantages over the usual 
“swing wing" design. It keeps the wing structure 
continuous across the pivot and makes the primary 
load on the pivot tension. With separate wing 
panels pivoted at the root, however, much greater 
loads develop on the pivots (F-6). Also, sweeping 
the wing panels back for high-speed flight 
displaces the center of lift rearward, compounding 
the normal rearward center-of-pressure shift at 


F-6 CONTINUOUS WING STRUCTURE; NO 
SENDING LOAD ON PIVOTS 



F-7 CENTROID OF LIFTING AREA NOT DISPLACED 
® Y ROTATION 



F-B STRAIGHT-BEAM STRUCTURE. 



0 VALUES OBTAINED BY THE PERTURBATION ANALYSIS 



FIGURE 38. - PIVOTING-WING VERSION OF THE LOCKHEED L-100 HERCULES 
CHANGE IN THE NATURAL FREQUENCIES VERSUS OBLICITY ANGLE a 
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TREATMENT OF THE PROBLEM BY SYMMETRY 

FIGURE 42. - CANTILIVERED PLATE IN TRANSVERSE VIBRATIONS 
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44. - POINT CLAMPED SQUARE PLATE, CONVERGENCE OF THE FUNDAMENTAL NATURAL FREOHENCY 
TOWARDS THAT OF A CANT I LIVE RED PLATE 
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FIGURE 45. - POINT CLAMPED SQUARE PLATE. 

CONVERGENCE OF THE SECOND NATURAL FREQUENCY TOWARDS THAT OF A C ANTI LIVE RED PLATE 
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FIGURE 47. - SIMPLY SUPPORTED SQUARE PLATE WITH A LINE OF INTERNAL SUPPORT 
CHANGE IN THE FUNDAMENTAL FREQUENCY WITH THE SUPPORT LENGTH. 













OF CONSTRAINT. 
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FIGURE 52. - CANTI LIVER BEAM 

CONVERGENCE OF THE INTERNAL SHEAR FORCE AT THE POINT OF CONSTRAINT. 



ABSTRACT 


The free vibrations of elastic structures of arbitrary complexity are 
analyzed in terms of their component modes. The method is based upon the 
use of the normal unconstrained modes of the components in a Rayleigh-Ritz 
analysis. The continuity conditions are enforced by means of Lagrange 
Multipliers. Examples of the structures considered are: beams with non- 

uniform properties, airplane structures with high or low aspect ratio 
lifting surface components, the oblique wing airplane as proposed by 
R.T. Jones, and plate structures. 

The method is also applied to the analysis of modal damping of 
linear elastic structures. Convergence of the method versus the number 
of modes per component and/or the number of components is discussed; the 
method is compared to more conventional approaches, to ad-hoc methods for 
some examples, and to experimental results. 



